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NOTATION 


The  body  ay  stain  of  axes,  radian  measure,  and  foot-pound-second  units  are  used 
throughout  the  paper  unless  specifically  stated  or  indicated  otherwise.  Basic  sign 
conventions  are  shown  in  Figure  1.  l.i  Section  1,  in  which  a  number  of  axis  systems 
are  considered,  the  subscripts  are  used  to  denote  quantities  referred  to  the  various 
systems  except  for  the  quantities  referred  to  the  body  system  of  axes.  The  subscripts 
for  these  quantities  are  omitted  for  convenience  except  to  ldcntifr  coordinates,  xb  , 
yb  .and  zb  . 


a,  A 


aT>  a*,  a 


a*i’  ati*  *1*1 


b 

b. B 
c 

c,  C 
C 


(Ccu)p 


perpendicular  distance  from  spring  to  knife  edge  (Fig.  13), 
ft 

polynomial  coefficients  (Section  3) 

cross-sectional  area  of  air-intake  duct  of  jet  engine  at 
entrance,  ft2 

cross-sectional  area  of  Jet -exhaust  duct  of  Jet  engine  at 
exit,  ft2 

longitudinal,  lateral,  and  normal  accelerations  of  the  air 
craft  at  the  center  of  gravity  relative  to  the  body  system 
of  axes;  positive  forward,  to  the  right,  and  up,  respec¬ 
tively,  g  units 

recorded  values  of  ax  ,  at  ,  and  an  ,  respectively; 
corrected  for  phase  lag  and  misalinement  but  not  for  loca¬ 
tion  relative  to  the  center  of  gravity,  g  units 

wing  span,  ft 

polynomial  coefficients  (Section  3) 
mean  aerodynamic  chord,  ft 
polynomial  coefficients  (Section  3) 
spring  couple  (Section  4),  ft  lb 

coefficient  of  axial  force  along  the  body  x-axis;  positive 
to  the  rear,  -X/qS 


Bcp 

phugoid  damping  coefficient,  V  + 


2C„ 


Bu  cos  acos/3 


cor  iribution  of  power  to  phugoid  damping  coefficient, 
Bc„  2C  „ 


-IV 


Bu  cos  a  cos/3 


x 


contribution  of  power  to  normal -force  coefficient 


3C, 

Ba 


CN« 


(CNa)p 
(CNa)h.  t. 


(CNa)r.  t. 


°Na 


CNq 


CNSe 

Cn-<Cn>s-<Cn>V«Vo 


Cn/3 


cn/j 


C 


nr 


variation  of  contribution  of  (C,)p  with  angle  of  attack 

variation  of  normal-force  coefficient  of  horizontal  tail 
with  local  angle  of  attack  at  the  tail;  coefficient  based 
on  horizontal -tail  area  and  local  dynamic  pressure. 


variation  of  coefficient  of  force  normal  to  vertical  tail 
with  vertical  angle  of  attack  of  vertical  tail;  coefficient 
based  on  vertical-tail  area  and  local  dynamic  pressure, 

\qr.t.Sv.  t.J 
3«v.t. 


3cn 

3  — 
2V 

3cn 

3^ 

2V 


longitudinal  phugoid  static-stability  derivative, 


3c„ 

V  ^  + 


2C„ 


3u  cos  a  cos 


3c, 

38. 


yawing -moment  coefficient  about  body,  stability,  wind,  and 
principal  z-axis,  respectively,  (yawing  moment)/qSb 


3c 

static  directional -stability  derivative, 

°P 


2V 


3c, 


n 


.  rb 

c?  — 

2V 


xiii 


g 

*1 


h 

H 


HrHrHz 


rn 


lvlrlz 

lx0‘lyo'lzo 


W  hs-  ^*8 
Ixr.Iyr.Izr 


I 


*c 


I 


i 

X 


I 


f 

z 


I 


xz 


k 


K 


acceleration  of  gravity. 


g  „  1 

=  -  sin#,  - 

V  sec 


g  1 

=  -  cos  6  sin  0,  — 

V  sec 


ft  /sec  2 


altitude,  ft 

angular  -momentum  vector  of  a  rotating  mass,  Ir|)Q,  lb  ft  sec 

angular  momentum  of  R  about  x  ,  y  ,  z  body  axes, 
respectively,  lb  ft  sec 

moment  of  inertia  of  rotating  mass  of  engine  about  its 
rotating  axis,  slug  ft2 

moments  of  inertia  of  aircraft  about  x  ,  y  ,  z  body  axes, 
respectively,  slug  ft2 

moments  of  inertia  of  aircraft  about  x  ,  y  ,  z  principal 
axes,  respectively,  slug  ft2 

moments  of  inertia  of  aircraft  about  x  ,  y  ,  z  stability 
axes,  respectively,  slug  ft2 

moments  of  inertia  of  aircraft  about  x  ,  y  ,  z  reference 
axes,  respectively,  slug  ft2 

moment  of  inertia  of  cradle  supporting  aircraft  (Section  4), 
slug  ft2 


product  of  Inertia  of  aircraft  referred  to  body  x-  and 
z-axes,  slug  ft2 

stability-augmentation-system  gain,  sec 
linear  spring  constant,  lb/ft 
correlation  constant  (Section  5) 
torsional  spring  constant,  2Kga2  ,  ft  lb/rad 


xv 


I 

L.M.N 


% 

L/3 

H 

L«a 


distance  as  defined  locally  at  tine  of  discussion,  ft 

rolling,  pitching,  and  yawing  moments  about  body  x  ,  y  , 
z  axes,  respectively,  ft  lb 

inertial  rolling,  pitching,  and  yawing  moments  about  the 
respective  body  axes,  ft  lb 


rolling,  pitching,  and  yawing  moments  due  to  gyroscopic 
action  of  rotating  mass  of  engine,  ft  lb 

rolling,  pitching,  and  yawing  moments  about  the  stability 
x  ,  y  ,  z  axes,  ft  lb 

rolling  acceleration  about  body  x-axis,  (rolling  moment)/I  , 
1/sec2 


Bl  qSb2 

—  =  Ci  -  ,  ft  lb  sec 

Bp  P  2V 

qSb2  1 
Clp  2VIX  '  sec 


1 

sec 


Bl  qSb2 

—  =  Cj_ - ,  ft  lb  sec 

Br  r  2V 


1 

sec 


L  +  ^  N 

r  r  _L 

1  2  ’  sec 

1  -  -£2- 

Vz 

Bl 

—  =  Cl/3qSb,  ft  lb 


qSb  1 


Bl 

—  =  C,$  qSb ,  ft  lb 

BS  °a 


xvi 


\ 


E*. 

L5„ 


M 

M, 


qSb  1 
*S®  Ix  ’  sec2 


Cjs^qSb,  ft  lb 

qSb  1 

Cl5r  »  2 

r  Ix  sec 

mass  of  airplane,  1/g  ,  slugs 

mass  rate  of  air  intake  of  jet  engi*:e,  slugs/sec 

mass  rate  of  jet  exhaust,  slugs/sec 

Mach  number 

ind  jated  Mach  number 


Ditching  acceleration  about  body  y-axis, 
(pitching  moment)/Iy  ,  1/sec2 


3m 


qSc2 


3q  =  2V  '  ft  lb  SCC 


=  Cffl 


qSc2  1 
"o  2VIy  '  sec 


3m  _  qSc 

— -  =  Ca,,  -  ,  lb  sec 

3u  u  V 


qSc  1 

Cmu  T  1  2 

sec' 


3m 

—  =  Cm(xqSc  ,  ft  lb 


qSc  1 
Cma  ~  ,  ^ 

Ty  sec 

3M  qSc2 

TT  =  cni(i -  •  ft  lb  sec 

di  a  2V 


qSc2 


1 


C|ni  2VIy  ’  sec 


Bm 

—  =  CraSeqSc  ,  ft  lb 


xvii 


qSc  _1_ 

Se  ly  '  sec2 

yawing  acceleration  about  body  z-axis, 
l/sec2 


3n  qSb2 

3p  "  %  2V 


ft  lb  Bee 


qf>6 
°P  2VT 


1 

Bee 


NP 

1 


1 

sec 


3n  qSb2 

37  =  C"r  lv~  *  “  lb  86C 


1 

sec 


N,  + 


^z- 


1  - 


xz 

V« 


1 

sec 


3n 

^  =  Cn/3qSb  ,  ft  lb 


1 

sec2 


3n 

—  =  Cn5r'qSb  .  ft  lb 


-ns 


qSb 

r  ~T7 


1 

sec2 


3n 

C„jaqSb  ,  ft  lb 


rcu  i 


(yawing  moment) /I, 


P.q.r 

Wr0 

PS'  V  r8 

P.q.r 

P0’  *0 
p 

Pi 

Pfpj 

PT 

p 

q 

*c 

qci 

^h.  t.  •  **v.  t. 

R 

R 

s 

sin 

S 

Sh.t.-Sv.t. 

t 


rolling,  pitching,  and  yawing  velocities,  respectively, 
about  body  axes,  rad/sec 

rolling,  pitching,  and  yawing  velocities,  respectively, 
about  principal  axes,  rad/sec 

rolling,  pitching,  and  yawing  velocities,  respectively, 
about  stability  axes,  rad/sec 

rolling,  pitching,  and  yawing  angular  accelerations, 
respectively,  about  body  axes,  rad/sec2 

rolling,  pitching,  and  yawing  accelerations,  respectively, 
about  principal  axes,  rad/sec2 

static  pressure,  lb/ft2 

indicated  static  pressure,  uncorrected  recorded  pressure, 
lb/ft2 

static  pressures  acting  across  inlet  of  air  intake  and 
exhaust,  respectively,  of  jet  engine,  lb/ft2 

stagnation  pressure,  qc  +  p  ,  lb/ft2 

period  of  short -period  oscillation,  sec 

dynamic  pressure,  t£v2  ,  lb/ft2 

impact  pressure;  dynamic  pressure  of  compressible  flow 
(Equations  (98)  and  (99)),  lb/ft2 

indicated  impact  pressure,  lb/ft2 

dynamic  pressures  at  the  horizontal  and  vertical  tail, 
respectively,  lb/ft2 

instantaneous  radius  of  turn  (Pig.  27),  ft 
reaction  force  (Section  4),  lb 

1 

Laplacian  operator,  a  +  iu; ,  - 

sec 

instrument  sensitivity 
wing  area,  ft2 

horizontal-  and  vertical -tail  areas,  respectively,  ft2 
time,  sec 

xix 


T 

thrust  due  to  power,  lb 

V* 

tine  required  for  absolute  value  of  transient  short -period 
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CONSIDERATIONS  IN  THE  DETER! I NATION  OP  STABILITY  AND  CONTROL 
DERIVATIVES  AND  DYNAHIC  CHARACTERISTICS  FRO!  FLIGHT  DATA 

Chester  H.Wolowicz 


1.  INTRODUCTION 

The  determination  of  stability  and  control  characteristics  from  flight  data  in  the 
form  of  derivatives  and  other  behavior  parameters  has  become  an  important  part  of 
flight  testing.  As  new  concepts  in  airplanes  are  developed  or  the  airplane  flies  in 
new  Mach  and  altitude  regimes,  there  is  the  need  to  verify  theory  and  wind-tunnel  data 
and  the  various  influences  on  stability  characteristics,  to  provide  information  not 
obtained  in  wind-tunnel  studies,  and  to  uncover  the  sources  of  discrepancies  between 
prediction  and  actual  flight  behavior.  Where  wind-tunnel  data  are  unavailable  or  where 
safety  of  flight  into  untested  regions  is  of  concern,  flight -determined  derivatives 
have  been  extrapolated  to  predict  airplane  behavior  prior  to  flight  into  these  regions. 

Because  of  the  exploratory  nature  of  many  of  the  investigations,  the  practical 
aspects  of  determining  derivatives  and  other  behavior  parameters,  such  as  oscillatory 
characteristics,  from  flight  data  are  very  important.  Experience  has  shown  that  a 
maximum  appreciation  and  understanding  of  the  practical  aspects  is  attained  when  back¬ 
ground  knowledge  includes  an  understanding  of  axis  systems,  transformations,  the 
equations  of  motions  and  the  limitations  of  the  equations,  techniques  used  to  determine 
the  mass  characteristics  of  the  airplane,  the  Installation  and  behavior  of  flight  test 
instrumentation,  flight  test  techniques,  and  the  theory  and  limitations  of  techniques 
used  to  determine  the  stability  and  control  characteristics  from  flight  data. 

Although  some  of  the  factors  mentioned  above,  such  as  axis  systems  and  transformations 
as  well  as  aspects  of  the  equations  of  motion,  may  be  found  in  textbooks,  the  treatment 
is  generally  not  oriented  toward  flight  testing.  Some  of  the  techniques  used  in  deter¬ 
mining  stability  and  control  characteristics  may  be  found  in  technical  reports;  however, 
limitations  of  the  techniques  occasionally  may  not  be  shown.  This  paper  attempts  to 
bring  all  the  factors  together  to  provide  a  ready  reference  of  pertinent  information. 

It  is,  in  fact,  a  greatly  expanded  version  of  AGARD  Report  224*. 

It  is  the  purpose  of  this  paper  to  discuss  the  various  factors  that  influence  the 
determination  of  stability  and  control  derivatives  and  other  behavior  characteristics 
from  flight  data.  Included  are  illustrations  of  the  application  of  flight  derivatives 
to  verification  of  predictions  and  to  determination  of  aeroelastic  effects,  stability 
criteria,  and  flight  guidance.  This  paper  is  Intended  not  only  for  the  practical 
engineer  who  is  working  with  flight  data  but  also  for  the  scientist  who  is  attempting 
to  develop  new,  sophisticated  analytical  techniques. 

*  Stability-Derivative  Determination  From  Flight  Data  hr  Chester  H.Wolowicz  and  Euclid  C.Hollem&n, 
October  1958. 
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Acknowledgment  of  Investigators  whose  work  has  directly  contributed  to  the  present 
paper  is  made  in  each  section.  It  is  recognized  that  many  noteworthy  works  of  other 
investigators  are  not  referenced. 


2.  AXIS  SYSTEMS  AND  COORDINATE  TRANSFORMATIONS 

2.1  Axis  Systems 

In  the  study  of  the  dynamics  of  the  airplane,  as  many  as  six  orthogonal  axis  systems 
may  be  used  simultaneously.  An  understanding  of  these  systems  or  reference  frames  and 
their  relation  to  the  aircraft  and  its  motions  at  various  flight  conditions  is  essential 
to  the  proper  analysis  of  flight  data.  Although  a  comprehensive  treatment  of  axis 
systems  may  be  found  in  Reference  1,  a  brief  treatment  of  the  axis  systems  is  presented 
in  this  section. 

2.1.1  Body  Systems 

The  body  axis  system  (xb,  yb,  z^  ia  body-fixed  with  its  origin  at  the  center  of 
gravity  of  the  airplane.  The  xb  axis  is  always  parallel  to  the  fuselage  reference 
line  and  when  the  center  of  gravity  is  in  the  plane  of  symmetry,  as  it  normally  is, 
both  the  xb  and  zb  axes  are  in  the  airplane' s  plane  of  symmetry,  as  shown  in 
Figure  1.  The  yb  axis  is  normal  to  the  plane  of  symmetry;  thus,  the  body  system  of 
axes  is  angularly  Invariant  with  respect  to  the  aircraft  structure. 

Because  of  its  angular  invariance  with  respect  to  the  aircraft,  the  body  axis  system 
is  an  excellent  frame  of  reference  for  mounting  flight  test  instruments.  The  orientation 
of  the  flight  test  instruments  and  their  consequent  output  relative  to  the  body  axes  - 
especially  the  linear  accelerometer  and  angular  rate  and  acceleration  sensors  -  make  it 
convenient  to  determine,  from  flight  data,  stability  and  control  parameters  with  respect 
to  this  reference  fraae.  Aside  from  convenience,  this  reference  frame  is  the  logical 
frame  about  which  to  orient  rates,  accelerations,  and  the  stability  and  control  para¬ 
meters  in  the  study  of  handling-quality  criteria,  inasmuch  as  the  orientation  of  the 
pilot  is  invariant  relative  to  this  frame. 

2.1.2  Stability  System 

Hie  stability  axis  system  (x8,  y8,  z8)  is  a  special  case  of  the  body  axis  system. 

Like  the  body  system,  the  x8  and  z8  axes  are  in  the  plane  of  symmetry  when  the 
center  of  gravity  is  in  this  plane,  and  parallel  to  the  plane  of  symmetry  when  the 
center  of  gravity  is  not  in  the  plane.  Unlike  the  body  system,  however,  the  x8  and 
z8  axes  are  angularly  variant  relative  to  the  fuselage  reference  line.  The  z8  axis 
is  perpendicular  to  the  resultant  velocity  vector  and  the  x8  axis  is  parallel  to  the 
component  of  the  resultant  velocity  vector  projected  onto  the  plane  oi  symmetry,  as 
shown  in  Figure  1. 

The  important  parametric  relationship  between  the  body  and  stability  axes  systems 
is  the  angle  of  attack,  a  ,  which  is  the  angle  between  the  x8  and  xb  axes  (Fig.  1). 

The  stability  axis  system  is  commonly  used  in  theoretical  subsonic  aerodynamics  and 
subsonic  wind-tunnel  force  and  moment  investigations.  It  is  also  employed,  on  occasion, 
in  place  of  body  axes  in  flight  test  investigations  of  longitudinal  stability  and 
control  characteristics. 
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2.1.3  Principal  System 

The  principal  axis  system  (xf,  yQ,  zQ)  defines  the  natural  axes  of  rotation  of  the 
aircraft.  They  are  the  axes  which  result  In  maximum  and  minimum  moments  of  inertia. 

The  orientation  of  this  axis  system  in  the  aircraft  is  a  function  of  the  mass  distri¬ 
bution  of  the  aircraft  and  will  remain  fixed  as  long  as  the  mass  and  mass  distribution 
remain  fixed.  When  the  lateral  distribution  of  mass  is  symmetrical  relative  to  the 
plane  of  symmetry,  which  is  generally  the  case,  the  y„  axis  will  coincide  with  the 
yb  axis,  and  the  x0  and  zQ  axes  will  lie  in  the  plane  of  symmetry,  as  shown  in 
Figure  1. 

The  inclination  of  the  xfl  axis  (Fig.  1)  to  the  x  axis  of  the  reference  axis 
system  (generally  body  axes  in  flight  test  investigations)  has  a  direct  bearing  on  the 
inertial  moments  experienced  about  the  reference  axes  as  reflected  in  the  product  of 
inertia  term  Ixz  is  the  equations  of  motion  and,  hence,  on  the  lateral  stability  of 
the  airplane. 

When  the  principal  axes  are  used  as  reference  axes,  as  they  occasionally  are  in 
theoretical  and  simulator  investigations,  they  are  used  to  simplify  the  equations  of 
motion  by  the  elimination  of  the  I__  term. 

In 

2.1.4  Wind  System 

The  wind  axis  system  is  related  to  the  resultant  velocity  vector  and  the  plane  of 
symmetry  of  the  airplane.  As  shown  in  Figure  1,  the  xw  axis  is  parallel  to  the 
resultant  velocity  vector  and  lies  in  the  transverse  plane  of  the  stability  axes 
(x8y8  plane).  Consequently,  the  zw  axis  is  coincident  with  the  z8  axis.  The  xw 
and  yw  axes  coincide  with  their  respective  counterparts  x8  and  y8  when  the 
aircraft  has  zero  sideslip. 

The  important  parameters  associated  with  the  wind  system  are  the  sideslip  angle, 

/3  ,  and  the  angle  of  climb,  y  .  By  basic  definition  the  angle  of  sideslip,  /3  ,  is 
the  angle  between  the  x„  axis  and  the  plane  of  symmetry  and  thus  lies  in  the  trans¬ 
verse  stability  axes  plane,  as  shown  in  Figure  1.  It  should  be  noted  that  not  all 
/8-sensors  necessarily  measure  this  (3  ;  this  will  be  discussed  in  Section  5  on 
"Instrumentation".  The  angle  of  climb  always  lies  in  the  vertical  plane  and  is  the 
angle  included  between  the  xw  axis  and  the  horizontal  plane. 

2  1.5  Spatial  Reference  System 

The  preceding  axis  systems  are  tied  in  with  the  plane  of  symmetry  of  the  airplane 
with  their  origins  at  the  center  of  gravity;  as  shown  in  Figure  1.  To  complete  the 
systems  of  axes  used,  at  least  one  Inertial,  space-fixed,  axis  system  is  required. 

In  dealing  with  general  motions  of  aircraft,  this  spatial  system  is  generally  earth- 
referenced  to  describe  the  motion  of  the  airplane  with  respect  to  time  for  short  time 
Intervals.  Such  a  situation  is  indicated  in  Figure  2,  which  shows  the  relationships 
of  the  various  axis  systems  previously  described  and  the  relationship  of  the  body  axis 
system  with  respect  to  the  spatial  reference  (xr,  yr,  zr).  Shown  in  the  figure  are 
flight  path  y  ,  angle  of  sideslip  /S  ,  angle  of  attack  a  ,  as  well  as  the  Euler 
orientation  angles,  </'  ,  9  ,  and  4>  of  the  airplane’ s  body  axes  relative  to  the  spatial 
axis  system.  This  is  shown  in  a  much  simpler  format  in  Figure  3.  The  sequence  of 
rotations  of  the  Euler  angles  is  important.  Generally,  the  sequence  of  rotation  is 
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\p  ,  6  ,  and  <f>  ;  this  means  that  the  airplane  is  initially  yawed,  then  pitched,  and 
finally  rolled. 

It  should  be  noted  that  y  -  9  -  a  only  when  the  aircraft  is  unbanked  (</>  =  0). 
2.1.6  Perturbation  Reference  Frmes 

In  using  perturbation  theory  in  stability  analysis,  Euler  angle  perturbations  may 
be  considered  to  be  superimposed  on  the  unperturbed  angles,  as  shown  in  Figure  4(a), 
with  the  result  that  the  perturbed  angles  are  */>  +  Ai/<  ,  9  +  A9 ,  and  4>  +  A</>  ,  or 
they  may  be  based  on  a  secondary  spatial  reference  frame  which  is  the  unperturbed 
airplane  axis  system  (xb0,  yb0.  zb0).  In  Figure  4(b)  the  unperturbed  body  axes  con¬ 
stitute  the  secondary  spatial  reference  frame  and  are  oriented  to  the  basic  spatial 
reference  frame  through  the  angles  \p  ,  9  ,  and  <f>  .  However,  the  perturbed  planes 
are  oriented  to  the  secondary  spatial  reference  plane  by  Ai//;  ,  A 9'  ,  and  A $>’  ,  which 
generally  are  not  the  same  as  Ai//  ,  A9  ,  and  Ac/>  . 


2.2  Coordinate  Transformations 

Coordinate  transformations  are  used  so  frequently  in  dynamic  studies  of  aircraft 
that  some  consideration  should  be  given  to  this  subject.  Literature  on  transformations 
is  extensive  and  ranges  from  the  classical  mathematical  treatments  (Reference  2,  for 
example)  to  engineering  applications  (References  3  and  4,  for  example).  At  this  time, 
the  most  pertinent  transformations  are  considered  to  serve  as  guidelines  for  other 
transformations  that  may  be  desired. 


2.2.1  Transformation  from  Earth  Reference 
Axes  to  Airplane  Axes 


Consider  Xr  ,  Yr  ,  and  Zr  as  generalized  vector  quantities  acting  along  the  co¬ 
ordinates  xr  ,  yr  ,  zr  ,  respectively.  Hie  transformed  vector  quantities  X  ,  Y  ,  Z 
acting  along  xb  ,  yb  ,  and  zb  axes,  respectively,  are  obtained  by  performing  three 
successive  rotations,  </;  ,  9  ,  and  c/>  ,  to  define  the  airplane’ s  orientation  with 
respect  to  the  reference  axes  xr  ,  yr  ,  and  zr  through  a  transformation  matrix  [l] 
as  follows 
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2.2.2  Transformation  from  Airplane  Axes  to  Earth  Axes 

Since  projection  from  airplane  axes  to  earth  axes  is  an  inverse  process  of  the 
preceding  transformation,  premultiplication  of  Equation  (la)  by  the  inverse  trans¬ 
formation  matrix  [l]*1  results  in 


V 

Y 
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=  [L]’1 

Y 

Zr 

Z 

(2a) 


However,  since  the  orthogonal  projections  on  the  airplane  axes  are  being  transformed 
to  orthogonal  projections  on  the  earth  axes,  the  inverse  of  the  transformation  matrix 
[l]  in  Equation  (lc)  is  the  same  as  its  transpose;  thus 
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2.2.3  Relationship  Between  Airplane  Rates  p  ,  q  , 
and  r  and  Euler  Rates  0  ,  #  ,  and  0 


(2b) 


It  should  be  recognized  from  Figure  5  that,  although  the  airplane  rate-vector 
quantities,  p  ,  q  ,  and  r  are  orthogonal,  the  Euler  rate-vector  quantities  are  not. 
Thus,  to  obtain  the  relationships  of  p  ,  q  ,  and  r  as  functions  of  0  ,  6  ,  and  0  , 
it  is  necessary  to  transform  0  ,  #  ,  and  0  to  components  along  xr  ,  yr  ,  and  zr 
axes  and  then  apply  Equation  (lc).  Die  first  transformation  is  accomplished  rapidly 
by  applying  Equation  (2b)  and  considering  each  Euler  quantity  as  a  special  case  of 
transformation  of  a  body  axis  quantity.  To  wit:  in  Equation  (2b)  both  0  and  8  are 
considered  zero  for  0  and  6  ,  and  0  is  considered  zero  for  0  .  Hence,  the  re¬ 
sulting  transformation  to  the  reference  axes  will  be 


V 

cos@cos0  -sin0  0 

0 

Yr 

= 

cos  8  sin  0  cos  0  0 

e 

_Zr 

-sin#  0  1 

0 

Substituting  Equation  (3a)  into  Equation  (lc)  results  in  the  following: 


x" 

P 

’l 

0 

-sin  8 

Y 

= 

q 

= 

0 

COS  0 

sin  0cos  8 

8 

Z 

r 

0 

-sin  0 

cos  8  COS  0 

0 

(3b) 


To  obtain  the  inverse  of  Equation  .(3b).,  it  is  necessary  to  solve  for  the  inverse  of 
the  transformation  matrix  since  0  ,  8  ,  and  0  are  not  orthogonal  and  hence  do  not 
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permit  the  use  of  the  transpose  for  the  inverse.  This  is  accomplished  by  solving  for 
the  inverse  matrix  tlJ  “ 1  in  the  relationship 


1  0  0 

[l]  [l]"1  =  0  1  0 


0  0  1 


After  solving  for  [l]"1  ,  the  inverse  of  Equation  (3b)  is  determined  to  be 


e\  =  o 


sin^tand  cos  $tan  9  p 

cos4>  -sin  4>  q 

sin<£sec#  cos<t>sec6  r 


2.2.4  Transformation  of  Euler  Angles  from  the  Body 
to  the  Stability  Axis  System 

If  two  different  rotation  series  give  the  same  starting  and  ending  orientation,  the 
matrices  representing  the  rotation  series  are  equal,  element  for  element,  in  the  two 
transformation  matrices.  Thus,  the  Euler  angles,  \pB  ,  9g  ,  and  4>B  ,  of  the  stability 
axes  can  be  derived  from  the  Euler  angles,  \ pb  ,  9b  ,  and  4^  ,  of  the  body  axes  by  the 
following  transformation  matrix  relationship 

tlJ8  =  CoJ  s  [Ll  b  (6) 

where  [l]b  is  the  transformation  matrix  of  Equation  (lc)  using  stability  axis  orienta¬ 
tion  angles  \pB  ,  #8  ,  and  <f>B  in  place  of  p  ,  6  ,  and  4>  ,  and  [l]  b  is  the  same 
transformation  matrix  using  body  axis  orientation  angles  »/»b  ,  6b  ,  and  4^  in  place 
of  4>  ,  9  ,  and  <f>  ,  if  the  same  successive  rotation  series  is  employed.  The  trans¬ 
formation  [a] 8  is  the  matrix  representing  the  transformation  from  the  body  to  the 
stability  axis  system,  or 


w8  =  0 


cos  a  0  sin  a 


-sin  a  0  cos  a 


Upon  performing  the  matrix  multiplication  shown  by  Equation  (6),  and  checking 
corresponding  elements  in  the  equated  results  to  obtain  the  most  feasible  elements  for 
the  desired  result,  the  following  relationships  are  arrived  at 


cos  a  sin  9b  -  sin  a  cos  9b  cos  4^ 

sin  4\,  cos  9b 
cos  9 . 


cos  a  cos  9b  sin  y,,  +  sin  a  (sin  ^  cos  4^  sin  9b  -  cos  \  sin  4^) 
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2.2.5  Transformation  of  Aerodynamic  Coefficients  to 
Various  Axis  Systems 

The  following  transformations  are  accomplished  readily  by  employing  Equation  (2b) 
and  replacing  </»  ,  6  ,  and  4>  in  the  equation  by  -/3  ,  a  ,  and  0  ,  respectively. 
Thus,  to  transform  from  body  to  stability  axes,  set  /3  =  0  ,  thereby  obtaining 


CD  =  Cc  cos  a  +  CN  sin  a 

«Vs  =  Cy 

CL  =  -Cc  sin  a  +  CN  cos  a 
=  Cj  cos  a  +  Cn  sin  a 
<C»>B  = 

(Cn)s  =  -Cj  sin  a  +  C„  cos  a  . 


(9) 


Similarly,  to  transform  from  body  to  wind  axes 


<c*>. 

=  -C 

c  cos  acos/3  +  Cy  sin/3  -  C, 

„  sin  a  cos  /3 

<Cy>w 

=  Cc 

cos  a  sin  p  +  Cy  cos  p  +  CN 

sin  asin/3 

<c2)w 

=  -c, 

L  =  Cc  sin  a  -  CN  cos  a 

<Cl>w 

=  c, 

cos  acos  /3  +  CB  sin  /3  +  Cn 

sin  a  cos  /3 

<Cm>w 

=  Cn 

cos/3  -  Cj  cos  asin/3  -  Cn 

sin  asin/3 

<Cn>w 

=  cn 

cos  a  -  Cj  sin  a  . 

(10) 


Also,  for  stability  to  wind  axes,  set  a  =  0  ,  obtaining 


(Cx)w  =  -CD  cos/3  +  (Cy)8  sin/3 
(Cy)w  =  CDsin/3  +  (Cy)g  cos/3 

<Cz>w  = 

(C,)w  =  (Ci)8  cos  /3  +  (Cm)8  sin/3 
<Cm>w  =  <cm)8  cos/3-  (C,)8  sin/3 
<Cn>.  =  <C„)8  • 


(il) 


To  transform  from  wind  to  stability  or  body  axes,  or  stability  to  body  axes,  use  is 
made  of  Equation  (lc). 
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2.2.6  Transformation  of  Derivatives 

The  transformation  of  derivatives  from  one  axis  system  to  another  goes  beyond  pure 
kinematic  transformations.  Longitudinal  derivatives  are  relatively  simple  in  their 
transformations;  lateral -directional  derivatives  are  more  complex  in  transformations. 
The  several  examples  will  illustrate  the  procedure  to  obtain  derivative  transformation. 
Influence  of  factors  such  as  power  is  not  considered  at  this  time. 

Transformation  of  longitudinal  derivatives  is  accomplished  by  direct  differentiation 
of  the  coefficient  equations.  This  is  possible  because  a  and  q  are  not  modified  by 
the  axis  system  used.  Fbr  example,  to  obtain  the  derivative  of  CL  with  respect  to 
a  in  a  transformation  from  body  to  stability  axes,  differentiate  the  equation  for 
CL  in  Equation  (9)  obtaining,  on  a  per  radian  basis, 

Cta  =  -CCasina  +  CNacosa- CD  .  (12) 

The  transformation  of  the  lateral-directional  derivatives  is  more  complicated,  inas¬ 
much  as  the  angular  rate  variables  r  and  p  are  affected  by  the  transformation.  At 
this  time,  sideslip  angle,  /3  ,  is  not  considered  to  be  affected  by  the  transformations 
because  of  its  definition;  however,  the  type  of  /3-sensor  used  in  flight  tests  -  whether 
it  be  a  vane,  floating  cone,  or  ball  nose  -  does  have  a  bearing  on  the  interpretation 
of  the  /3  readout  and  the  meaning  of  the  derivatives  with  respect  to  the  sensed  /3  . 
This  is  discussed  in  Section  5. 

Ooirider  the  transformation  of  lateral -directional  derivatives  from  the  stability 
to  thr  body  axis  system.  Transformation  of  the  yawing  and  rolling  moment  equations  is 
accomplished  by 


N  =  (N)„cosa  +  (L)_sina 

a  a 

L  =  (L)s  cos  a  -  (N)s  sin  a  , 
where  L  and  N  represent  rolling  and  yawing  moments,  respectively. 


(13) 


However, 


[r„b  /3b  p„b  *1 

(Cn/3)s/3  +  (Cnr)8  +  (Cn^)s—  +  (Cnp)  g  +  (Cng) sS  l  gSb 

[r  b  /3b  p_b  j 

(<V s/3  +  (Clr)s  -J-  +  (Cn/3)s  —  +  (cip)s  -^-  +  (Cls)s8|  gSb 


It  will  be  necessary  to  express  Ars  and  Apg  in  Equation  (14)  as  functions  of 
Ar  and  Ap  using  the  transform 


rg  =  r  cos  a  -  p  sin  a 


PB  =  p  cos  a  +  r  sin  a 
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\ 


Upon  substituting  Equation  (15)  into  (14)  and  Equation  (14)  into  (13),  and  regrouping 
terns. 


N 

gSb 


and 


L 

gSb 


[(Cn^s  cos  a  +  (0^)3  8^01/3  + 

+  |(Ccr)S  cos2  a  +  (Cjp)8  sin*  a  +  (Cnp  +  Cjr)3  sinacojaj  ^ 

r  1 4b 

+  KCn>g)8  COSO  +  (C IjpB  Sin  Oil  —  + 

+  J^Cnp)8  cos2  a  -  (Cjr)8  sin2  a  -  (Cnr  -  Cip)8  sin  acos  aj  ^ 
+  |(Cnj)s  cos  a  +  (Clj)ssinaJs 

[< 


pb 

2V 


(Cl^)a  cos  a  -  (Cn;3)sSin  aJ/3  + 

+  cos2  a  +  (Cnr)s  sin2 a  -  (CDp  +  Cir)8  sinacos  aj 

+  J^Cj^gis  cos  a  -  (Cnp  a  sin  cTJ  ~  + 

+  [<Clr)8  cos2  a  -  (Cnp)8  sin2  a  -  (C„r  -  Cjp)8  sin  acos  aj  — 


( 


+  (Clj)s  cos  a  -  (Cns)s  sin 


a  8  . 


(16a) 


(16b) 


Summaries  of  transformations  of  aerodynamic  derivatives  from  stability  to  body  axis 
system,  and  vice  versa,  are  given  in  Tables  I  and  II. 


2.2.7  Transformation  of  Moments  of  Inertia  from 
One  Axis  System  to  Another 

Although  this  topic  iB  covered  in  applied  mechanics  literature,  an  illustrative 
example  is  given  as  a  refresher.  Also  Included  are  tables  of  transformations  for 
ready  reference. 

To  obtain  IXs  in  terms  of  body  axes  quantities,  use  is  made  of  the  fundamental 
relation 


^8  =  *  ■ 


(17) 
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Substituting  the  following  transform 


into  Equation  (17), 
x  cos  a  +  z  sin  a 
z  cos  a  -  x  sin  a 


(18) 


and  expanding, 

1  xg  =  [/(yb  +  zb)  <*“]  cos201-*-  [/(yb  +  xb>  <*“]  sin2  a  -  2^xbzbdmj  sinacosa 

=  Ix  cos2  a  +  I2  sin2  a  -  2  Ixz  sin  acos  a 
=  t  (Iz  +  Ix)  -  i  (I2  -  Ix)  cos  a  -  Ixz  sin  2  a  . 


3.  EQUATIONS  OF  MOTION 

The  equations  of  motion  of  an  airplane  as  found  in  texts  on  aircraft  dynamics  (such 
as  Reference  5)  and  as  normally  presented  in  the  technical  literature,  ilthough  prosaic 
in  appearance,  do  contain  complexities  in  the  significance  of  the  individual  terms. 

The  following  discussion  is  intended  to  acquaint  the  reader  with  the  scope  of  the 
complexities  which  may  be  encountered  and  which  should  be  recognized  and  managed  in 
dealing  with  the  equations  of  motion.  An  understanding  of  this  matter  is  important 
in  applying  the  equations  to  derivative  determination  from  flight  *ata. 

3. 1  Inertial  Quantities 

In  all  considerations  of  the  inertial  portions  of  the  equr  >ns  of  motion,  the  axis 
system  used  has  a  direct  bearing  on  the  expressions  for  inertial  forces;  the  degree 
of  asymmetry  of  the  mass  distribution  of  the  aircraft  and  the  magnitude  and  violence 
of  the  aircraft  motl'V'g  affect  the  format  of  the  expressions  for  inertial  moments. 

It  in  assumed,  for  the  purposes  of  this  paper,  that  the  aircraft  behaves  as  a  rigid 
body.  Where  aeroelasticity  is  a  factor,  it  is  assumed  that  proper  precautions  will 
have  been  taken  to  provide  assurance  that  the  rigid-body  concept  will  provide  a  good 
degree  of  approximation 

Inertial  quantities  arise  from  the  inherent  action  of  the  aircraft  whose  various 
components  act  as  a  rigid-body  assembly  and  from  the  rotating  masses  attached  to  the 
aircraft. 

3.1.1  Inherent  Aircraft  General  Inertial  Force  Expressions 

Inasmuch  as  our  interest  lies  in  the  analysis  of  flight  test  data  oriented  to  the 
body  axis  system,  the  inert) al  force  expressions  applicable  to  this  axis  system  and 
for  all  attitudes  of  flight  are 


Yi 


m(u  +  qw  -  rv) 

(19a) 

m(v  +  ru  -  pw) 

09b) 

m(w  -  qu  +  pv)  . 

(19c) 
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If  the  stability  axis  system  were  employed  as  the  reference  instead  of  the  body  axis 
system,  qw  -  pw  =  w  =  0  ,  inasmuch  as  there  is  no  linear  velocity  component  along  the 
z-stability  axis, 

3.1.2  General  Inertial-Moment  Expressions 

For  the  general  case  where  the  principal  axes  are  asymmetric  to  the  various  planes 
of  the  reference  axes,  the  inertial -moment  expressions  are 

Lt  -  Ixp  +  I^irp  -  q)  -  Ixz(r  +  pq)  +  Iyz(r2  -  q2)  +  (Iz  -  Iy)qr  (20a) 

=  Iyq  +  Iyz(rq  -  r)  -  I^r  +  qr)  +  Ixz(p2  -•  r2)  +  dx  -  Iz)rp  (20b) 

Nj  =  Izr  +  Ixz(qr  -  p)  -  Iyz(q  +  rp)  +  I^q2  -  p2)  +  (Iy  -  Ix)pq  .  (20c) 

Fortunately,  situations  involving  general  asymmetry  of  the  aircraft  are  rare.  Normally, 
the  vehicle  will  have  a  mass  distribution  symmetrical  relative  to  the  xz-body  plane  of 
symmetry,  with  the  result  that  the  principal  y-axis  coincides  with  the  y-body  axis. 
Under  such  circumstances,  1^  =  Iyz  =  0  and  the  general  inertial -moment  expressions 
reduce  to  the  following  normally  employed  form: 


Li 

=  hi 

-  Ixz(r  +  pq)  +  (Iz  -  Iy)qr 

(21a) 

Mi 

=  IyQ 

+  Ixz<P2  -  r2>  +  dx  -  Iz)rP 

(21b) 

Ni 

=  lZi 

-  Ixz(p  -  qr)  +  (Iy  -  Ix)  pq  . 

(21c) 

The  inertial  expressions  in  Equations  (21a,  b,  c)  are  nonlinear  and  thus  not  suitable 
for  use  in  the  derivation  of  closed-form  stability  equations.  However,  they  are  re¬ 
quired  in  analog  or  digital  computer  study  of  the  motion  of  the  aircraft  in  general 
or  violent  maneuvers  and  in  the  analog  matching  of  flight  data  from  such  maneuvers 
in  attempts  to  determine  the  effective  values  of  the  stability  and  control  derivatives 
for  the  maneuver. 

In  violent  maneuvers,  the  terms  involving  pq  and  rp  are  particularly  important. 
These  terms,  as  well  as  qr  ,  are  gyroscopic  terms.  Modern  high-performance  aircraft 
tend  to  have  low  values  of  Ix  compared  to  Iy  and  Iz  ,  with  the  result  that  gyro¬ 
scopic  action  represented  by  (Ix  -  Iz)rp  and  (Iy  -  Ix)pq  ,  in  particular,  has  been 
responsible  for  the  uncontrollable,  catastrophic  roll -coupling  behavior  of  at  least 
one  jet  aircraft  after  a  deliberate  high  roll  rate  input. 

When  the  motions  of  the  aircraft  are  small  or  gradual,  the  inertial -moment  expressions 
may  be  simplified  to 

=  lxp  -  Ixzr  (22a) 

=  Iyq  (22b) 

Ni  =  V  "  JxzP  • 


(22c) 
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3.1.3  Small-Perturbation  Inertial  Expressions 

The  classical  approach  to  the  study  of  aircraft  dynamic  stability  and  control  involves 
the  use  of  small  disturbances  (perturbations).  Restricting  the  motions  to  small  de¬ 
viations  from  steady -state  conditions  allows  the  elimination  of  non-linear  terms  from 
the  inertial  expressions.  Such  motions  are  useful  in  defining  the  stability,  control, 
and  handling  qualities  of  the  aircraft,  and  the  pilot  effort  or  autopilot  character¬ 
istics  required  to  control  the  motion.  It  has  been  found  that  the  use  of  small- 
perturbation  theory  gives  good  results  and  permits  the  development  of  analytical 
expressions. 

To  arrive  at  the  small -perturbation  inertial  expressions,  replace  the  individual 


acceleration  and  velocity  terms  in  Equations  (19a,  b,  c)  and  (22a,  b,  c)  by  accelera¬ 
tions  and  velocities  made  up  of  disturbances  superimposed  on  equilibrium  conditions 
so  that  u  .  etc.,  is  replaced  by  u  +  Au  ,  etc.,  respectively;  expand  the  product 
terms;  neglect  the  second-order  quantities  (ArAu,  for  example);  and  subtract  the  initial 
conditions  from  the  final  resulting  conditions.  The  resulting  small-perturbation 
inertial  expressions  are 

AXi 

=  m!Au  +  wAq  +  qAw  -  rAv  -  vAr] 

(23a) 

AYi 

=  mlAv  +  uAr  +  rAu  -  wAp  -  pAw] 

(23b) 

•H 

NJ 

< 

=  m[Aw  -  uAq  -  qAu  +  pAv  +  vAp] 

(23c) 

and 

ALi  =  -  XxzAf 

(24a) 

Ama  =  IyAq 

(24b) 

ANi  =  . 

(24c) 

Equations  (23a,  b,  c)  show  that  lateral-directional-mode  perturbations  Av  ,  Ar  , 
and  Ap  appear  in  the  longitudinal -mode  equations  Axa  and  Azt  ,  and  that  the 
longitudinal -mode  perturbations  Au  and  Aw  appear  in  the  lateral -directional -mode 
equation  Azt  .  This  coupling  of  the  two  modes  can  normally  be  minimized  to  permit 
practical  use  of  the  uncoupled  practical  approximation  of  Equations  (23a,  b,  c)  shown 
in  Equations  (25a,  b,  c).  This  minimization  is  achieved  in  flight  test  maneuvers  such 
as  elevator  pulses  for  perturbation  of  the  longitudinal  mode  and  rudder  or  aileron 
pulses  for  perturbation  of  the  lateral -directional  mode  initiated  during  steady  wings- 
ievel  or  steady  turn  flight. 


AXj  =  ml^u  +  wAq]  (25a) 

Ay1  =  m[Av  +  uAr  -  wAp]  (25b) 

AZj  =  m&^w  -  uAq  -  qAu]  .  (25c) 
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3.2  Gyroscopic  Couples  of  Rotating  Masses 

Spinning  masses  mounted  on  the  aircraft  -  such  as  propellers  and  rotating  elements 
of  engines  -  possess  angular  momentum  relative  to  the  reference  (body)  axes  and  pro¬ 
duce  gyroscopic  couples  on  the  aircraft  which  could  be  significant,  as  was  the  case  on 
the  X-5  airplane6.  Normally,  the  gyroscopic  couples  are  negligible;  however,  the  advent 
of  vertical -rising  aircraft  with  tilting  engines  and  the  increase  in  size  of  propulsion 
units  on  high-performance  aircraft  make  it  inadvisable  to  arbitrarily  ignore  this 
coupling. 

For  a  rotating  mass  having  a  rotating  axis  in  or  parallel  to  the  xz-plane  of  symmetry 
but  at  an  angle  #rm  to  the  x-body  axis  (Fig.  6),  it  can  be  shown  from  the  moment  of 
momentum  relation,  co  x  fl  and  H  =  IrmQ  ,  that  the  gyroscopic  couple  about  each  of  the 
body  axes  is 


Lrm  =  «Hz  “  rHy  =  -1 rm ^  sin  6rm 

(26a) 

Mrm  =  rHi  "  PHZ  =  Irmn(rcos0rni  +  P sin 

(26b) 

Nrm  =  PHy  "  <1HX  =  cos  . 

(26c) 

These  rotating  mass  contributions  are  added  to  the  inertial  moments 
Equations  (20a,  b,  c),  (21a,  b,  c),  and  (22a,  b,  c). 

expressed  by 

For  small  perturbations  of  the  aircraft,  the  perturbations  of  tlr. 
resulting  from  the  rotating  mass  are  expressed  by 

.  gyroscopic  couples 

ALrm  =  -IrmnAc»sin0riii 

(27a) 

AMrm  =  Irmn  <Ar  008  9ra  +  AP  sin  0rm> 

(27b) 

ANrm  =  -IrmnAtlcoserm  • 

(27c) 

These  perturbations  are  added  to  Equations  (24a,  b,  c)  when  significant,  in  which 

case,  Equations  (24a,  b,  c)  will  become  inter -dependent  because  of  the  coupling  of  the 

longitudinal -mode  and  lateral -directional -mode  moment  equations.  It  should  be  noted 

that,  if  6rm  were  variable,  the  above  relations  in  Equations  (27a,  b,  c)  would  have 

required  further  expansion  and  introduced  an  additional  degree  of  freedom  in  the  form 

of  A#  . 
r  m 

3.  3  Gravitational  Force 

The  gravity  force  will  not  contribute  to  the  moment  equations  as  long  as  the  origin 
of  the  axis  system  is  at  the  centre  of  gravity. 
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3.3,1  Components  of  Gravity  Force 

With  the  gravity  force  R  acting  along  the  zr  axis,  the  expressions  for  the 
components  of  gravity  force  acting  along  the  body  axes  are  readily  deduced  from  Figure 
4(a)  to  be 


-W  sin  9 

(28a) 

W  cos#  sin<£ 

(28b) 

Rcos^cos0  . 

(28c) 

These  components  are  subtracted  from  the  inertial -force  equations  (19a,  b,  c). 


3.3.2  Small  Perturbations 

Small  perturbations  of  the  components  of  the  gravity  force  may  be  based  on  Euler 
angle  perturbations  superimposed  on  the  unperturbed  angles  using  the  same  basic 
reference  frame,  or  on  Euler  angle  perturbations  relative  to  a  secondary  spatial 
reference  frame  made  up  of  the  unperturbed  aircraft  axis  system  as  discussed  in 
Section  2.1.6.  and  as  shown  in  Figure  4(b).  In  this  second  approach,  unperturbed 
body  axes  are  used  as  the  secondary  spatial  reference  when  interest  is  primarily  in 
perturbations  of  body-oriented  flight  test  data. 

Using  the  first  approach,  replace  9  and  <t>  in  Equations  (28a,  b,  c)  by  9  +  A# 
and  <p  +  A 4>  ,  respectively;  expand  the  resulting  trigonometric  functions,  consider 
cosA_~  l  ,  sin  A  A_ ,  and  A_A_~  0  ;  and  subtract  the  initial  conditions  from 

the  result.  The  resulting  small -perturbation  expressions  are 


AX* 

=  -wA$  cos  9 

(29a) 

AYg 

-  W(A#  cos  9  cos  4>  -  At?  sin  9  sin<£) 

(29b) 

AZg 

=  -W(A<£  cos  0 sin 0  +  A#  sin#cos<£)  . 

(29c) 

In  the  second  approach,  using  the  unperturbed  body  axes  as  the  reference  and 
Ai p'  ,  A 9'  ,  and  A$'  (Fig.  4(b))  as  the  Euler  angles  of  the  perturbations,  the  per¬ 
turbations  of  the  components  of  gravity  are  obtained  by  using  Equations  (lc)  and 
(28a,  b,  c).  In  Equation  (lc)  the  generalized  quantities  Xr  ,  Yr  ,  and  Zr  are 
replaced  by  the  expressions  for  Xg  ,  Yg  ,  and  Zg  ,  respectively,  as  given  in 
Equations  (28a,  b,  c);  and  the  Euler  angles  \p  ,  9  ,  and  4>  are  replaced  by  Ai/>' , 

A 9'  ,  and  A ,  respectively.  The  generalized  quantities  Xb  ,  Yb  ,  and  Zb  in 
Equation  (lc)  are  now  equal  to.  (Xg  +  Axg),  (Yg  +  Ayg),  and  (Zg  +  Azg),  respectively. 

By  subtracting  the  initial  conditions  (Equations  (28a,  b,  c))  from  the  resulting 
perturbed  equation  after  considering  cosA_~  i  ,  sinA_^A_  ,  and  A_A_Oi  0  , 
the  perturbation  expressions  for  this  second  approach  will  have  the  following  form: 

Axg  =  H  (cos^  sin0  Av/» '  -  cos #cos0  A#')  (30a) 

AYg  =  W  (sinAv/>'  +  cos  9cos<t>  A<£')  (30b) 

Azg  =  -W(sinf9A#'  +  cos  <9sin0  A<£')  . 


(30c) 
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Hie  advantage  in  using  Equations  (30a,  b,  c)  instead  of  Equations  (29a,  b,  c)  is 
that,  for  small  perturbations  during  highly  banked  as  well  as  wings-level  flight, 


bxp'  ~  jAr  dt 

Ola) 

A 6'  ~  jAq  dt 

(31b) 

A </>  ~  jAp  dt  . 

(31c) 

To  apply  such  simple  integrations  to  A</»  ,  A0  ,  and  A<£  in  Equations  (29a,  b,  c) 
requires  that  <p  and  6  be  small. 

Both  Equations  (29a,  b,  c)  and  (30a,  b,  c)  show  coupling  of  the  longitudinal  and 
lateral  modes.  In  both  sets  of  equations,  the  longitudinal  modes  (Axg  and  Azg) 
are  uncoupled  from  the  lateral -mode  perturbations  by  performing  a  longitudinal  pulse 
when  initial  conditions  are  steady-state.  In  performing  a  lateral -directional  pulse 
from  steady-state  conditions,  the  lateral -mode  expression  (30b)  is  inherently  un¬ 
coupled  from  longitudinal  perturbations,  whereas  expression  (29b)  shows  interaction 
of  the  longitudinal  perturbation  A0  which  is  excited  by  the  lateral -directional 
pulse. 

When  banked  and  climbing  fl ight  are  being  considered,  it  may  be  surmised  from  the 
preceding  that  Equations  (30a,  b,  c)  are  more  amenable  than  Equations  (29a,  b,  c)  to 
theoretical  stability  analysis  and  for  analysis  of  flight  data  when  longitudinal  or 
lateral  pulses  are  applied  from  initial  steady-state  conditions. 

3.4  Aerodynamic  Derivatives 

In  stability  and  control  investigations  based  on  flight  data,  the  previously  dis¬ 
cussed  inertial,  gyroscopic,  and  gravitational  quantities  are  normally  equated  to 
aerodynamic  parameters  only.  This  is  done  primarily  to  facilitate  the  analysis  of 
flight  data.  However,  in  doing  this,  the  parameters  are  no  longer  pure  aerodynamic 
parameters,  inasmuch  as  they  will  have  been  modified  by  influences  arising  from  power 
and  aeroelasticity  as  well  as  possible  other  sources.  Generally,  these  influences 
can  be  accounted  for  and  the  pure  aerodynamic  parameter  arrived  at. 

Inasmuch  as  the  equations  are  set  up  under  the  principle  of  super-position  of 
influences,  situations  may  be  encountered  in  which  the  accuracy  of  the  results  obtained 
from  the  equations  will  deteriorate.  This  is  of  particular  concern  where  very  rapid 
control  inputs  are  encountered.  Also,  inasmuch  as  the  aerodynamic  parameters  are  in 
the  form  of  derivatives,  care  must  be  exercised  not  to  exceed  the  validity  of  the 
derivative. 

Finally,  there  are  some  limitations  in  combining  several  of  the  derivatives, 

Cnr  -  Cn^g  ,  for  example. 

Consideration  is  given  at  this  time  to  the  above-mentioned  factors  which  have 
significance  in  the  utilization  of  aerodynamic  derivatives  in  the  equations  of  motion 
and  in  the  determination  of  the  derivatives  from  flight  data.  For  convenience,  the 
conventional  derivatives  are  tabulated  overleaf. 
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Longitudinal  Derivatives 


Lateral-Directional  Derivatives 


3.4.1  Significance  of  the  Derivatives 

The  aerodynamic  derivative  provides  the  slope  of  the  curve  of  the  aerodynamic  force 
or  moment  coefficient,  as  the  case  may  be,  with  respect  to  an  independent  variable  - 
other  independent  variables  being  considered  constant  -  at  a  particular  value  of  the 
variable.  In  analog  simulation  studies,  nonlinear  curves  are  reduced  to  straight-line 
segments,  each  segment  being  valid  only  for  an  incremental  range  of  the  independent 
variable. 

In  the  inverse  problem  of  obtaining  derivatives  from  flight  data,  the  derivative  is 
valid  only  for  the  incremental  range  of  disturbance  oi  the  independent  variable,  at  the 
steady-state  condition,  used  in  determining  the  derivative.  An  example  involving  a 
nonlinear  variation  of  Cn  with  /3  is  shown  in  Figure  7.  In  this  example,  the 
origin,  0,  represents  steady  state  and  (A/3)t  and  (A/3)z  represent  two  disturbance 
ranges  of  the  variable.  It  will  be  noticed  that  the  derivative  obtained  may  differ 
appreciably  in  magnitude  because  of  the  nonlinearity  of  the  curve  in  the  disturbance 
ranges  and  2  . 
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?. 4.  2  Unsteady  Flow  Effects 

In  dealing  with  the  derivative  concept  in  accounting  for  the  influence  of  independent 
variables  on  an  aerodynamic  force  or  moment  coefficient,  for  example 

,  „  „  Arb  A5b  Apb 

Acn  =  ^A^  +  C„r  —  +  C„j  +  Cnp  —  +  cnsrASr  , 

it  is  assumed  that  each  derivative  contributes  to  the  total  as  though  it  acted  alone 
and  that  the  aerodynamic  force  and  moment  coefficients  are  functions  of  the  instan¬ 
taneous  values  of  the  disturbance  displacements  and  velocities,  control  angles,  and 
their  derivatives.  Further,  the  derivatives  are  based  on  the  variation  of  the  co¬ 
efficients  under  near -steady -state  conditions  of  the  variable.  Although  the  deriva¬ 
tive  concept  of  treating  aerodynamic  force  and  moment  perturbations  has  generally 
worked  well,  the  application  to  situations  of  rapidly  changing  independent  variables 
(unsteady  flow  conditions),  as  in  the  case  of  a  very  rapid  control  displacement  or  a 
sharp-edged  gust,  does  not  necessarily  give  correct  answers.  This  is  due  to  apparent 
mass  effects  of  the  air,  whose  inertia  will  not  produce  instantaneous  changes  in 
circulation  and  consequently  causes  aerodynamic  lag.  This  is  illustrated  in  Figure  8, 
which  shows  the  variation  of  CN  as  a  function  of  nondimensional  time,  c/2V  ,  as  a 
result  of  a  step  gust.  The  derivative  concept  would  show  a  constant  slope  curve, 
whereas  the  actual  variation  of  CN(t)  would  show  a  lag  at  the  initial  instance  of 
the  step  gust  input. 

When  an  aircraft  is  oscillating  sinusoidally,  the  lift  will  follow  the  sinusoidal 
variation  in  angle  of  attack  but  will  be  of  smaller  magnitude  and  there  will  be  a 
phase  difference  between  the  lift  and  angle  of  attack.  This  unsteady  flow  effect  is 
a  function  of  reduced  oscillating  frequency,  wc/2V  ,  as  well  as  Mach  number.  Although 
the  magnitude  of  Cna  is  not  normally  affected  appreciably  for  normal  airplane 
osci’liting  frequency  conditions,  the  phase  lag  may  bring  about  a  large  change  in  Cn^  . 
This  may  be  of  considerable  importance  in  pitch  damping  of  tailless  aircraft  (Ref. 7). 

In  general,  all  the  aerodynamic  derivatives  behave  in  a  similar  manner.  Thus,  it 
is  seen  that  attempts  to  use  the  derivative  concept  in  analog  simulators  involving 
very  rapid  changes  of  the  independent  variables  can  lead  to  errors;  conversely,  deter¬ 
mination  of  derivatives  from  flight  data  requires  awareness  of  the  maneuvering  or 
unsteady  flow  factors  mentioned  which  can  influence  the  magnitude  of  the  derivative. 


3,4.3  Derivatives  with  Respect  to  u 

Aerodynamic  derivatives  with  respect  to  u  are  of  concern  when  phugoid  modes  are 
being  investigated.  Because  this  mode  is  often  overlooked,  these  derivatives  are 
generally  unfamiliar.  Thus,  some  consideration  is  given  to  them  at  this  time  for 
future  reference  as  needed.  Consider  -Z  =  CNqS  .  Differentiation  with  respect  to 
u  shows 


Bz 

Bu 


BcN  BV 

3?  ®  + 


(32) 
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where  3v/c>u  -  l/(cos  acos/?)  from  u  =  Vcos/3cosa  .  The  Az  due  to  change  In  u 
is  now  expressed  as 


(33) 


Since 


is  more  than  simply  the  variation  of  the  normal  coefficient  with  respect  to  velocity, 
u  .  it  can  fittingly  be  called  the  effective  aerodynamic  derivative  of  Cj,  with  respect 
to  u  ,  or  Cnu  .  Similarly,  the  effective  aerodynamic  derivatives  of  CB  and  Cc 
with  respect  to  u  are 


cos  a  cos  fi ) 


and 


or  Chu  and  Pc„  .  respectively. 


d.U.U  Derivatives  with  Respect  to  q  and  a  , 
and  r  and  $ 

It  is  customary  in  reporting  flight-determined  derivatives,  wherein  transient 
oscillations  are  used  in  the  analysis,  to  pair  the  derivatives  varying  with  respect 
to  q  and  &  and  those  varying  with  respect  to  r  and  $  .  For  example 


Aqc 

Adc 

Aqc 

Cmq_2V 

+  cm  ■  ~ 

“  2V 

(C^ 

+  Cn^) 

2V 

and 

Arb 

A/3b 

Arb 

cnr - 

r  2V 

+  Cn^l 7  ~ 

(Cnr 

-  CB/g) 

2V 

Inasmuch  as  the  phenomenon  involving  d  is  different  from  that  involving  q  ,  and 
the  phenomenon  involving  r  is  different  from  that  of  $  ,  the  pairing  is  valid  only 
when  small-perturbation  transient  oscillations  of  a  maneuver  are  involved  and  satisfy 
the  linearized  equations  of  motion.  In  addition,  although  the  pairing  works  well  for 
the  longitudinal  equations  whether  or  not  stability  or  body  axes  are  employed,  the 
validity  of  the  paring  for  the  lateral-directional  equations  is  dependent  on  the  use 
of  the  stability  system  of  axes;  if  body  axes  are  used,  the  pairing  of  r  and  $ 
derivatives  is  permissible  at  low  angles  of  attack. 


In  performing  a  small-perturbation  longitudinal  transient  oscillation,  the  center 
of  gravity  of  the  aircraft  tends  to  move  along  the  flight  path  as  though  it  were  not 
disturbed;  consequently,  the  amplitude  ratio  |Aql/lAdl  is  similar  to  1.0  and  the 
vector  quantities  Aq  end  Ad  are  approximately  in  phase.  Thus  Aq  can  be  sub¬ 
stituted  for  Ad  .  In  the  case  of  a  lateral-directional  (Dutch  roll)  transient 
oscillation  relative  to  the  stability  axis  system,  the  aircraft,  in  tending  to 
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maintain  its  center  of  gravity  along  the  flight  path  as  though  it  were  not  disturbed, 
will  experience  Ar  c^-A/3  ,  inasmuch  as  r  and  /3  are  now  referred  to  the  sane  axis 
system.  Thus,  the  amplitude  ratio  |Arl/!A/?t  is  similar  to  1.0,  but  the  phase  relation 
is  approximately  180°;  consequently,  the  sign  of  the  A/3  derivative  is  changed  to  minus 
in  pairing  r  and  /3  derivatives.  In  dealing  with  the  body  axis  system,  |Ar|/|A/3| 
and  $1^3  can  differ  appreciably  from  1.0  und  lo0°,  respectively,  at  high  angles  of 
attack. 

It  is  reiterated  that  pairing  the  derivatives  is  valid  only  for  the  special  con¬ 
ditions  mentioned.  On  the  other  hand,  it  has  not  been  possible  to  solve  for  the  & 
derivatives  independent  of  the  q  derivatives,  and  /3  derivatives  independent  of  r 
derivatives,  from  flight  data  with  any  degree  of  consistency  and  confidence. 

3.4.5  Power  Effects 

The  propulsive  system  may  have  a  significant  influence  on  the  stability  as  well  as 
the  trim  of  the  airplane.  Its  force  and  moment  contributions  to  the  equations  of 
motion  may  be  presented  as  derivatives  in  the  equations.  If  the  power  contributions 
are  not  accounted  for  by  their  own  derivatives,  they  will  be  reflected  in  the  magnitudes 
of  the  aerodynamic  stability  and  control  derivatives  which  will  then  become,  in  essence, 
effective  derivatives.  A  comprehensive  treatment  of  power  effects  is  complex  and  beyond 
the  scope  of  this  paper.  Only  major  effects  are  considered,  to  show  how  propulsion 
system  derivatives  contribute  to  the  effective  values  of  the  aerodynamic  derivatives. 

It  is  essential  at  this  time  to  emphasize  an  important  point  regarding  consideration 
of  the  effect  of  power  on  stability.  True  inherent  stability  of  the  aircraft  with 
power  on  can  only  be  evaluated  by  keeping  the  settings  of  the  engine  and  propeller 
controls  fixed  during  the  maneuver.  Any  maneuver  that  entails  alteration  of  the  pro¬ 
pulsive  system  controls  during  the  maneuver  will  not  provide  a  true  index  of  the 
stability  from  an  analysis  of  the  time  history  of  the  maneuver. 

Influence  of  propellers :  Influences  of  propellers  consist  of  direct  propeller  effects 
and  also  indirect  effects  due  to  the  propeller  slipstream  on  the  wing-fuselage  and  the 
tail  surfaces. 

Direct  propeller  effects:  Direct  propeller  effects,  as  shown  in  Figure  8,  consist 
of  a  direct  thrust  T  acting  along  the  thrust  axis,  and  a  transverse  forve  (Y)p  , 
as  well  as  a  normal  force  (-Z)p  ,  perpendicular  to  the  thrust  axis  in  the  plane  of 
the  propeller  disk.  The  thrust  T  is  a  primary  function  of  a.  and  V  .  Quantitative 
determination  of  the  normal  and  transverse  forces  (-Z)p  and  (Y)p  may  be  accomplished 
by  solving  for  (CNa)p  and  (Cy^)p  ,  as  discussed  in  References  8  and  9.  Actually, 
the  derivatives  are  of  more  concern  for  the  purposes  of  this  paper  than  the  actual 
magnitudes  of  the  forces. 

The  contributions  of  the  direct  propeller  effect  (Fig. 9)  on  longitudinal  and  lateral 
stability  are  reflected  in 

(Craa)P  =  CTa^+(CNa)p^  (35) 

c  c 

and 

=  ~(cyp)p  "S'  • 


(36) 
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It  is  opportune  to  note  the  influences  of  the  direct  propeller  effects  on  CLa 
when  performing  a  transformation  from  the  body  to  the  stability  axis  system.  The  net 
effective  CL  and  CD  of  the  aircraft,  in  the  absence  of  angular  rates  and  for  fixed 
controls,  can  be  expressed  as 


Aero  Power 


Cj 

t 

cos  a 

-sin  a 

sin  otp 

COS  Op 

_CD 

sin  a 

cos  a 

-  cos  ap 

sin  ap 

(37) 


where  the  direct  thrust  is  assumed  to  be  vectored  parallel  to  the  body  x-axis  and 
<xp  =  ex  +  ep  ,  Differentiating  Equation  (37)  with  respect  to  a  for  Cija 


CLa  =  (Cnacosa-  CC(xsina)  -  Cp  +  CTasinap  +  (CNa)pcosa  ,  (38) 


where  CD  is  the  effective  value  as  shown  in  Equation  (37). 


Study  of  Equation  (38)  shows  that  power  increases  the  effective  CLq,  of  the  aircraft. 
On  low-performance  aircraft,  the  power  effect  is  generally  negligible. 


Propeller  slipstream:  The  propeller  slipstream  influences  the  distribution  of  the 
aerodynamic  forces  on  the  aircraft  structure  as  a  result  of  (i)  the  increase  in  local 
velocity  over  the  structure  due  to  and  in  the  propeller  slipstream,  and  (ii)  upwash  and 
downwash  effects  of  the  rotating  slipstream  of  the  propeller.  The  slipstream  can  be 
stabilizing  or  destabilizing,  depending  upon  the  direction  of  rotation  of  the  propeller 
and  the  position  of  the  tail  relative  to  the  rotating  slipstream.  Analytical  techniques 
to  quantitatively  account  for  the  propeller  slipstream  effects  on  the  stability  of  the 
aircraft  have  not  been  satisfactory.  Generally,  powered  models  are  used  to  provide 
engineering  data  on  new  designs. 


Influence  of  jet  enpines:  The  jet  engine  has  the  counterpart  of  the  effects  that 
were  shown  for  the  propeller.  It  provides  a  direct  thrust,  shows  normal  and  transverse 
force  effects  at  the  entrance  of  the  intake  duct,  and  -  depending  on  geometry  -  is 
capable  of  influencing  the  equilibrium  and  stability  of  the  aircraft  by  inflow  of  air 
into  the  jet  exhaust.  Unlike  the  propeller,  the  influence  of  the  jet  engine  on  the 
tail  surfaces,  and  hence  the  stability  of  the  airplane,  is  amenable  to  analytical 
techniques  to  quantitatively  account  for  these  effects. 


The  thrust  produced  on  the  aircraft  equipped  with  a  jet  engine  is  equal,  as  shown 
in  Figure  10,  to  the  vectorial  change  in  momentum  of  the  air  and  fuel  passing  through 
the  engine  plus  the  resultant  of  the  pressure  forces  acting  across  the  inlet  and  outlet 
areas.  Where  the  intake  and  exhaust  are  in  line  with  the  thrust  axis  and  the  x-body 
axis 


T  =  CrqS 


mjVj  -  maV  cos  Op  +  (pjAj  -  Pj Aj )  . 


(39) 
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A  change  in  direction  of  the  momentum  vectors  at  intake  or  exhaust  relative  to  the 
x-body  axis  brings  into  being  forces  normal  to  the  body  x-axis.  Where  the  jet  exhaust 
is  parallel  tc  the  x-body  axis,  the  component  of  the  normal  force  in  the  xz-plane  of 
symmetry  is  expressed  by 


<~Z>p  =  <CN>P  5S 

=  maV  sin  Op 


(40) 


Similarly,  a  normal  force  in  the  transverse  plane  is  in  evidence  during  a  sideslip,  or 


<*)p  =  (Cy)pqS 

=  -maV  sin 


(41) 


A  jet-induced  inflow  toward  the  jet  axis  at  the  tail  may  affect  the  stability  of 
the  aircraft  if  design  precautions  have  not  been  taken.  As  a  result  of  the  jet  exhaust 
spreading  out  behind  the  engine,  a  turbulent  mixing  of  the  air  outside  of  the  jet 
stream  with  the  Jet  exhaust  takes  place  along  the  boundary  of  the  jet  stream  (see 
Figure  11).  The  drawing  i*-  of  the  air  from  outside  the  stream  is  jet-induced  inflow. 

A  horizontal  tail  located  in  this  Jet-induced  inflow  field  will  be  subjected  to  jet- 
induced  downwash  angles.  Thus,  the  angle  of  attack  of  the  tail  would  be  modified  and 
pitching  moments  would  be  created  that  would  affect  the  stability  as  well  as  the 
equilibrium  of  the  aircraft. 

The  quantitative  effects  of  the  jet-induced  downwash  at  the  tail  can  be  calculated 
by  using  the  theory  developed  by  Ribner10.  This  theory  allows  for  curvature  of  the 
Jet  due  to  angle  of  attack  of  the  aircraft.  It  is  also  applicable  to  determination  of 
jet-induced  sidewash  of  the  vertical-tail  surfaces  at  asymmetric  power  conditions  or 
during  sideslip. 

In  the  absf>"ce  of  suitable  design  precautions,  such  as  boattailing  of  the  exhaust 
to  shield  the  il  surfaces  from  the  inflow  effect,  the  change  in  pitching  and  yawing 
moments  resulting  from  the  jet-induced  downwash  (Ao^  t  )p  and  sidewash  (Ao^  t  )  , 

respectively,  can  be  expressed  by 

<AMh.t.>p  =  (CNa)h.t.(A«h.t.)p^tSt*h.t.  <42> 

<ANv.t.>p  =  -<CN«>v.t.<AcV.t.)pqtStxv.t.  (43) 

where  xh  t  and  xy  t  are  negative  values  with  tails  aft  of  the  center  of  gravity. 


22 


The  variations  of  forces  and  moments  due  to  the  jet  engine  are  primarily  functions 
of  a  ,  yS  ,  and  V  ,  assuming  that  control  settings  are  constant.  From  the  preceding 
it  may  be  readily  deduced  that 


2(CNa)P 


m.V  cos  a„  3 


qS 


a  a 


qS 


(44a) 


z„  m.V  cos  ^  x„  3a„  «A.t.*h.t.  ^“h.t.*! 


^  zr  mav  cos  cl  x„  o a_ 

^Coa^P  “  Cr  -=■  +  - =£ - F  TT  +  'CNa^h- t. 


c  o  a 


qSc 


(44b) 


m„V  cos  B 


^(Cy^p  = - - (CNa)v.t. 


qS 


qSb 


(45a) 


2«Vp  =  - 


mftv  cos/3  xD  %  OA.t^v.t.  B<a7.t.)t 

— ^-f-(CNa)v.t.  - - ^ 


(45b) 


where  xp  is  positive  when  the  air  intake  is  forward  of  the  center  of  gravity  and 
xy  t  is  negative  wii;h  vertical  tail  aft  of  the  center  of  gravity,  and 


(Ccu)p 


(cmu)p 


2C 


cos  acos  /3; 


(46) 
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where 


<Cm>p 


n:aV  3ln  ap 
qS 


(CNa)v.  t.O^Cly.  t.  )p 


qSc 


(48) 


end 


3C t  Zp  +  masin  o^,  Xp 
3u  c  qS  c 


(49) 


3.4.6  Aeroelastic  Effects 

The  preceding  discussions  assumed  that  the  aircraft  was  rigid.  This  assumption  was 
permissible  in  the  past;  however,  modern  aircraft  flying  at  high  speed  under  dynamic- 
pressure  conditions  are  subject  to  degrees  of  flexibility  of  component  parts  which 
cannot,  at  times,  be  ignored  and  which  affect  the  stability  of  the  aircraft11"16.  The 
contribution  of  aeroelastic  deformation  to  derivatives  is  dependent  primarily  on 
aircraft  geometry  and  dynamic  pressure  as  well  as  structural  rigidity  and  Mach  number. 
Aeroelastic  phenomena  may  be  considered  in  two  separate  parts:  static  and  dynamic 
eeroelastic  effects. 


23 


When  aerodynamic  loading  takes  place  at  a  sufficiently  slow  rate  in  comparison  to 
the  natural  frequency  of  vibration  of  the  pertinent  part  of  the  structure  to  permit 
the  assumption  of  static  deformation  cf  the  structure,  the  influence  of  aeroelasticity 
can  be  accounted  for  by  modifying  the  derivatives.  Illustrations  of  steady-state 
distortions  which  have  been  serious  in  the  past  are  aileron  reversal  and  wing  divergence. 
Today,  such  factors  as  thinner  wings  and  more  flexible  fuselages  have  magnified  the 
effects  of  structural  flexibility  on  stability  and  control  of  aircraft. 

If  the  aerodynamic  loading  frequency  were  to  approach  the  structural  frequency  of 
the  pertinent  component,  me  structural  deformation  would  produce  perturbations  in  the 
aerodynamic  forces  and  moments  which  have  to  be  accounted  for  by  the  introduction  of 
additional  appropriate  derivatives  in  the  equations  of  motion  and  the  introduction  of 
additional  equations,  which  would  be  elasticity  equations. 

3.  I*.  7  Other  Effects 

The  preceding  discussion  has  included  major  factors  which  influence  analysis  and 
account  for  discrepancies  between  wind-tunnel  and  flight  data;  however,  it  does  not 
account  for  all  factors.  Other  factors  could  include  jet  pluming,  flow  separations 
associated  with  movements  of  shock  waves,  and  fuel  sloshing.  Since  one  never  knows 
what  phenomena  will  occur,  it  is  imperative  to  have  an  open  mind  in  trying  to  account 
for  discrepancies  in  comparisons  of  data. 

3. 5  Summary  of  the  Equations  of  Motion 

The  various  dynamic  relations  which  have  been  discussed  are  pertinent  to  an  under¬ 
standing  of  the  equations  of  motion  and  the  conditioning  of  data  to  the  equations.  The 
influence  of  power  and  structural  flexibility  on  the  various  aero'yr.&mic  parameters 
(coefficient  and  stability  derivatives)  was  stressed,  and  it  was  pointed  out  that  the 
net  result  cf  these  influences,  or  modifiers,  was  the  emergence  of  an  effective  aero¬ 
dynamic  parameter. 

It  is  easily  recognized  that  the  introduction  into  the  equations  of  motion  of  each 
individual  modifier  to  the  aerodynamic  parameters  would  result  in  a  cumbersome  set  of 
equations.  It  is  more  practical  to  let  the  normally  accepted  stability  symbol  (Cn^, 
for  example)  represent  the  effective  value  than  to  list  all  modifiers.  In  so  doing, 
one  should  be  aware  of  the  various  sources  which  contribute  to  the  magnitude  of  the 
effective  parameter  in  order  to  properly  account  for  these  contributions  during  an 
analog  investigation,  or  other  study,  in  which  wind-tunnel  and  calculated  data  are 
used.  On  the  other  hand,  in  the  inverse  problem  of  determining  coefficients  and 
derivatives  from  flight  data,  a  discrepancy  in  trends  as  well  as  magnitude  between 
wind-tunnel  and  flight  data  will  suggest  possible  influences  from  sources  not  accounted 
for  by  tunnel  data. 

3. 5.1  General  Equations 

The  following  assumptions  are  made  with  regard  to  the  equations  of  motion  summarized 
in  Table  IV: 

(i)  The  airplane  behaves  as  a  rigid  body,  in  that  the  moments  of  inertia,  inclina¬ 
tion  of  principal  axes,  etc.,  are  not  affected  significantly. 
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(ii)  The  airplane  is  symmetrical  about  the  xz-plane  with  regard  to  geometry  and 
mass  distribution. 

(iii)  The  axes  of  rotating  elements  on  the  aircraft  are  fixed  in  a  direction  relative 
to  the  body  reference  axes. 

(iv)  The  earth  is  flat.  Aircraft  speeds  are  assumed  to  be  insufficient  to  include 
earth  curvature  in  the  equations. 

(V)  The  forcing  frequency  of  a  disturbance  is  sufficiently  far  removed  from  the 
natural  frequency  of  the  pertinent  part  of  the  structural  components  to  permit 
the  disturbance  to  be  considered  as  a  static  load  and  the  effect  of  deformation 
to  be  accounted  for  by  modification  of  the  aerodynamic  parameters. 

(vi)  Each  aerodynamic  parameter  is  an  effective  parameter,  in  that  it  includes  all 
sources  contributing  to  its  net  value. 

Although  listed  in  Table  IV  for  completeness,  experience  has  shown  Cyp  ,  Cyr  , 

Cy^j  and  CCq  and  Cc^  to  be  normally  negligible. 

3.  5.2  Small -Perturbation  Equations 

The  general  equations  of  motion  in  Table  IV  are  suitable  for  analog  and  digital 
programing  which  involves  large  disturbances  and  nonlinear  terms;  they  are  not  suitable 
for  analytical  purposes.  For  such  purposes,  it  is  necessary  to  linearize  the  equations 
at  least  to  an  engineering  degree  of  accuracy.  This  is  accomplished  by  restricting 
their  applications  to  small  perturbations,  as  has  been  discussed  previously.  In 
addition,  the  perturbations  are  referred  to  a  secondary  spatial  reference  frame, 
discussed  in  Section  2.1.6,  which  is  the  unperturbed  airplane  axis  system  shown  in 
Figure  4(b).  Using  the  secondary  reference  system  for  small  perturbations  permits 
the  use  of  Equations  (31a,  b,  c),  which  simplifies  analysis  and  extends  the  validity 
of  the  linearized  perturbation  equations  to  maneuvers  involving  high  pitch  attitude 
and  large  bank  angles. 

The  uncoupled,  linearized  perturbation  equations  are  shown  in  Table  V  in  a  format 
which  generally  constitutes  the  basis  for  application  to  derivative  determination. 

The  assumptions  listed  for  the  general  equations  of  motion  are  also  valid  for  the 
equations  in  this  table.  In  addition,  it  is  assumed  that  the  maneuvers  are  such  as 
to  minimize  the  errors  in  the  g  terms  arising  from  the  approximation  of  the  gravity 
terms  shown  in  Equations  (40a,  b,  c).  Also,  it  is  assumed  that  the  gyroscopic  couples 
of  rotating  elements  are  not  significant,  which  may  not  always  be  the  case.  The 
equations  are  complete  within  the  limits  of  the  assumptions,  and  analysis  would  reveal 
all  modes  of  longitudinal  and  lateral  motions. 

Omission  of  the  longitudinal  firce  equation  and  the  terms  in  the  longitudinal 
equations  (50a,  b,  c)  would  remove  the  phugoid  mode  from  the  analysis  of  the  longitudinal 
motions,  leaving  only  the  short-period  mode.  This  short-period  format  of  the  long¬ 
itudinal  equations  is  the  one  usually  employed.  Although  the  small -perturbation 
equations,  shown  in  Table  V  are  frequently  used  in  the  format  shown  to  develop  relations 
for  derivative  determination,  it  is  also  desirable  to  list  the  equations  in  an  opera¬ 
tional  format  as  Laplacian  transforms  with  Laplace  operator  s  . 
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Using  Laplace  transforms  enables  the  dynamic  properties  of  the  airplane  to  be  defined 
by  a  series  of  transfer  functions  relating  the  various  responsive  motions  of  the  airplane 
to  disturbing  inputs.  The  transfer  functions  are  extensively  used  in  stability  and 
control,  handling  qualities,  and  automatic  flight  control  investigations  to  assess  the 
effects  of  configuration  changes,  the  effects  of  particular  stability  derivatives,  and 
the  effects  of  changes  in  automatic  control  systems.  They  are  also  helpful  in  obtaining 
stability  derivatives  from  flight  data. 

With  zero  initial  conditions  and  inputs  due  only  to  control  deflections,  the  Laplace 
transforms  of  the  small-perturbation  equations  of  motion  take  on  the  operational  forms 
shown  in  Table  VI  as  Equations  (62a,  b,  c)  and  (63a,  b,  c).  The  notations  Xp  ,  Ma  , 
etc.,  shown  in  the  equations,  are  a  convenient  means  of  listing  the  parameters. 

3.6  Determination  of  the  Roots  of  the  Determinant  of  the 
Lateral -Directional  Small -Perturbation  Equations 

The  following  discussion  regarding  the  determination  of  the  roots  of  the  determinant 
of  the  lateral -directional  small -perturbation  equations  is  based  on  Reference  17. 

Although  the  main  points  are  brought  out  at  this  time,  recourse  should  be  made  to  the 
reference  for  more  detailed  considerations. 

3.6.1  The  Determinant 

Using  the  Laplace  transform  format  of  the  lateral -directional  equations  (Equations 
(63a,  b,  c)  in  Table  VI),  the  determinant  of  these  equations  may  be  expressed  in  either 
of  two  formats,  as  follows: 

(i)  When  expressed  as 


As“  +  Bs3  +  Cs2  +  Ds  +  E  =  0  , 


then 

A  =  KK  - 1 

B  =  (Lp  +  IjNp)  +  (Nr  +  l'Lr>  -  (I  - 

C  =  (NpLr  -  NrLp)  -  (Lp  +  I^Np)^  -  (Nr  +  l$r)Yp  - 
-  t(l  -  Ij  sin  <x)Ng  -  (sin  a  -  i')!^] 

D  =  (N/p  -  NpLyj)  +  (NjXp  -  NpLr)Y0  +  +  l£g)  - 

~  82(^3  +  I^N/j)  +  .  -  Njii^g)  sin  a 

E  =  g^LgNp  -  LpNg)  -  gz(N/SCr  -  . 


(65) 


(66) 
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(il)  When  expressed  as 


s4  +  bs3  +  cs2  +  ds  +  e  =  0  , 


(67) 


then 


b  = 

c  =  -  (N'l;  -  n;l')  +  L'Yg  +  NyY^g  +  -  Vp  sin  a 

d  =  -  (NgLp  -  N'Lp  -  (R;l'-  n#)?^  -  +  • 

~  (N^L;  -  N£Lg)  sin  a 

e  =  -  61(L>gNp  -  L#&)  +  b2(N^  -  N^)  , 


where  the  primed  values  are  equal  to 


N,  +  ilL 

R/  =  — - rr  and 

1  1  -  I'l' 

x  z 


Ei 


Lj  + 

1  -  I'l'  ’ 
x  z 


(68) 


(69) 


The  determination  of  the  roots  of  the  determinant  is  dependent  upon  the  modes  of 
motion  of  the  aircraft.  The  modes  may  be: 

(a)  Lateral  phugoid  (coupling  of  spiral  and  roll  modes)  and  Dutch  roll. 

(b)  Spiral  divergence,  roll  subsidence,  and  oscillatory  (Dutch  roll). 


3.6.2  Determination  of  the  Roots  when  Lateral 
Phugoid  and  Dutch  Roll  Modes  Exist 


The  determinant  (Eqn.(67))  can  be  approximated  by  the 


following  biquadratic 


(70) 


in  which  the  first  and  second  quadratics  represent  the  Dutch  roll  and  lateral  phugoid 
modes,  respectively. 

Two  sets  of  conditions  must  be  satisfied  if  Equation  (70)  is  to  be  applicable17: 

(a)  The  approximate  nature  of  Equation  (70)  requires  that  e/c2  «  1  and  bd/c2  «  1 
to  assure  validity  of  the  equation. 

(b)  It  is  necessary  that  d2  -  4ed  <  1  in  order  that  the  lateral  phugoid  exist. 

Reference  17  points  out,  on  the  basis  of  limited  experience,  that,  for  values  of  e/c2 
of  approximately  0.05  or  less,  bd/c2  can  be  as  large  as  0.25  and  d2/4ec  as  low  as 
0.005  without  compromising  the  engineering  accuracy  of  Equation  (70).  Thus,  the 
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equation  applies  if 


e 


bd 

-7  <  0.25  . 
cz 


0.005 


4ec 


(71) 


The  second  quadratic  in  Equation  (70)  expresses  the  lateral  phugoid  very  simply; 
thus, 


from 


and 


^ph 


2^phc^>ph  —  “ 


(72) 


The  first  quadratic  in  Equation  (70)  is  unwieldy.  It  is  simpler  to  determine  the 
characteristics  of  the  Dutch  roll  mode  by  the  following  factored  form  of  determinant 


(s?  +  2£aJns  +«z)(s2  +  2Cph«nphs  +"nPh>  =  0  •  (73> 

Expansion  of  tiils  determinant  and  comparison  with  Equation  (67)  shows  that 


b  =  2&a  +  2Cphwnph 
c  =  wn  +  ^pt  +  <2^n>  <2^prf°nph> 
d  =  OCph^phK  +  <2H>"nph 
e  =  ‘tfptftf  • 


(74) 


Since  «nph  and  (2Cpllwnph)  are  obtained  from  Equation  (72),  and  (2 £a>n)  can 
now  be  determined  from  Equations  (74)  or 


2£a'n  =  b  -  2Cphc^nph 
wn  =  c  -fJnph  “  (2H>(2£phwnph)  • 

* 


(75) 


3.6.3  Determination  of  the  Roots  when  Spiral  Divergence, 
Roll  Subsidence,  and  Dutch  Roll  Modes  Exist 


When  the  spiral  divergence,  roll  subsidence,  and  Dutch  roll  modes  constitute  the 
lateral-directional  characteristics  of  the  ai  rplane,  which  i*  normally  the  case,  the 
determinant  as  represented  by  Equation  (67)  may  be  factored  in  the  following  terms 
characterizing  these  modes: 


(76) 
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The  coefficients  b  ,  c  ,  d  uid  e  (Equation  (68))  in  terms  of  the  factors  of 
Equation  (76)  are 


.  r  1  1 

b  =  2  +  —  + 


n  T  T 
*r  *s 


C  =  CO*  +  2&  (■—  +  —  ]  +  — — 

n  n  \tr  v  Vs 


,  A  A 


V. 


e  -  a)_ 


T„T 


R  S 


(77) 


When  the  spiral -mode  root,  1/TS  ,  is  much  less  than  the  roll  subsidence  root, 
1/Tr  ,  as  it  usually  is,  the  coefficients  b  ,  c  ,  d  ,  and  e  may  be  approximated 
to  a  good  degree  of  accuracy  by 


b  ~  2  lco  +  — 

tr 


c  -  <  +  =" 
*R 


d  ~  'ol  —  . 

n  -1 

XR 


(78) 


Eliminating  lAR  in  Equation  (78) 


c  ~  (24c<U  -7  +  a>l 

CO*  u 

n 


d  ~  (2 ico  )  +  —  . 

n  or 


(79) 


Eliminating  2£&>n  in  Equations  (79)  provides  an  accurate  solution  of  w*  within  the 
limitation  that 

1  «± 

Ts  tr 

or 


=  0  . 


(w„)3  -  c(cu*)z  +  bd {col',  ~  d? 


(80) 
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Eliminating  cojj  in  Equations  (79)  to  solve  for  2£<x>2  within  the  limitation  that 
1/TS  «  1/Tr  results  in 

(2 £o>n)3  -  2b(2C^n)2  +  (c  +  b2)(2^n)  +  (d  -  cb)  =  0  .  (81) 

The  roll-subsidence  root,  1/TR  ,  may  now  be  obtained  from  coefficient  b  or  d 
In  Equation  (78)  or 


1 


=  b  -  2£o.'n  . 


(82) 


The  spiral -divergence  root.  1/TS 
coefficient  expressions  in  Equations 


,  may  now  be  approximated  from  any  one  of  the 
(77),  such  as 


(83) 


4.  MASS  CHARACTERISTICS 

The  airplane  mass  characteristics  -  weight,  location  of  the  center  of  gravity, 
moments  of  inertia,  and  inclination  of  principal  axis  -  significantly  affect  airplane 
motions.  Errors  in  the  knowledge  of  these  quantities  are  reflected  directly  in  the 
flight-determined  derivatives  and  may  govern  the  validity  of  the  derivatives  in  com¬ 
parisons  with  wind-tunnel  data.  Although  possible  inaccuracies  in  the  knowledge  of 
the  inertia  characteristics  must  be  given  serious  consideration  in  comparisons  of 
flight -determined  derivatives  with  wind-tunnel  data,  these  derivatives  have  been  used 
effectively  in  flight-guidance  simulator  studies. 

The  weight  and  horizontal  location  of  the  center  of  gravity  are  always  determined 
experimentally.  Inasmuch  as  the  vertical  location  of  the  center  of  gravity,  moments 
of  inertia,  and  location  of  the  principal  axis  are  difficult  to  determine  experimentally, 
manufacturer1  s  estimates  are  usually  relied  upon.  These  estimates  are  considered  to 
be  cf  sufficient  accuracy  for  most  work  involving  flight  tests.  If  more  precise  data 
are  required,  they  should  be  determined  by  using  experimental  techniques. 

It  would  be  highly  desirable  to  determine  all  of  the  mass  characteristics 
experimentally.  This  is  not  always  feasible  because  of  the  lack  of  proper  facilities. 
Large,  flexible  aircraft,  such  as  the  Boeing  B-52,  offer  practical  problems,  in  that 
experimentally  determined  rolling  moments  of  inertia  with  wings  drooped  would  not  be 
representative  of  flight  conditions.  The  following  discussion  of  the  experimental 
determination  of  mass  characteristics  of  aircraft  is  intended  to  serve  as  a  guideline 
in  setting  up  suitable  facilities  for  use  with  most  categories  of  aircraft. 


4.1  Weight  and  Center -of-Gravity  Location 


The  weight  and  longitudinal  position  of  the  center  of  gravity  relative  to  the 
horizontal  reference  line  of  the  airplane  for  the  empty  and  gross  weight  conditions 
can  be  obtained  easily  by  leveling  the  airplane  on  suitable  scales  or  electronic 
weighing  cells.  With  weighing  cells,  two  of  the  cells  (Rj  and  R2)  are  usually  located 
at  the  wing  jackpoints  and  the  third  cell  (R3)  is  located  at  some  convenient  distance, 
I  ,  forward  or  aft  of  the  wing  Jackpoints.  The  horizontal  position  of  the  center  of 
gravity  relative  to  the  Jackpoints  is  then  determined  from 


AZ 


R3Z 
2r  ' 
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For  aircraft  operating  on  conventional  fuels,  the  variation  of  the  center  of  gravity 
with  fuel  consumption  can  usually  be  defined  adequately  by  weighing  the  airplane  at 
several  fuel  levels,  providing  there  is  a  predetermined  sequence  or  mode  of  operation 
in  obtaining  the  fuel  from  the  various  fuel  cells.  When  the  aircraft  is  equipped  with 
fuel  cells  from  which  the  fuel  can  be  drawn  selectively,  the  center  of  gravity  position 
becomes  a  function  of  the  sequence  in  drawing  off  the  fuel  from  the  various  cells  as 
well  as  the  weight  of  the  fuel.  In  some  instances,  it  has  been  found  necessary  to 
account  for  fuel-tank  shape  and  airplane  attitude.  Where  hazardous  fuels  are  used, 
the  center  of  gravity  is  determined  experimentally  for  the  no-fuel  condition  only;  the 
effect  of  fuel  on  the  center  of  gravity  position  is  calculated.  The  horizontal  location 
of  the  center  of  gravity  is  experimentally  obtained  at  least  to  within  0.01  mean  aero¬ 
dynamic  chord,  which  is  considered  adequate  for  derivative  determination. 

During  flight  tests,  the  center  of  gravity  is  obtained  by  observing  the  total  amount 
of  fuel  consumed  and  subtracting  it  from  the  takeoff  weight.  Reference  to  a  chart 
showing  the  variation  of  weight  with  center  of  gravity  provides  the  desired  answer. 


An  accurate  knowledge  of  the  vertical  location  of  the  center  of  gravity  is  pertinent 
to  the  experimental  derivative  studies,  insofar  as  experimental  determination  of  moments 
of  inertia  and  comparison  with  wind-tunnel  data  are  concerned.  The  vertical  center  of 
gravity  can  be  obtained  by  static  or  oscillatory  techniques.  For  the  static  test 
techniques,  the  airplane  is  placed  in  various  pitch  or  roll  attitudes.  For  the  roll 
approach  (Fig.  12),  the  airplane  is  mounted  in  a  horizontal,  wings-level  attitude  on 
knife  edges  alined  with  respect  to  each  other  in  the  plane  of  symmetry  of  the  aircraft. 
By  rolling  the  airplane  to  various  attitude  angles  and  measuring  the  reaction  R1  , 
moment  arm  ,  and  the  roll  angle  <t>  ,  using  a  clinometer,  the  vertical  position  of 
the  center  of  gravity  is  obtained  from  the  equation 


RiYi  ~  WcZcsln(^ 
W  sin<£ 


(85) 


For  rigid  aircraft  of  the  order  of  15,000  lb,  and  under  carefully  controlled  conditions, 
the  vertical  position  is  considered  to  be  determinable  to  within  1  inch. 


To  determine  the  vertical  position  of  the  center  of  gravity  from  free-oscillati-  a 
tests,  any  one  of  several  techniques  may  be  used.  The  simplest  technique  consists  of 
changing  the  equivalent  torsional  spring  constant  for  pitching  or  rolling  moment  of 
inertia  tests.  For  rolling-oscillation  tests  with  the  setup  shown  in  Figure  13  and 
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with  small  damping  effects  -  a  necessary  condition  for  successful  tests  -  the  equations 
of  motion  for  the  two  spring  conditions  are 

dx  +  Ixc  t  mz2  +  mcz 2)<j6j  +  (K^  -  Wz  -  W^)^  =  0  (86a) 

(Ix  'r  IXc  +  mz2  +  mczc)<jt>2  +  (Kt2  -  Wx  -  Wcz c)<f>2  =  0  .  (86b) 

Considering  4>1  =  Acosojjt  and  -  BcosoJ2t  ,  it  is  found  upon  solving  Equations 
(86a,  b)  for  z  ,  the  vertical  distance  from  the  knife  edge  to  the  center  of  gravity, 
that 


w[l  -  (Pj/p2)2]  w 

The  equivalent  torsional  spring  constant,  Kt  ,  may  be  changed  from  Ktt  to  Kt2  by 
changing  the  linear  springs  or  the  distance  a  which  Is  perpendicular  to  the  spring 
(see  Figure  13).  The  change  in  linear  springs  is  probably  the  more  desirable  approach. 

Inasmuch  as  the  rolling-oscillation  test  setup  discussed  constitutes  an  inverted 
pendulum,  it  is  imperative  that  the  equivalent  torsional  spring  constant,  Kt  ,  be 
greater  than  Wz  +  Wcz„  for  stability  of  setup.  Also,  the  accuracy  of  the  results 
depends  upon  avoiding  secondary  spring  actions  of  tiebacks  and  structural  flexibility, 
which  could  inadvertently  result  in  a  lower  effective  spring  constant  than  expected 
because  of  an  equivalent  sories  action  of  the  secondary  unwanted  spring  action  with 
the  intended  spring. 

4.2  Moments  of  Inertia 

The  moments  of  inertia  of  an  airplane  are  usually  calculated  during  the  design 
phase  and  are  based  on  estimated  weights  and  centroid  locations  for  various  parts  of 
the  aircraft.  These  calculated  moments  of  inertia  are  considered  to  be  adequate  for 
most  analyses  when  the  results  are  to  be  used  in  simulator  studies.  However,  should 
experimental  determination  of  the  inertia  be  required,  methods  are  available  (see 
References  18  to  21).  The  methods  are  generally  restricted  to  rigid  aircraft  and  to 
aircraft  whose  weight,  as  well  as  the  safety  precautions  of  the  experiment,  will  permit 
pivoting  the  aircraft  on  knife  edges  and  suspending  it  from  overhead  cables. 

Schematic  representation  of  typical  methods  for  determining  the  rolling  and  pitching 
moments  of  inertia  are  illustrated  in  Figures  i3  and  14,  respectively.  Equation  (86) 
is  applicable  to  the  determination  of  rolling  moments  of  inertia  in  accord  with  Figure 
13,  with  consideration  given  to  the  proper  interpretation  of  the  lengths  z  and  zc 
to  the  mountings  shown.  In  Figure  14,  cradle  weight  is  zero.  The  yawing  moment  of 
inertia  may  be  safely  determined  from  a  cable-suspension  method  used  to  determine  the 
inclination  of  the  principal  axis  (Figures  15  and  16),  which  is  discussed  subsequently. 

Unless  precautions  are  taken  in  every  detail  of  an  experimental  setup,  difficulties 
may  be  encountered  because  of  flexibility  of  experimental  components,  which  will  alter 
the  effective  spring  constant,  Kt  ,  or  modify  the  free-oscillation  pivotal  point 
relative  to  the  center  of  gravity  of  the  aircraft.  In  one  instance  of  determining  the 
pitching  moment  of  inertia  when  the  aircraft  was  supported  at  the  wing  jackpoints  and 
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oscillated  with  the  spring  at  the  nose,  the  wing  section  which  had  been  considered  rigid 
was  observed  to  flex  as  the  aircraft  oscillated.  This  flexing  caused  the  axis  of 
rotation  to  shift  forward  and  downward  from  the  line  through  the  jackpoints. 

A  common  fault  is  the  use  of  flexible  cables  as  tiebacks  for  the  springs  and  con¬ 
nection  from  the  spring  to  the  aircraft.  Under  no  condition  should  flexible  connections 
be  used,  inasmuch  as  they  constitute  springs  in  series  with  the  actual  intended  springs 
employed;  thus,  the  system  from  tieback  to  aircraft  represents  a  much  softer  spring 
than  intended.  It  should  also  be  noted  that,  on  some  aircraft,  attaching  the  spring  to 
the  aft  portion  of  the  fuselage  would  be  an  error,  since  the  aft  portion  of  the  fuselage 
would  constitute  a  relatively  flexible  structure  and  alter  the  effective  spring  constant. 

Serious  errors  can  also  result  when  knowledge  of  the  center-of-gravity  location  is 
inaccurate  and  when  the  line  of  action  of  the  spring  from  the  attach  point  to  the  air¬ 
craft  is  not  perpendicular  to  the  plane  formed  by  the  axis  of  rotation  and  the  point  of 
spring  attachment  on  the  aircraft  (Pig. 13). 

Generally,  the  inertia  characteristics  are  determined  for  no-fuel  conditions  because 
fuel  sloshing  tends  to  bring  in  a  beat  action  in  the  oscillatory  motions.  When  deter¬ 
mination  is  attempted  with  fuel  onboard,  the  difference  in  oscillatory  modes  between 
the  sloshing  fuel  and  the  aircraft  should  be  as  large  as  is  practical,  with  due  regard 
to  safety  of  the  setup,  to  minimize  the  beat  action  and  permit  determination  of  the 
natural  frequency  of  oscillation  of  the  aircraft. 

Measuring  the  inertias  of  very  large  aircraft  is  difficult  and  is  compounded  with 
flexible  aircraft.  Such  measurements  are  not  in  the  realm  of  the  methods  discussed. 

A  unique  facility  designed  to  enhance  the  feasibility  for  determining  the  moments  of 
inertia  of  large  aircraft  about  all  three  axes  is  located  at  the  US  Air  Force  Flight 
Test  Center,  Edwards,  California,  USA.  Its  capabilities  cover  a  weight  range  from 
30,000  to  300,000  b  and  moments  of  inertia  from  250,000  to  10  x  106  slug  ft2. 

The  facility  enables  the  determination  o'  aircraft  moments  of  inertia  from  measure¬ 
ments  of  changes  in  pendulum  characteristics  resulting  from  the  addition  of  an  aircraft 
to  a  freely  oscillating  platform.  The  basic  elements  of  the  facility  consist  of  the 
platform,  a  control  console  for  activating  varicus  systems  which  ready  the  platform 
for  oscillation,  and  an  instrumentation  console  for  regulating  the  amplitude  and 
measuring  the  period  of  the  oscillations.  The  platform  is  an  integral  cruciform 
structure  110  ft  long  and  80  ft  wide,  with  its  loading  surface  flush  with  the  surround¬ 
ing  floor  space.  The  apparatus  employs  special  hydrostatic  bearings  (identical  to 
those  used  in  the  200  in.  Paloraar  telescope)  to  support  the  platform,  which  is  lockable 
in  two  axes  with  oscillation  about  the  axis  of  interest. 

To  contend  with  the  problem  of  aircraft  flexibility,  stiffening  jacks  are  used  to 
support  the  aircraft  structure.  As  a  result  of  the  stiffening  operation,  flexibility 
effects  are  considered  to  be  less  than  4%  in  roll  and  2 %  in  pitch. 

The  experimental  error  in  the  methods  discussed  is  of  the  order  of  ±5%  or  less. 


4.3  Inclination  of  Principal  Axis 

The  inclination  of  the  principal  axis  of  the  airplane  is  one  of  the  more  difficult 
quantities  to  determine  experimentally.  An  error  of  1/4°  in  the  value  of  the  inclina¬ 
tion  of  the  principal  axis  can  significantly  affect  some  of  the  derivatives.  The 
method  of  Reference  22  is  considered  accurate  to  1/6°. 
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Ibis  method  consists  of  finding  the  direction  of  the  restoring-moment  vector  which 
produces  no  rolling  moment  relative  to  the  body  x-axis  during  the  yawing  oscillations 
of  the  airplane  as  a  spring  mass  system  while  suspended  by  means  of  a  cable  attached 
to  a  hoisting  sling.  Figure  15  shows  schematically,  and  Figure  16  shows  photogra¬ 
phically,  a  general  arrangement  of  the  setup.  The  airplane  is  suspended  at  a  horizontal 
pitch  attitude,  and  yaw  restraint  is  provided  by  two  sets  of  springs  whoBe  lines  of 
action  lie  in  a  common  plane.  The  springs  should  provide  a  pure  couple  action.  Hie 
restoring-moment  vector  acts  normal  to  the  plane  of  the  springs.  The  springs  may  be 
attached  to  short,  rigid  mounting  brackets  located  below  the  wings  equidistant  from  the 
plane  of  symmetry  or  to  brackets  mounted  below  the  fuselage  ahead  of  and  behind  the 
center  of  gravity.  In  this  respect,  the  wing  mounting  arrangement  is  most  convenient 
and  less  time-consuming.  It  is  essential  that  the  springs  provide  a  pure  couple  action. 


As  the  airplane  oscillates  in  yaw  with  various  inclinations  of  the  plane  of  the 
spring  couple  (angle  8  in  Figure  15),  some  coupling  is  present  between  yaw  N  and 

op 

roll  L  ,  which  results  in  a  certain  amount  of  rolling  oscillation.  This  is  shown  in 
the  following  equations  where  the  subscript  r  denotes  the  reference  attitude  of  the 
airplane: 


IxrPr  —  lxrzjXr  -  L 

(88) 

Izr*'r  “  IxrZrPr  =  N  • 

(89) 

At  some  one  value  of  8gJ}  ,  however,  the  rolling  motion  accompanying  the  yawing  motion 
is  zero  (Ipl/lrl  =  0)  .  In  this  situation  the  preceding  equations  reduce  to 


*xrzr^r  =  k 
Izrrr  =  N  . 


However,  as  shown  in  Figure  15, 


(90) 


Hence 


tan  8 


sp 


tan  8 
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(91) 


(92) 


Inasmuch  as  the  inclination  of  the  principal  axis  is  given  by  the  wellknown  expression 


tan  2e 
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xrzr 


Izr  " 


substitution  of  Equation  (92)  for  IXZr  in  Equation  (93)  gives 


2IZf  tan  S8p 


(93) 


(94) 


tan  2e 
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The  value  of  IZr  is  determined  as  a 


However,  IXr  must  be  determined  from  other  tests. 

Figure  17  shows  a  typical  experimental  plot  of  the  variation  of  |p|/|rj  with 
for  determining  the  value  of  5SI)  at  which  |p|/|r|  is  zero.  In  obtaining  the  tests 
points  shown  in  the  figure,  the  flight  test  roll-  and  yaw-rate  gyros  mounted  in  the 
airplane  were  used  to  obtain  oscillograph  records  for  determining  |p|/|r|  from  the 
tr&nsient  oscillations. 

The  measured  values  of  moments  of  inertia  relative  to  the  reference  axes  and  the 
determined  inclination  of  the  principal  axes  mzy  be  used  to  determine  the  principal 
moments  of  inertias,  IXQ  and  IZo  ,  by  using  the  following  equations 


I*o 

=  IXr  -  IXrzr  tane 

(93) 

xzo 

=  Izr  +  Ixrzr  tan  e  . 

(97) 

Although  no  mention  was  made  of  the  effects  of  air  mass  on  the  experimental  values 
of  moments  of  inertia,  the  effects  should  be  considered  and  corrections  applied  if 
necessary.  Reference  23  provides  formulas  to  correct  for  air-mass  effects. 

Formulas  for  transferring  moments  of  inertia  from  one  set  of  axes  to  another  were 
presented  in  Section  3. 1. 


byproduct  of  the  test  by  using 
C  cos  S 

'V.  °P 


cor 


(95) 


5.  INSTRUMENTATION 

Basic  to  an  analysis  of  flight  data  is  the  instrumentation.  Considerable  instrumen¬ 
tation  research  has  been  in  progress  and  many  flight  test  Instruments  have  been  deve¬ 
loped  to  improve  the  linearity  of  response,  resolution,  dynamic  response  characteristics, 
readability,  ruggedness,  and  reliability  of  calibrations  over  varying  operating  con¬ 
ditions  and  extended  periods  of  time.  In  addition,  the  application  of  the  instruments 
requires  knowledge  of  mounting  accuracy,  sources  of  error  in  the  flight  records,  and 
methods  of  correcting  the  errors.  Inadequate  appreciation  of  the  instrument  character¬ 
istics,  mounting  accuracy,  and  possible  influence  of  sources  of  error  serves  as  a 
detriment  to  the  successful  application  of  new  techniques  of  analysis  as  well  as  a 
detriment  to  the  analysis  by  approximate  methods. 

In  the  following  discussion,  sufficient  guidelines  are  presented  to  show  the  care 
required  in  the  selection,  installation,  and  calibration  of  instruments  to  minimize 
errors  in  the  analysis  of  flight  data.  Individual  instruments  may  differ  from  one 
organization  to  another  and  the  degree  of  sophistication  in  instruments  and  recorders 
will  vary  with  the  individual  investigation;  however,  the  principals  of  operation  of 
the  sensors  are  generally  the  same. 
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5.1  Kach  Nimber,  Altitude,  and  Dynamic  Pressure 

Accurate  determination  of  Mach  number  is  of  fundamental  importance  in  fligit  testing 
high-speed  aircraft.  The  principal  methods,  discussed  in  detail  in  References  24  and 
25,  are  based  upon  the  following  relationship  for  subsonic  conditions  (M  <  1.0) 

IT  *  >"> 


1  - 


-  i  =  (i  +  o.ai)7/2  -  l 


(98) 


For  supersonic  conditions  (M  >  1.0),  the  equation  is  modified  to  include  the  loss  in 
total  pressure  behind  the  shock  wave 


_£ 

P 


1/C/-D 


-  1  =  1.2M2 


5.76M2  \5/z 


5.6M2  -  0.8J 


-  1 


(99) 


The  impact  pressure  qc  and  the  static  pressure  p  are  measured  by  using  a  pitot- 
static  head  and  pressure  recorders.  The  maximum  Mach  number  as  well  as  dynamic 
pressure  which  can  be  determined  by  using  pitot -static  heads  is  of  the  order  of  3.5. 
digher  speeds  are  primarily  dependent  upon  inertial  platforms  and  radar.  Dynamic 
pressures  at  Mach  numbers  is  the  approximate  range  of  2.5  to  8.0  can  be  determined 
through  the  use  of  a  spherical  flow-direction  sensor  and  a  total-  (stagnation)  pressure 
technique. 


5.1.1  Pitot-Static  Head  (M  <  3.5) 

Much  research  has  been  done  on  various  types  and  configurations  of  total -pressure 
heads  to  reduce  angularity  effects26' 27.  The  type  shown  in  Figure  18  is  used  widely. 
This  head  has  an  external  cylindrical  shape,  a  cylindrical  chamber,  and  a  10°  slant 
profile.  It  is  insensitive  (zero  error)  to  angle-of-attack  from  -5°  to  20°  and  up  to 
10°  of  sideslip.  The  error  is  less  than  \%  in  the  angle-of-attack  range  from  -10°  to 
25°  and  ±10°  sideslip. 

The  arrangement  of  the  static-pressure  orifices  on  the  head  has  been  found  to  be 
pertinent,  in  increasing  the  range  of  insensitivity  of  the  orifices  to  flow  angularities. 
The  arrangement  used  has  been  determined  from  tests  of  orifice  configurations20'  29. 

The  two  identical  sets  or  arrangements  shown  in  Figure  18  are  each  circumferential, 
with  four  orifices  on  the  top,  six  on  the  bottom,  and  one  on  the  bottom  centerline 
behind  the  others.  The  two  sets  of  static-pressure  orifices  are  used  to  provide  for 
separate  pressure  systems.  One  set  of  static  orifices  is  used  for  the  pilot’s 
instruments,  the  other  set  for  flight  test  recording  instruments  to  minimize  the  time 
lag  of  response  that  would  be  encountered  with  a  common  system.  The  arrangement  of 
the  orifices  in  each  set  provides  an  increased  range  of  insensitivity  to  angle-of- 
attack;  however,  it  is  not  as  insensitive  to  sideslip.  Large  static-pressure  errors 
are  encountered  at  sideslip  angles  greater  than  3°,  Since  constant  sideslip  angles 
are  seldom  encountered,  the  static -pressure  data  c*n  be  readily  faired. 
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Installation  of  the  pitot-static  head:  Installation  of  the  pitot-static  head  requires 
consideration  of  the  coaplicated  flow  field  of  the  airplane,  which  is  a  function  of  the 
airplane  configuration  as  well  as  Mach  number  ana  attitude.  Errors  in  pitot-static- 
head  readings  resulting  from  this  flow  field  are  referred  to  as  petition  errors.  The 
static-pressure  orifices  are  particularly  affected  hy  position  erro:'*,  at  subsonic 
speeds;  thus,  precautions  are  taken  to  mount  the  pitot-static  head  as  far  ahead  of  the 
airplane  as  is  practical. 

Of  the  various  types  of  installations  of  the  pitot-static  head  -  such  as  nose  boom, 
wing  boom,  and  fuselage  -  the  nose-boom  installation  is  the  most  suitable  for  minimizing 
position  errors.  In  this  installation,  shown  in  Figure  19,  the  head  is  mounted  on  a 
boom  extending  as  far  ahead  of  the  nose  of  the  airplane  as  is  practical.  As  reported 
in  Reference  29,  the  amount  of  error  in  Mach  number  due  to  position  error  in  the  static- 
pressure  measurements  can  be  related  to  certain  physical  measurements  on  the  airplane. 
This  is  shown  in  Figures  20  and  21,  which  are  reproduced  from  Reference  30.  In  Figure 
20,  the  error  in  Mach  number  due  to  static-pressure  error  is  plotted  as  a  ratio  of 
boom  length  to  the  maximum  effective  fuselage  diameter  for  subsonic,  transonic,  and 
supersonic  speeds.  In  Figure  21,  the  variation  in  Mach  number  error  with  Mach  number 
is  plotted  for  two  airplanes  having  boom-length-to-fuselage-diaoeter  ratios  of  0.60 
and  0.95.  Above  a  Mach  number  of  1.05,  the  position  error  drops  to  zero.  The  Mach 
number  at  which  the  position  error  drops  to  zero  is  dependent  upon  the  nose-boom 
geometry  and  is  the  Mach  number  at  which  the  shock  wave  ahead  of  the  airplane  crosses 
over  the  static-pressure  orifices. 

Wing -boom  installations  of  the  pitot-static  head  are  subject  to  several  disadvantages, 
including  possible  susceptibility  to  the  shock  wave  caused  by  the  wing  as  well  as  the 
shock  wave  caused  by  the  fuselage.  This  complicates  the  calibration  and  makes  it  more 
difficult  for  the  pilot  to  fly  at  the  desired  Mach  number  in  the  regions  where  the 
shoc't  waves  are  in  the  vicinity  of  the  orifices.  Wing  booms  are  usually  more  sensitive 
to  sideslip  and  subject  to  more  lag  in  response  because  of  the  longer  tubing  required. 

Fuselage  installations  of  the  head  are  subject  to  position  errors,  which  are  diffi¬ 
cult  to  estimate. 

Calibration:  Calibration  of  the  pitot-static  head,  fortunately,  involves  only  the 
determination  of  the  position  error  for  the  static  pressure  -  the  total  pressure  is 
not  affected  by  position  error.  Various  methods  that  have  been  used  include  the  pacer 
method,  the  fly-by  (tower-pass)  method,  and  modifications  of  the  basic  radar-photo - 
theodolite  method24.  The  pacer  method  requires  the  use  of  a  pacer  airplane  with  a 
calibrated  system  and  special  flights  for  calibration  purposes.  The  fly-by  method 
requires  lg  flight  at  extremely  low  altitudes  past  an  instrumented  course.  This 
latter  method  not  only  requires  special  flights,  but  is  hazardous  and  limited  to 
Mach  numbers  of  about  0.8. 

The  radar -phototheodolite  method  has  the  advantage  of  providing  calibration  data 
during  routine  research  flights.  The  method  makes  use  of  a  radiosonde  unit  to  measure 
static-pressure  and  temperature  variations  of  the  atmosphere  with  altitude.  It  also 
requires  ground  equipment  consisting  of  a  radar  unit,  a  phototheodolite,  a  chronograph, 
and  three  cameras.  One  of  the  cameras  photographs  the  radar  scope  and  gives  the  slant 
range;  the  target  camera  gives  the  correction  to  the  elevation  scales;  and  the  third 
camera  gives  the  elevation  scale.  The  airplane  itself  is  equipped  with  a  radar  beacon 
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to  assist  in  tracking.  The  three  cameras  and  the  airplane’s  internal  records  are 
synchronized  by  means  of  the  chronograph.  Hie  radar-phototheodolite  unit  determines 
the  range  and  elevation  angle  of  the  airplane  from  which  the  true  geometric  altitude 
of  the  airplane  is  determined  .within  ±100  ft)  as  a  function  of  time. 

A  cross  plot  of  the  radar -phototheodolite  data  (airplane  altitude  vs.  time)  with  the 
radiosonde  results  (free-stream  static  pressure  vs.  altitude)  provides  a  plot  of  true 
free-stream  static  pressure  as  a  function  of  time.  Since  the  time  base  of  the  air¬ 
plane*  s  indicated  static -pressure  records  is  synchronized  with  the  radar-phototheodolite, 
a  comparison  of  the  airplane’ s  indicated  static-pressure  records  with  the  cross  plot 
provides  the  position  error,  Ap  ,  of  the  static  head.  The  corrected  static  pressure 
may  now  be  obtained  from  the  relation  p  =  p^  +  Ap  .  The  true  Impact  pressure  is  now 
determined  from  qc  =  pT  -  p  =  qci  -  Ap  . 

True  Mach  lumber:  True  Mach  number  is  determined  from  tables  of  qc/p  as  functions 
of  Mach  number  based  on  Equations  (98)  and  (99).  The  indicated  Mach  number,  M1  ,  as 
determined  from  i>£  ,  qci  ,  and  the  tables,  is  plotted  against  the  corrected  Mach  number, 
M  ,  tc  provide  a  calibration  curve,  such  as  shown  in  Figure  22,  for  the  pitot-static- 
head  installation  on  the  airplane.  Generally,  calibration  data  points  for  four  or 
five  flights  are  used  before  the  calibration  curve  is  finalized.  The  scatter.  Am  , 
in  calibration  points  is  usually  within  ±0.01  at  subsonic  and  supersonic  speeds  and 
within  ±0.02  at  transonic  speeds. 

Pressure  altitude:  Altitude  is  generally  expressed  in  terms  of  "‘pressure  altitude”, 
which  is  the  altitude  in  the  standard  atmosphere  tables  corresponding  to  the  corrected 
static  pressure.  The  corrections  for  a  given  pitot-static  pressure  system  are  obtained 
in  the  form  Pj/p  vs.  M  .  The  curve  for  this  relationship  is  derived  from  the  Mach 
number  calibration  of  the  system  and  the  position  error  for  the  static  pressure  deter¬ 
mined  as  a  ratio  of  the  true  static  pressure  by  the  following  equations  from  Reference 
25: 

When  M  <  1.0  , 

Ap 

P 

When  M  >  1.0  , 

Ap 

T  *  V5.8M2  -  0.8  "  'AT/' 

In  routine  tests,  the  pressure  ratio  Pj/p  is  divided  by  pt  to  obtain  p  ,  which 
is  used  to  determine  the  pressure  altitude. 

Dynamic  pressure:  The  dynamic  pressure,  q  ,  is  determined  from  the  simple  relation 

q  =  0,7pM2  . 
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5.1.2  Use  of  Spherical  Flow-Direction  Sensor  to 
Obtain  Dynamic  Pressure 

In  the  absence  of  true  Mach  number,  such  as  when  flight  Is  beyond  the  practical 
limit  of  the  pitot-static  tube  (M  ~  3.5),  a  technique  has  been  evolved  to  obtain  the 
dynamic  pressure,  in  the  higher  supersonic  and  hypersonic  regions,  directly  from  the 
total -pressure  port  of  a  spherica1  flow-direction  sensor31.  The  flow-direction  sensor, 
described  in  more  detail  in  Section  S.3,  is  a  movable  sphere  mounted  at  the  nose  of 
the  airplane  to  form  a  "ball  nose".  The  total-pressure  port  vectors  into  the  resultant 
velocity  at  the  sensor. 

Inasmuch  as  q  =  0.7  pM2  and,  from  the  Rayleigh  pitot  formula, 


—  =  f (M)  , 
Pr 


2 

y/(y+0 

zyn2  -  (7-  1) 

(y  +  dm2 

7  +  1 

i/(y+i) 


the  dynamic  pressure  can  be  expressed  as 

3  =  0.7p[f(M)]p 


—  =  F(M)  . 
PT 


A  plot  of  q/pT  versus  M  (Pig. 23(a))  shows  that  this  ratio  varies  only  about  5%  in  the 
Mach  range  above  2.5.  As  l  result  of  this  small  variation  in  q/PT  at  the  higher  Mach 
numbers,  it  was  suggested  that  an  "indicated”  dynamic  pressure,  qA  ,  could  be  expressed 
as 


Figure  23(b)  shows  the  ratio  of  indicated  to  true  dynamic  pressure,  q^/q  ,  for  two 
values  of  K  .  Using  K  =  0.  526  ,  q  is  5%  high  at  M  =  2. 1  and  2.5%  low  at  M  =  7  . 

5.1.3  Pressure-Recording  Instruments 

Selection  of  the  pressure -recording  instruments  an  !  their  ranges  for  a  given  instal¬ 
lation  depends  on  the  altitude  and  Mach  number  range  i  .er  which  a  specified  attainable 
Mach  number  accuracy  is  desired.  When  tests  are  to  be  conducted  at  one  altitude,  it 
is  no  problem  to  select  a  pressure-recording  instrument  to  provide  the  requisite 
accuracy.  For  tests  conducted  over  a  large  range  of  altitudes,  the  requisite  accuracy 
may  be  attained  by  using  a  combination  of  limited-range  instruments.  Considerations  of 
the  pressure  time  lag  require  that  the  Instrument  volume  remain  as  small  as  possible, 
thus  necessitating  an  evaluation  of  instrument  accuracy  with  consideration  for  the 
errors  caused  by  the  added  time  lag  of  multiple-instrument  installation. 


39 


\ 


The  lag  In  response  at  the  recorder  servo  or  pilot  display,  as  the  case  nay  be,  can 
be  calculated  from  the  following  formula  (from  Reference  32),  which  takes  into  account 
the  sense  line,  Instrument  volume,  and  pressure 


128m0  l(VB 1) 
77pD'S 


(107) 


where 

X  =  lag  in  response,  sec 
M0  =  viscosity  of  the  fluid,  lb  sec/in*2 
I  =  length  of  sense  line,  in 
Vol  =  instrument  volume,  in’3 
P  =  mean  pressure  in  sense  line,  lb/in2 
D  =  diameter  of  sense  lir.e,  in  . 

Several  ranges  of  instruments  are  available  for  both  the  static-pressure  and  total- 
pressure  recorders.  For  the  static-pressure  recorders,  the  lower-range  instruments 
require  temperature  calibration.  Hysteresis  and  friction  errors,  and  temperature 
errors,  should  be  within  ±1/2%  of  range  or  better. 

5.2  Control  Position  Transmitters 

Control  position  transmitters,  commonly  referred  to  as  CPT  units,  sense  the  control- 
surface  deflections  and  must  be  accurate  and  sensitive  enough  to  measure  small 
deflections.  Transmitters  of  the  sliding  contactor  type  change  the  ratio  of  resistance 
in  two  arms  of  a  Wheatstone  bridge  circuit.  Any  variation  in  the  resistance  of  the 
arms  unbalances  the  circuit  and  causes  current  to  flow  to  the  recording  galvanometer 
(see  Figure  24). 

In  a  properly  installed  system,  the  phase  lag  between  the  transmitter  and  the  re¬ 
corder  should  be  negligible.  The  errors  due  to  hysteresis,  zero  shift,  temperature, 
accelerations,  or  vibrations  should  also  be  negligible. 

The  transmitters  are  firmly  mounted  at  the  contro]  surfaces  to  eliminate  the  effect 
of  control-system  deformations.  The  spanwise  location  of  the  transmitter  gives  an 
approximate  spanwise  surface  deflection. 

Zero  checks  are  made  before  and  after  each  flight  to  detect  any  zero  shift  in  the 
galvanometer  recording  system. 

5.3  Angle-of- Attack  and  Sideslip 
5.3.1  Vane-Type  Flow-Direction  Sensors 

Of  the  various  types  of  flow-direction  devices  for  sensing  angle  of  attack  and 
sideslip  up  to  a  Mach  number  of  approximately  3.0,  good  accuracy  and  reliability  is 
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obtained  with  a  counterbalanced,  freely  turning  vane  mounted  on  a  nose  boom,  which 
also  serves  as  a  mount  for  the  pitot-static  head  (Fig. 19).  Each  vane  is  directly 
connected  to  a  synchro  transmitter  within  the  boom,  which  !e  electrically  connected 
to  a  synchro  receiver  mounted  in  a  recorder  located  withir.  the  airplane.  It  should 
be  noted  that,  although  the  a-vane  measures  the  aerodynamic  a.  ,  the  /3-vane  measures 
a  p  referenced  to  the  body  axis  system  of  the  airplane. 

Inherent  accuracy:  Hysteresis  in  the  system  is  practically  nil.  Friction  introduces 
an  error  of  less  than  ±0. 1°.  In  an  optical  recorder  system,  the  unbalance  of  a  balanced 
optical  recorder  element  may  cause  a  trace  deflection  equivalent  to  0.05°  per  g  of 
acceleration.  Temperature  has  no  direct  effect  on  sensitivity.  Natural  frequency  and 
damping  of  the  system  should  be  of  the  order  to  10  c/s  and  0.65,  respectively,  to  pro¬ 
vide  flat  response  to  within  ±1%  for  sinusoidal  inputs  up  to  6  c/s. 

Mounting:  Hie  angle -of -attack  and  angle -of -sideslip  vanes  are  mounted  on  a  nose  boom 
extending  forward  as  far  as  possible  to  minimise  the  effects  of  upwash  and  shock  wave. 

In  this  respect,  vanes  are  mounted  1%  maximum  fuselage  diameters  ahead  of  the  airplane 
when  feasible.  Figure  25,  reproduced  from  Reference  30,  shows  the  theoretical  effects 
of  upwash  from  the  nose  boom  and  fuselage  at  low  speeds.  Wing  upwash  was  not  considered. 

Hie  boom  and  its  mount  should  be  sufficiently  stiff  to  minimize  deflections  due  to 
inertia  and  air  loads.  Particular  care  must  be  taken  to  aline  the  longitudinal  axis 
of  the  boom  with  the  longitudinal  body  axis  of  the  airplane  and  the  vane  struts  to  the 
boom  so  that  the  angle-of-attack  and  angle-of-sideslip  vane  struts  are  parallel  to  the 
body  y  and  z  axes,  respectively.  Hie  rear  vane  is  for  sideslip  and  projects  verti¬ 
cally  downward.  Hie  recorder  is  mounted  in  any  convenient  location. 

Field  checks:  Hie  final  calibration  of  a  transmitter-recorder  combination  is  made  in 
place  on  the  airplane  with  the  aid  of  a  calibration  fixture  that  provides  an  accurate 
alinement  of  the  vane  with  the  boom  and  the  zero  of  the  calibration  quadrant.  Calibrations 
should  be  made  in  increments  of  about  2°  to  detect  nonlinearities.  Calibration  should 
be  performed  both  before  and  after  flight. 

Hie  vanes  should  be  given  periodic  checks  for  alinement  with  their  pivotal  shafts 
and  for  friction.  An  extension  of  the  chordline  of  the  vane  should  be  within  0.005  in. 
of  alinement  with  the  center  of  the  shaft. 

Correction  of  recorded  data:  The  angle-of-attack  and  angle-of-sideslip  vanes  measure 
local  flow  direction.  The  effects  of  boom  bending  due  to  inertia  and  air  loads,  flow 
components  resulting  from  angular  velocities,  flight-path  curvature,  and  upwash  due  to 
the  boom,  fuselage,  and  wings  introduce  errors  in  the  measured  flow  angles  with  respect 
to  the  true  airplane  angle-of-attack  or  sideslip.  In  addition,  phase  lag  and  dynamic 
amplification  of  the  sensing -recording  system  introduce  additional  errors  in  the  re¬ 
cording  of  the  vane  indications.  The  magnitude  of  each  effect  must  be  investigated 
and  corrections  made  to  the  recorded  data  wherever  pertinent  to  the  analysis  for  deter¬ 
mination  of  derivatives. 

Bending  of  the  boom  results  in  errors  in  vane  indications,  inasmuch  as  the  vane  is 
referenced  to  the  axis  of  the  boom.  As  pointed  out  in  Reference  30,  deflections  due  to 
aerodynamic  loading  have  been  negligible;  however,  where  very  long  booms  are  used  or 
extremely  long,  flexible  fuselages  are  being  dealt  with,  bending  corrections  may  be 
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determined  from  calculated  aerodynamic  loading31.  Boom-bending  error  restilting  "on 
inertia  loads  is  accounted  for  through  static  deflection  calibration  of  the  boon. 

Upwash  error  resulting  from  the  boom,  fuselage,  and  ving  is  generally  considered 
negligible  and  within  the  accuracy  of  the  methods  of  analysis  employed.  This  may  not 
necessarily  be  true.  In  a  boom-vane  installation  on  a  large  bomber  where  the  vane  was 
one  fuselage  diameter  ahead  of  the  nose,  the  upwash  error  in  angle-of-attack  at  sub¬ 
sonic  speeds  was  cf  the  order  of  4%.  On  other  large  aircraft  having  fuselages  of 
larger  cross  section,  the  Influence  was  much  larger.  Upwash  error  due  to  the  boom 
itself  may  be  measured  by  wind-tunnel  calibration  of  the  system28’34.  The  effect  of 
upwash  at  subsonic  speeds  r.t  the  vane  due  to  the  fuselage  can  be  calculated  by  the 
method  of  Reference  35.  The  effect  of  upwash  at  subsonic  speeds  due  to  the  wing  can 
be  calculated  by  the  equations  in  Reference  36  for  unswept  wings  and  the  methods  of 
Reference  37  for  swept  wings.  At  supersonic  speeds,  the  wing  and  fuselage  do  not 
contribute  any  upwash  effects  to  the  vane. 


The  angle-of-attack  sensor  is  also  subject  to  pitch-rate  effects  of  flight-path 
curvature.  Corrections  for  flight-path  curvature  (Pig. 26)  may  be  significant  at  sub¬ 
sonic  speeds;  whereas,  pitch -rate  corrections  may  be  significant  from  the  subsonic 
through  the  low  supersonic  speed  range.  Corrections  for  flight-path  curvature  affect 
magnitude  primarily;  whereas,  corrections  for  pitch  rate  affect  phase  angle  primarily. 
This  is  illustrated  in  Figure  27,  which  shows  the  graphical  time-vector  determination 
of  the  absolute  amplitude  of  the  corrected  angle  of  attack  as  a  ratio  of  the  indicated 
amplitude  for  an  aircraft  performing  small-perturbation,  free-osciilation  maneuvers  at 
a  Mach  number  of  0.8  at  40,000  ft.  The  presentation  considers  only  corrections  for 
flight-path  curvature  and  pitch-rate  affects  and  is  based  on  the  equation 
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The  solution  shows  the  Influence  of  the  flight -path  curvature  to  be  of  the  order  of  3%. 


Plight-path  curvature  in  yaw  has  a  negligible  effect  on  the  sideslip  vane.  Correction 
for  yaw-rate  and  roll -rate  effects  should  be  considered.  The  approximate  expression  for 
correcting  the  indicated  sideslip  for  angular-rate  effects  is 
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5.3.2  Spherical  Hypersonic  Flovi-Direction  Sensor 

The  spherical  flow-direction  sensor  shown  in  Figure  28  was  designed  to  replace  the 
a  and  /3  vane-type  sensors  at  the  higher  supersonic  Mach  numbers  and  dynamic  pressure 
where  the  combined  temperature  and  aerodynamic  loads  exceed  the  limitations  of  the  vane 
type  sensor38.  The  spherical  sensor  is  a  null -seeking,  hydraulically  operated,  electro' 
nlcally  controlled  servo-mechanism.  It  has  pressure  measurements  as  its  sole  sensing 
inputs.  It  operates  on  the  principle  that  when  two  static  ports  are  located  on  the 
great  circle  of  a  sphere,  a  null  reading  will  result  when  the  bisector  of  the  Included 
angle  of  the  .‘■■wo  static  ports  is  parallel  to  the  fluid  stream  Immediately  in  front  of 
the  sensor.  The  rotation  of  the  bisecting  line  relative  to  a  reference  gives  the 
inclination  of  the  fluid  stream  relative  to  the  reference. 
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When  the  spherical  sensor  is  in  the  zero  position  (axis  alined  with  the  airplane), 
the  a-ports  are  42°  above  and  below  the  reference  line  in  the  vertical  plane  of  symmetry 
of  the  aircraft  and  the  ^J-ports  are  42°  on  either  side  of  the  reference  line  in  the 
transverse  plane. 

The  sphere  constitutes  the  outer  gimbal  of  a  two-gimbal  pivot  syst'ia  in  which  the 
outer  gimbal  is  pivoted  to  the  inner  gimbal  whose  pivotal  axis  is  fixed  and  is  normal 
to  the  plane  of  symmetry  of  the  airplane.  As  the  sensing  sphere  seeks  null  readings 
in  each  of  its  two  sets  of  static -pressure  ports,  the  gimbals  rotate  about  their  res¬ 
pective  axes.  The  inner  gimbal,  rotating  about  its  fixed  axis,  which  is  normal  to  the 
plane  of  symmetry,  sweeps  an  angle  a  in  the  plane  of  symmetry.  The  outer  gimbal, 
whose  pivotal  axis  is  mounted  on  the  inner  gimbal  and  remains  in  the  plane  of  symmetry 
at  all  times,  sweeps  an  a*igle  /3  in  a  plane  which  is  perpendicular  to  the  plane  of 
Symmetry;  this  plane  is  tL*i  transverse  plane  of  the  stability  axis  system  of  the 
aircraft.  The  a  and  /3  angles  nicked  off  by  synchros  are  the  aerodynamic  a  and 
/3  angles  of  the  airplane. 

The  inherent  accuracy  of  the  spherical  sensor  is  of  the  order  of  ±0.5°  or  better 
for  dynamic  pressures  in  excess  of  20  lb/ft2.  At  high  angles  of  attack  in  excess  of 
approximately  26°  at  low  dynamic  pressures  of  about  40  lb/ft2  and  less,  the  a  indica¬ 
tions  are  subject  to  large  errors,  possibly  due  to  flow  interference  of  the  lip  on  the 
collar  of  the  housing. 

5.4  Angular  Velocities  and  Accelerations 

The  angular  velocity  and  angular  acceleration  relative  to  any  one  axis  can  be  sensed 
by  individual  sensors  or  sensed  and  recorded  in  a  convenient  packaged  unit.  Figures 
39(a)  and  39(b)  show  the  details  of  the  angular-velocity  aspect  of  a  NACA  designed, 
packaged  unit  which  includes  the  recorder.  The  angular  acceleration  sensing  and  re¬ 
cording  involves  a  relatively  small  extension  of  this  unit.  The  operation  of  the  unit 
depends  upon  the  precessional  force  of  a  restrained  gyro  motor  when  the  unit  is  sub¬ 
jected  to  an  angular  rate  about  an  axis  which  is  perpendicular  to  both  the  axis  of 
rotation  of  the  gyro  motor  and  the  axis  of  rotation  of  the  gimbal  rings.  The  gyro¬ 
scopic  element  is  the  rotor  of  a  synchronous  motor.  The  sensitive  element  is  res¬ 
trained  hy  a  precision  helical  spring.  The  moving  system  is  damped  by  rotating  an 
aluminium  disk  In  the  field  of  a  strong  permanent  magnet.  The  angular -velocity  measure¬ 
ment  is  made  by  optically  recording  on  the  film  the  angular  displacement  of  the  gimbal. 
Sensitivity  of  the  angular-velocity  recorder  can  be  adjusted  by  rotating  the  actuator 
arm  along  the  mirror  staff  tail. 

Angular  acceleration  is  obtained  by  differentiating  the  gimbal  motion.  The  differ¬ 
entiation  is  accomplished  by  mounting  a  coil  in  a  magnetic  field  ard  driving  it  from 
the  damping  shaft  so  that  it  rotates  with  speed  proportional  to  the  angular  velocity 
of  the  gimbal.  The  output  voltage,  which  is  proportional  to  the  angular  acceleration, 
is  recorded  on  the  film  by  a  self-contained  reflecting  galvanometer. 

5.4.1  Inherent  Accuracy 

In  well -designed  angular-velocity  systems,  the  reading  accuracy  is  of  the  order  of 
0. 5%  of  full  scale  or  better;  the  errors  due  to  friction  and  hysteresis  are  less  than 
1%  of  full  scale,  and  the  change  in  sensitivity  from  large  changes  in  temperature  should 
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be  as  small  as  possible.  Errors  due  to  linear  accelerations  of  5g  should  be  less 
than  1%.  The  sensor  should  provide  flat  response  characteristics  within  ±1%  for  all 
anticipated  impressed  frequencies.  The  phase  lag  (time  lag)  is  a  function  of  damping 
ratio  and  undamped  natural  frequency  of  the  sensor. 

Recorded  angular  accelerations  are  subject  to  the  errors  found  in  the  angular- 
velocity  record.  In  the  NACA  acceleration-velocity  packaged  unit,  additional  errors 
are  introduced  by  the  acceleration  recording  galvanometer;  inasmuch  as  the  angular- 
acceleration  pickup  is  a  differentiation  device,  the  response  and  phase  lag  of  the 
accelerometer  and  velocity  portions  of  the  unit  are  similar. 

5.4.2  Mounting  and  Corrections 

It  is  iuportant  that  the  instrument  mounting  be  rigid.  Although  small -amplitude, 
high-frequenqy  vibrations  may  not  be  apparent  on  the  velocity  trace,  the  vibrations 
can  introduce  considerable  noise  in  the  acceleration  trace. 

Angular-velocity  gyros  are  subject  to  coupling  errors  caused  by  an  interference 
(airplane)  angular  velocity  about  the  spin  axis  of  the  gyro  rotor.  Care  should  be 
exercised  in  orienting  the  instrument  during  mounting  so  as  to  subject  its  spin  axis 
to  the  minimum  interference  angular  velocity.  A  mathematical  study  of  the  coupling 
error  is  presented  in  Reference  39.  Tan  interference  angular  velocity  (also  known 
as  the  q  rate)  affects  the  sensitivity  of  the  instrument,  the  undamped  natural  fre¬ 
quency,  and  the  damping  ratio.  The  extent  of  the  errors  is  a  function  of  the  glmbal 
tilt,  which,  itself,  is  a  function  of  the  gyro  sensitivity  in  spring -restrained  in¬ 
strument  and  the  magnitude  of  '.ue  interference  angular  velocity.  This  is  illustrated 
in  Figure  40  for  an  angular-ve  ocity  unit  having  a  static  sensitivity  of  0.256  radian 
per  radian  per  second.  A  <  .-crease  in  sensitivity  would  reduce  the  coupling  error; 
however,  a  decrease  is  not  always  desirable.  To  minimize  the  coupling  error  for  any 
one  instrument,  the  axes  should  be  oriented  as  follows; 
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Alinement  of  the  sensing-recording  units  should  be  within  ±0.2°  of  correct  orientation 
with  relation  to  the  body  axes.  Undetected  misalinement  has  been  known  to  result  in 
erroneous  values  of  highly  pertinent  derivatives,  which  resulted  in  misleading  results 
in  analog-simulated  rolling  characteristics.  In  any  correction  for  misalinement,  it 
is  pertinent  that  the  recorded  values  be  corrected  for  phase  lag  of  the  instrument  prior 
to  insertion  in  the  correction  equations.  Simultaneously,  the  response  of  the  instru¬ 
ment  should  be  checked,  if  there  is  any  appreciable  deviation  from  the  damping  ratio 
of  0.65,  to  ascertain  the  percentage  error  in  magnitude  of  the  indicated  quantity  due 
to  the  dynamics  of  the  instrument.  Mlsalinements  in  the  mounting  of  the  unit  may  be 
accounted  for  by  using  the  equations  shown  in  Figure  31. 


5.5  Linear  Accelerations 

In  general,  flight  testing  is  done  with  the  bean-type  linear  accelerometers  which 
are  available  as  single-component  or  three -component  units.  Drag  determination  is 
frequently  mode  with  single-component  units30.  The  beam-motion  restraining  force  is 
generally  supplied  by  a  pair  of  opposed  helical  springs.  The  sensitivity  and  undamped 
natural  frequency  are  dependent  upon  the  springs  used. 

5.5.1  Inherent  Accuracy 

In  properly  designed  beam-type  linear  accelerometers,  sensitivity  and  zero  changes 
from  random  causes  are  less  than  0.  5%  of  full  scale.  The  sensor  should  have  a  damping 
ratio  of  0.65  and  a  sufficiently  high  undamped  natural  frequency  to  provide  flat  res¬ 
ponse  characteristics  within  ±1%  for  impressed  frequencies  up  to  60%  of  the  undamped 
natural  frequency  of  the  sensor.  Each  linear  accelerometer  is  affected  by  an  inter¬ 
acting  acceleration  acting  along  the  beam.  The  effect  is  generally  small  but  should 
net  be  arbitrarily  ignored. 

5.5.2  Mounting  and  Corrections 

The  instrument  should  be  moulted  as  close  to  the  center  of  gravity  of  the  airplane 
as  possible.  It  should  be  rigidly  fastened  on  a  rigid  mounting  attached  to  the  primary 
structure  of  the  airplane  to  avoid  or  at  least  minimize  extraneous  vibratory  accelera¬ 
tions.  It  should  be  alined  to  within  ±0.2°  of  correct  orientation  with  relation  to 
all  three  reference  axes.  When  the  instrument  is  not  mounted  at  the  center  of  gravity 
of  the  airplane,  corrections  of  the  indicated  readings  to  the  center  of  gravity  must 
be  made  by  using  the  expressions  shown  in  Figure  32.  The  equations  for  normal  accelera¬ 
tion.  an  ,  and  transverse  acceleration,  at  ,  can  be  linearized  and  corrections  thus 
simplified  by  mounting  the  accelerometers  in  the  plane  of  symmetry  along  the  x-axis. 

5.6  Phase  Lag  and  Response 

Since  several  individually  recorded  quantities  are  utilized  in  the  determination  of 
various  derivatives,  it  is  important  that  the  phase-lag  (time-lag)  characteristics  of 
each  recording  instrument  be  taken  into  consideration.  For  systems  where  all  the 
quantities  can  be  nscorded  on  electrical  galvanometers,  it  is  generally  possible  to 
equalize  the  individual  phase  lags  by  proper  choice  of  the  frequency  response  of  the 
recording  system.  Where  this  is  not  possible,  as  in  the  use  of  certain  cf  the  self- 
recording  NASA  Instruments,  phase-lag  corrections  must  be  considered  and  applied  to 
bring  all  pertinent  quantities  into  correct  time  relationship. 

Phase-lag  corrections  must  be  applied  before  making  any  corrections  for  misalinement. 
Corrections  for  misalinement  must  be  made  before  correcting  the  vane  and  linear -accelero¬ 
meter  records  to  the  center  of  gravity  of  the  airplane. 

Because  of  the  nature  of  the  control  inputs,  phase-lag  corrections  can  be  applied 
simply  by  shifting  the  data  time  scale40,  as  in  the  determination  of  control  derivatives, 
or  by  correcting  phase-angle  relationships,  as  in  the  time-vector  method  of  analysis. 

This  is  accomplished  by  determining  the  undamped  natural  frequency  of  the  airplane 
from  free-oscillation  maneuvers  and  Figure  33.  Amplitude  corrections  are  not  required, 
since  the  instruments  have  flat  response  characteristics. 
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When  the  instruments  are  not  sufficiently  damped  to  provide  flat  response  character¬ 
istics,  corrections  to  the  magnitudes  of  the  recorded  quantities  may  be  determined 
from  Figure  34. 

5.7  Ranges  and  Sensitivity 

Instruments  used  for  studies  of  general  handling  qualities  have  relatively  low 
sensitivities  in  order  to  accommodate  the  normal  flight  range  and  are  used  for  approxi¬ 
mate  evaluation  of  derivatives  in  conjunction  with  these  studies.  For  accurate  evalua¬ 
tion  of  the  derivatives,  using  small  disturbance  maneuvers,  sensitive  gyros  anu  accelero¬ 
meters  are  installed  to  supplement  or  replace  those  used  for  the  handling-qualities 
studies.  The  ranges  and  sensitivities  of  the  instruments  are  usually  selected  after 
studying  flight  test  records  of  small -perturbation  maneuvers  performed  over  a  Mach 
number  range  during  pilot  familiarization  flights  when  the  airplane  is  equipped  with 
general-purpose  flight  test  instruments.  The  increase  in  sensitivity  of  any  one  in¬ 
strument  must  be  accomplished  with  discretion,  inasmuch  as  an  optimum  sensitivity  is 
attained  beyond  which  any  increase  may  result  simply  in  a  false  sense  of  accuracy. 

Table  VII  shows  the  characteristics  of  instruments  which  are  desirable  for  derivative 
investigations  for  one  high-performance  airplane  when  the  pulsed  free-oscillation 
maneuver  is  employed.  The  listed  instrument  natural  frequencies  are  more  than  adequate 
to  maintain  flat  response  characteristics  during  forced  portions  of  the  maneuver  up  to 
the  anticipated  maximum  frequencies  for  all  recorded  quantities. 

5.8  Pulse  Code  Modulation  (PCM)  Data-Acquisition  Systems 

In  the  preceding  considerations  of  instrumentation,  emphasis  was  placed  on  factors 
that  affect  the  accuracy  of  individual  sensors.  Self-contained  sensor-recorded  units 
are  compact,  reliable,  and  accurate.  The  use  of  sensors  wired  to  remote  recorders 
can  introduce  degradation  in  the  accuracy  of  the  overall  sensor-recorder  system;  how¬ 
ever,  such  systems  are  used  to  keep  the  instrumentation  volume  to  a  minimum  where 
space  is  a  prime  factor  and  a  large  number  of  parameters  are  involved.  As  the  number 
of  sensed  and  recorded  parameters  increases,  the  time  lag  in  the  recovery  of  the  data 
for  the  user  increases.  In  flight  test  investigations  where  the  bulk  of  the  instru¬ 
mentation  is  a  serious  problem  or  where  the  number  of  parameters  recorded  may  constitute 
a  serious  time  lag  in  the  recovery  of  the  data  for  the  user,  a  sophisticated  data- 
acquisition  system  is  available  to  alleviate  these  problems.  This  system,  originated 
to  fulfill  the  needs  of  the  space  industry,  in  which  transducers  of  superior  quality 
are  used,  is  capable  of  handling  the  data  to  reasonable  accuracy  (0.2%  to  1%).  The 
system,  referred  to  as  the  PCM  system,  converts  the  analog  signal  from  the  sensor  to 
digital  format  and  records  the  digitized  data  on  tape  on  a  time-sharing  basis. 

Figure  35(a)  shows  a  schematic  drawing  of  an  airborne  PCM  system.  The  analog  signals 
from  the  sensors  go  to  a  PCM  encoder  to  convert  the  signal  to  an  identification  coded, 
digitized  format.  The  coded,  digitized  signals  are  then  recorded  in  parallel  on  an 
onboard  tape  recorder  on  a  time-sharing  basis.  To  recover  the  data,  the  taped  signals 
are  processed  through  a  PCM  decommutation,  which  identifies  (unscrambles)  the  individual 
ser.sor  signals,  to  a  format  computer  to  provide  real-time  data  outputs  in  the  form  of 
strip  charts  or  oscillograph  readouts  for  an  immediate  look  at  the  data.  The  real-time 
data  are  also  transmitted  to  a  general-purpose  computer  which  tabulates,  plots,  or 
performs  complex  manipulation  of  the  data  in  engineering  units. 


46 


Where  weight  is  a  serious  factor,  Figure  35(b)  shows  a  schematic  drawing  of  the 
PCM  system  using  telemetry.  The  main  differences  between  the  telemetered  and  airborne 
PCU  systems  involves  the  transmission  of  the  coded,  digitized  signals  in  series  to  the 
decommutator  (instead  of  parallel  to  the  recorder),  which  provides  time  synchronization 
of  the  signals  before  the  signals  are  taped.  In  processing  the  data,  the  format  com¬ 
puter  properly  identifies  the  individual  data  channels  for  real-time  data  out;nit. 

As  stated  earlier,  the  PCM  system  is  a  sophisticated  operation.  One  Installation  at 
the  NASA  Flight  Research  Center,  Edw  ds,  California,  is  designed  to  handle  15,400  data 
samples  per  second  from  77  to  a  maximum  of  800  data  sourcee 


6.  FLIGHT  TEST  TECHNIQUES 

Determination  of  the  flight  test  techniques  to  be  used  in  obtaining  stability  and 
control  derivatives  from  flight  data  is  governed  by  a  number  of  factors,  including  the 
methods  of  analysis  to  be  employed.  Successful  mathematical  methods  of  analysis  have 
been  limited  to  the  linearized  form  of  the  equations  of  motion  and  thus  restrict  the 
maneuvers  to  small  perturbations.  Inasmuch  as  stability  derivatives  are  functions  of 
angle  of  attack  and  Mach  number  and,  to  some  extent,  aeroelBsticity  of  the  airframe, 
the  controlled  variables  are  Mach  number,  load  factor,  and  pressure  altitude.  For 
safety  of  flight,  the  investigation  of  the  stability  and  control  characteristics  is 
usually  initiated  with  a  gradual  buildup  of  maneuvers  at  high  altitude  where  the 
natural  frequency  and  damping  of  the  airplane  are  lower  than  at  low  altitudes  and 
thus  permit  better  control.  It  is  desirable,  when  feasible,  to  have  the  maneuvers 
performed  with  the  airplane  weight  within  such  limits  over  the  derivative-determination 
phase  of  the  flight  test  program  that  the  effects  of  changes  in  centei -of -gravity 
position  and  moments  of  inertia  will  be  negligible. 

The  important  factors  to  be  considered  in  flight  testing  for  stability  and  control 
derivatives  are  discussed  in  the  following  sections. 

6. 1  Mach  Number  and  Altitude 

Flight  test  maneuvers  are  generally  performed  at  lg  initial  conditions  at  constant 
Mach  number  and  altitude.  Normally,  some  variations  in  these  quantities  are  accepted 
if  the  resultant  change  in  dynamic  pressure  is  not  more  than  5%  over  that  portion  of 
the  maneuver  encompassed  in  the  analysis.  In  regions  where  large  Mach  number  effects 
exist  (Fig. 36),  tests  should  be  conducted  at  close  Mach  number  intervals  with  more 
rigid  requirements  at  constant  Mach  number  and  altitude.  Failure  to  trim  the  aircraft 
to  the  desired  Mach  number  and  to  nu-  ;rf  a  that  Mach  number  during  the  maneuver  in 
regions  of  rapidly  varying  characteristics  may  produce  a  scatter  of  data  and  an 
erroneous  analysis. 

The  very  nature  of  flight  testing  requires,  for  expediency,  plotting  the  results 
of  analysis  as  a  function  of  Mach  number,  with  each  curve  representing  a  constant- 
altitude  condition.  Figure  37,  taken  from  Reference  41,  shows  the  influence  of 
altitude  on  flight  test  data  on  one  supersonic  aircraft. 


\ 
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6.2  Angle  of  Attack  and  Load  Factor 

Hie  variation  in  airplane  characteristics  with  angle  of  attack  is  determined  by 
performing  maneuvers  at  different  altitudes  with  lg  trim  conditions  prevailing  prior 
to  tbe  perturbation,  or  at  constant  altitude  witb  tbe  maneuver  performed  during  a 
stabilized  constant -g  pushover  or  turn.  It  should  be  noted  tbat  it  is  difficult  to 
obtair  good  maneuvers  during  stabilized  turns;  exceptional  piloting  skill  is  required. 
Figures  37  and  38  sbow  the  influence  of  load  factor  on  stability  characteristics.  In 
instances  where  tbe  aeroelasticity  of  the  structure  is  nil  (dynamic  pressure  effects 
are  nil),  a  combination  of  tbe  two  techniques  will  result  in  the  determination  of  the 
variation  of  tbe  derivatives  over  an  extended  range  of  angle  of  attack.  Should  aero¬ 
elasticity  of  the  structure  be  a  factor  to  contend  with,  the  results  from  tbe  two 
techniques  will  differ  for  tbe  same  angle  of  attack,  Mach  number,  and  center  of  gravity. 

6. 3  Aeroelasticity 

Aeroelastic  deformation  of  the  structure  assumes  Increasing  significance  as  tbe 
aircraft  increases  in  size  and  slenderness  and  operates  at  increasing  dynamic  pressures. 
Supersonic  transport  designs  are  flexible  in  order  to  keep  the  structural  weight  down, 
tbe  payload  bigh,  and  the  range  capability  a  maximum.  To  apply  theoretical  flexibility 
corrections  to  rigid  wind-tunnel  data  for  comparisons  with  flight  data  provides  an 
intuitive  basis  in  ascertaining  flexibility  effects.  When  such  comparisons  are  em¬ 
ployed  and  a  definite  disagreement  is  evident  in  the  comparison  in  regard  to  level  and 
trends  of  tbe  stability  and  control  parameters  as  a  function  of  Mach  number,  it  may 
become  difficult  to  locate  the  source  of  tbe  discrepancy  -  wind-tunnel  data  or  predicted 
flexibility  corrections.  Thus,  a  more  positive  approach  is  required  to  assess  flexibi¬ 
lity  effects. 

The  stability  and  control  derivatives  should  be  essentially  invariant  for  a  rigid 
airplane  as  long  as  Mach  number,  angle-of -attack,  and  tbe  center  of  gravity  are  constant 
(assuming  Reynolds  number  effects  to  be  a  minor  factor).  Thus,  any  direct  approach  to 
investigating  flexibility  effects  based  on  flight  data  should  show  the  variation  of 
the  stability  parameters  -  obtained  at  the  same  Mach  number,  angle-of-attack,  and  center 
of  gravity  -  as  a  function  of  dynamic  pressure.  Although  Mach  number  and  center-of- 
gravity  control  is  straightforward,  the  angle-of-attack  is  a  problem. 

The  location  of  the  angle-of-attack  sensor  exposes  the  sensor  to  errors  resulting 
from  structural  deformations,  in  addition  to  the  other  sources  discussed  in  Section  5. 3. 
Hence,  it  is  more  judicious  to  use  the  life  coefficient  CL  in  lieu  of  angle-of-attack 
a  .  Thus,  from  a  practical  point  of  view,  a  direct  Investigation  of  aeroelastic 
effects  should  be  based  on  a  comparison  of  flight  data  for  different  dynamic-pressure 
conditions  obtained  at  the  same  Mach  number,  lift  coefficient,  and  center  of  gravity. 

An  effective,  flexible,  and  simple  flight-planning  procedure  to  determine  the  flight 
test  conditions  as  a  function  of  weight  and  altitude  to  provide  constant  M  ,  CL  ,  and 
center  of  gravity  can  be  achieved  by  using  a  nomograph  such  as  that  in  Figure  39.  In 
this  nomograph  W  ,  M  ,  CL  ,  and  q  are  variables,  and  center  of  gravity  is  constant. 
The  nomograph  is  based  on  the  following  two  basic  relations  for  lg  flight; 

ff  =  CLqS  (110) 

q  =  0. 7pMz  . 


and 


(111) 
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It  assuaes  the  weight  distribution,  which  could  influence  structural  deformation,  to 
be  essentially  constant.  Inasmuch  as  CL  is  a  constant  for  any  one  Mach  number  con¬ 
dition  being  investigated,  the  following  expression  is  readily  derived  from  the  above 
equations  and  constitutes  the  basis  for  the  nomograph: 


(112) 


The  subscripts  1  and  2  denote  the  initial  and  compatible  second  condition.  It  will 
be  noticed  that,  for  any  one  initial  weight  Wj  at  altitude  hj  (as  typified  by 
pressure  p'j),  the  vehicle  will  have  to  be  at  a  weight  W2  at  altitude  h2  to  maintain 
the  same  CL  at  the  selected  constant  Mach  number. 


To  illustrate  the  use  of  the  nomograph,  consider  an  aircraft  to  have  a  weight  of 
411  x  103  lb  at  the  time  a  stability  maneuver  was  performed  at  Mach  2.34  at  55  x  103  ft. 
These  initial  conditions,  which  nave  been  shotted  on  Figure  39,  show  the  dynamic  pressure 
to  be  730  lb/ft2.  If  it  is  desired  to  perform  the  next  stability  maneuver  at 


q2  =  450  lb/ft2  ,  the  intersection  of  q2  (450)  and  the  constant  Mach  line  (2.34) 
determines  the  new  altitude,  h2  ,  to  be  65  x  103  ft.  The  intersection  of  the  constant- 


altitude  line  with  the  constant 


C,  ,  center-of -gravity  line  extended  from  condition 
1  determines  the  weight  (W2  =  252  x  103  lb)  required  to  provide  the  same  lit  and  CL 


at  condition  2  as  was  present  at  the  time  of  the  stability  maneuver  at  condition  1 


(center-of-gravity  being  constant). 


The  nomograph  is  invaluable  in  systematic  flight  planning  for  determinat  on  of  aero- 
elastic  effects.  It  permits  on-the-spot  changes  in  planned  flight  conditions.  It  also 
accentuates  the  large  changes  in  weight  required  to  obtain  significant  changes  in 
dynamic  pressure  to  assure  aeroelastic  flight  data  which  will  be  outside  the  area  of 
experimental  error  of  uncertainty. 


6.4  Control  Inputs 

Die  method  of  analysis  selected  governs  the  control  input.  The  magnitude  and  duration 
of  the  input  influence  the  magnitude  of  the  perturbation.  In  the  case  of  an  aerodynamic 
coefficient  that  is  highly  nonlinear  with  respect  to  an  independent  variable,  different 
magnitudes  of  the  perturbation  may  result  in  different  magnitudes  of  the  derivative  of 
the  coefficient  in  analyzing  flight  data.  Thus,  in  comparing  flight  results  with  wind- 
tunnel  data,  it  is  essential  that  the  wind-tunnel  value  of  the  derivative  be  based  not 
only  on  the  same  trim  condition  but  also  on  the  same  magnitude  of  perturbation  as  the 
flight  data. 

Where  nonlinearity  of  the  coefficients  is  not  a  factor  and,  in  lieu  of  increase  of 
instrument  scale  factor,  larger  perturbations  of  the  independent  variables  are  desired 
to  provide  more  accurate  readability  of  the  records,  larger  control  inputs  or  co^nlex 
control  inputs  may  be  used.  Figure  40  shows  the  increase  in  amplitudes  of  recirded 
quantities  resulting  from  a  change  in  control  input. 


6. 5  Maneuvers 

Maneuvers  performed  for  determination  of  stability  and  control  derivati,  ■?  from 
flight  data  should  be  compatible  with  the  requirements  of  the  method  of  analysis  to  be 
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employed.  Current  practical  methods  of  analysis,  whether  they  involve  approximate 
equations  solving  for  individual  derivatives  or  comprehensive  techniques  solving  a 
number  of  derivatives,  have  limitations  in  their  utility;  as  a  result,  different  types 
of  maneuvers  are  employed  within  the  range  of  their  individual  limitations  to  obtain 
the  derivatives.  As  a  generality,  it  might  be  said  that  typical  handling-quality 
maneuvers  are  employed  in  the  determination  of  derivatives  wherein  analytical  techniques 
are  used.  Included  are  longitudinal  elevator-pulse  maneuvers,  pullups  and  push-overs, 
pullups  and  releases,  rudder-pulse  and  aileron-pulse  maneuvers,  constant -heading  side¬ 
slips,  recovery  from  sideslip,  and  rudder-fixed  rolls. 

When  flight  maneuvers  applicable  tc  analytical  technique  for  derivative  determination 
are  not  available  or  usable,  the  airplane  response  to  random  inputs  is  analyzed  to  give 
limited  stability  data.  This  is  accomplished  effectively  with  the  aid  of  an  analog 
computer,  using  a  technique  involving  the  matching  of  analog  and  flight  time  histories. 

6.5.1  Pulse  Maneuvers 

The  simple  pulse  maneuver,  shown  in  Figure  41  for  a  longitudinal  perturbation,  is 
the  current  mainstay  for  derivative  determination.  Normally,  for  this  maneuver  the 
airplane  is  trimmed  at  the  desired  angle-of -attack,  altitude,  and  Mach  number,  and  a 
free  oscillation  is  initiated  by  an  abrupt  pulse  -  an  elevator  pulse  for  longitudinal 
oscillation,  a  rudder  or  aileron  pulse  for  lateral -directional  oscillations.  The 
resulting  free-oscillation  of  the  aircraft  is  allowed  to  damp  out  with  the  controls 
held  fixed  at  the  initial  trim  setting.  With  an  irreversible  control  system,  this  is 
easily  accomplished  by  releasing  the  controls.  On  tailless  aircraft,  even  small  in¬ 
advertent  control  inputs  during  the  free  oc ;i llat ion  can  significantly  affect  the 
damping  and,  hence,  the  damping  derivatives.  Moderate  inadvertent  control  inputs  can 
affect  the  period  of  oscillation,  as  well  as  the  damping,  and  then  influence  the  static 
derivative  results  as  well. 

Free  oscillations  are  also  initiated  by  release  of  controls  at  the  end  of  a  side¬ 
slip  maneuver  and  at  the  end  of  pullup  and  push-over  maneuvers. 

In  investigating  the  effects  of  angle  oi  attack  and  load  factor  when  utilizing  the 
pulse  maneuver  in  an  elevated  g  turn,  the  application  of  the  pulse  technique  is 
limited  by  the  difficulty  of  performing  a  good  maneuver.  Difficulty  has  been  exper¬ 
ienced  during  the  maneuver  in  holding  the  proper  bank  angle  to  maintain  constant  load 
factor  and  Mach  number.  With  a  conventional  control  system,  exceptional  piloting 
skill  is  required  to  maintain  fixed  control  during  the  airplane  oscillations  at 
elevated  g  .  The  use  of  the  airplane  damper  as  a  device  for  applying  a  known  deflec¬ 
tion  signal  to  excite  the  desired  unaugmented  oscillations  (Fig.  42)  offers  a  means  of 
improving  the  quality  of  the  data  for  elevated  g  conditions  as  well  as  lg  conditions. 

In  well -performed  pul  «e  maneuvers  and  lightly  damped  oscillations,  it  is  possible 
to  determine  a  2-secon-i  to.rioi  to  within  0.02  second.  Good  accuracy  in  damping  can  be 
measured  for  doping  ratios  less  chan  0.2.  The  accuracy  of  period  and  damping  measure¬ 
ments  becomes  rather  pc<  r  for  damping  ratios  greater  than  about  0.3. 
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6.5.2  Constant-Heading  Sideslip  Maneuvers 

In  the  absence  of  pertinent  aid  applicable  pulse -maneuver  data  or  in  an  effort  to 
complement  such  data,  the  constant -heading  sideslip  maneuver  can  he  used  to  determine 
the  weathercock  and  effective  dihedral  derivatives  Cn/3  and  Cj^  ,  provided  control - 
effectiveness  derivatives  are  available  from  other  maneuvers. 

Because  of  frequent  loose  usage  of  terminology,  the  expression  “steady  sideslip”  is 
used  when  “constant -heading  sideslip”  is  meant.  Actually,  a  sideslip  can  be  accom¬ 
plished,  as  shown  in  Figure  43,  as  a  wings-level  sideslip  in  which  yaw  rate  and,  hence, 
a  changing  heading  is  involved,  as  a  constant -hea  ing  sideslip  in  which  a  constant 
linear  flight  path  is  maintained  (r  =  r  =  0),  or  as  a  combination  of  thpjie  two  varia¬ 
tions  of  sideslipping  maneuvers.  The  distinctions  in  the  variation  of  the  sideslip 
maneuver  affect  the  parameters  involved  in  the  analysis  of  the  flight  data  and  the 
format  of  the  equations  employed. 

It  is  difficult  to  perl-  rm  the  sideslip  maneuver  as  a  steadily  increasing  sideslip 
at  a  constant  heading  without  experiencing  angular  rate  and  acceleration  transients. 

A  more  successful  approach  to  the  maneuver  is  to  increase  the  sideslip  in  increments 
in  order  to  damp  out  the  angular  rates  at  each  increment  before  proceeding  to  the 
next  increment.  Although  this  manner  of  accomplishing  the  maneuver  involves  more  time, 
it  is  justified  by  the  refinement  and  resulting  usable  data. 

6.5.3  Pullup  and  Push-Over  Maneuver 

This  maneuver,  or  any  one  of  its  variations,  is  intended  primarily  for  handling- 
qualities  investigations.  However,  the  control -effectiveness  parameter,  Cms  ,  can 
be  mathematically  determined  from  the  initial  phases  of  the  maneuver.  The  maneuver 
is  useful  also  in  determining  the  other  longitudinal  derivatives  by  analog-matching 
techniques. 

6.5.4  Recovery -From-Sideslip  Maneuver 

This  maneuver  has  been  valuable  for  determining  lateral -directional  derivatives  by 
the  analog-matching  technique.  Good  conditioning  is  achieved  by  first  reducing  rudder 
input  to  half  the  value  present  at  the  end  of  a  constant-heading  sideslip  and  then 
releasing  it.  This  maneuver  is  considered  in  more  detail  in  Section  7.8.4. 

6.5.5  Elevated-g  Turn  Maneuver 

The  use  of  this  maneuver  in  derivative  determination  was  discussed  in  Section  6.5.1. 

6.5.6  Roll  Maneuver  (Rudder-Fixed) 

This  maneuver  lends  itself  to  the  determination  of  Cjp  and  Ci8  ,  even  though  it 
is  primarily  a  handling-qualities  maneuver.  In  its  execution,  the  roll  is  initiated 
by  an  abrupt  aileron  step  input.  The  initial  phase  of  the  maneuver,  up  to  maximum  roll 
rate,  is  the  useful  portion  for  derivative  analysis.  The  initial  phase  involves  neg¬ 
ligible  sideslip,  an  essential  factor  in  its  utility  for  derivative  analysis. 
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6.6  General  Consents 

Hie  maneuvers  discussed  constitute  those  commonly  used  in  mathematical  analysis  of 
flight  data  for  derivative  determination  wherein  approximate  expressions  for  deter¬ 
mining  individual  derivatives  or  a  more  comprehensive  technique,  such  as  the  graphical 
time-vector  method,  is  employed.  Many  of  the  approximate  expressions  and  the  time- 
vector  method  are  dependent  upon  control -fixed  free-oscillatlon  data  which  are  not 
usable  when  damping  is  high,  thus  leaving  a  vacuum  for  mathematical  analysis  of  suitable 
data.  Least  squaring  of  the  equations  of  motion  has  not  been  too  successful,  inasmuch 
as  proper  conditioning  of  the  motions  is  difficult  to  establish  and  the  requisite 
accuracy  of  the  recorded  data  appears  to  be  lacking.  In  the  absence  of  suitable  mathe¬ 
matical  techniques,  recourse  is  made  to  analog  matching  of  higher  damped  oscillations 
and  response  to  random  inputs. 

At  times,  it  is  desirable  to  perform  maneuvers  for  power-off  as  well  as  power-on 
conditions  to  investigate  the  influence  of  inflow  effects  of  jet  exhausts  and  possibly 
other  jet-exhaust  effects.  This  may  not  be  operationally  feasible  for  jet  engines. 
Jet-exhaust  effects  of  rocket-engine  aircraft  have  been  studied  by  performing  free- 
oscillation  maneuvers  just  prior  to  and  immediately  following  power  cutoff.  Only 
limited  ranges  of  the  records  were  usable  for  the  power-off  oscillations  because  of 
the  decelerations  and  changes  in  altitude. 

The  analysis  of  data  of  a  complete  flight  program  for  th^  determination  of  stability 
derivatives  can  be  tedious  and  exacting.  The  number  of  computations  necessary  for  an 
effective  analysis  of  the  data  makes  it  apparent  that  systematic  procedures  are  helpful. 
Tabulation  forms,  such  as  shown  in  Table  VIII,  that  include  many  pertinent  flight 
quantities  havt  proved  to  be  helpful. 


7.  ANALYSIS  OF  FLIGHT  DATA 

Of  the  many  methods  proposed  Irr  determination  of  stability  and  control  derivatives, 
only  a  few  are  practical  for  a  relatively  rapid  determination  of  the  derivatives  using 
approximate  equations.  The  limitations  of  these  equations  must  be  known  in  order  to 
avoid  improper  application:.  Of  the  more  comprehensive  techniques  of  analysis  proposed, 
the  graphical  time-vector  method  appears  to  be  the  most  practical  and  provides  reliable 
results  within  the  limits  of  its  applications.  When  analytical  techniques  are  not 
applicable,  analog  matching  of  flight  data  has  proven  to  be  a  practical  technique  for 
determining  derivatives  from  flight  data.  In  the  following  sections,  the  preceding 
techniques  are  discussed  at  some  length.  Comments  on  other  detailed  methods  are  also 
included. 

Inasmuch  as  flight-test  instruments  are  referenced  to  the  body-fixed  axes,  the 
derivatives  are  considered  with  respect  to  these  axes.  Conversion  of  the  derivatives 
from  the  body  to  the  stability  system  of  axes,  if  required,  is  accomplished  by  the 
equations  listed  in  Section  2.2. 

7.1  Fundamentals  of  the  Time-Vector  Approach 

Inasmuch  as  some  of  the  approximate  equations  are  based  on  time-vector  considerations, 
it  is  opportune  to  briefly  discuss  time-vector  properties.  Time-vector  methods  of 
analysis  make  use  of  the  time-invariance  of  the  amplitude  and  phase  relations  between 


52 


the  degrees  of  freedoa  of  an  exponentially  daaped  sinusoidal  oscillating  system  (second- 
order  linear  system)  and  the  differential  and  Integrals  of  the  degrees  of  freedoa  to 
determine  the  values  of  these  amplitude  and  phase  relations,  or  to  determine  the  con¬ 
stants  of  the  system  of  equations. 

Consider  the  daaped,  transient,  sinusoidal,  small -perturbation  oscillation  of  the 
rolling  degree  of  freedom.  This  simple  system  is  described  by 

Ap  +  2^EAp  +  w*A <£'  -  o  .  (113) 

The  solution  to  this  equation  is 

M>'  =  !A<£'le  5  ntcos«ndt  ,  (114) 


■here 


^nd  =  aV^1  “  £2>  • 


Differentiating  with  respect  to  time  t  , 


Ap  =  |A</>'|cu  e  n  £cos(&)nrtt  +  tt)  + /(I  -  £2)  cos[<iJn.t  + 


=  IA<£'|a>ne  iW,|t  cos^fyjyt  +  ^  , 


■here  $d  is  the  damping  angle 


£ 

=  tan’1  - 77  . 

d  v(i  -  C2) 


Similarly 


Ap  =  e  n  cos  (wnrit  +  77  +  2$rt) 


Equatior.3  (114),  (116),  and  (118)  show  that  the  amplitudes  of  these  equations  shrink 
at  the  same  rate  and  the  phase  relationship  between  the  amplitudes  is  time  invariant. 
The  amplitudes  of  the  first  and  second  derivatives  of  A 4>'  are  equal  to  the  amplitude 


of  A<£'  multiplied  by  the  undamped  natural  frequency, 


and  by 


respectively. 


The  phase  of  the  derivatives  is  a  function  of  the  damping  angle,  5>d  ,  which  1^  a 


function  of  the  damping  ratjlo, 
the  displacement  vector  A <£' 


displacement  vector  (180  +  2$d). 


£  .  As  shown  in  Figure  44,  velocity  vector  Ap  leads 
by  (90  +  $d),  and  the  acceleration  vector  Ap  leads  the 


Where  more  than  one  degree  of  freedom  is  involved  in  the  damped,  sinusoidal,  tran¬ 
sient  oscillation  system,  and  the  frequency  is  common  to  all  the  freedoms  involved, 
the  instantaneous  absolute  values  of  the  rotating  vectors  may  be  considered  as  ratios 
(referred  to  as  amplitude  ratios)  and  the  phase  relations  of  the  ratios  established. 
These  ratios  of  the  rotating  vectors  and  their  corresponding  phase  angles  are  time 
invariant.  As  a  result,  the  instantaneous  value  of  any  one  degree  of  freedom  may  be 
readily  determined  if  the  characteristics  of  any  one  of  the  motions  are  known  and 
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the  amplitude  ratio  and  phase  angle  relative  to  the  characteristic  motion  are  known. 
For  example,  if  the  known  characteristic  motion  is 

-  £o>  t 

Ar  =  Ar  e  “cob  (Wndt)  (119) 

and,  if  |A/3| / |Ar I  ,  IAp|/|Arl  ,  .  and  $r>r  are  known,  then 

|Ar!e  cos  (&>ndt  +  $gr)  (120) 

lArle  C  1  cos  (&>ndt  +  $pr)  .  (121) 


The  time  invariance  of  the  amplitude  ratios  and  their  phase  angles  permits  the  re¬ 
presentation  of  any  one  of  the  linearized  equations  of  motion  by  vectors.  Fbr  example, 
by  substituting  Equations  (119),  (120),  and  (121)  and  tl'd  differentials  of  Equations 
(119)  and  (121)  into  the  linearized,  small -perturbation,  rolling-moment  equation,  the 
following  format  is  obtained,  using. the  Ar  vector  as  the  reference  for  the  amplitude 
ratios  and  phase  angles 


I,  |AP>  Jxz  IA*'  IapI  b 

— i - Z$-  — — - -Ci  ■  —  —  Z$ 

qSb  lArl  «>r  qSb  I  Ar  I  rr  |Ar|  2V  >r 


(Clr-Cig) 


lArl  _b_ 

lArl  2V  " 


0  ,  (122) 


where 

|Ap|  ^  |Ap|  |Af  I 

lArl  =  ^  lArl  *  |Ar|  =  "n  ’ 


|Arl 

|Arl 


1  ,  and  $rr  =  0  .  (123) 


The  vector  properties  described,  plus  the  requirement  that  the  vector  polygon  re¬ 
presenting  any  one  equation  must  close,  make  possible  the  determination  of  two  unknown 
derivatives  in  any  one  equation.  The  accuracy  with  which  the  unknown  derivatives  are 
determined  is  dependent  not  only  on  the  accuracy  of  the  amplitude  ratios  used  but  also 
on  the  accuracy  of  the  phase  angle  and  the  sensitivity  of  the  unknown  derivative  to 
small  errors  in  the  phase  angles. 

It  should  be  noted  that  the  introduction  of  cross-coup? ing  terms  into  the  equations 
of  motion  would  result  in  nonlinear  equations  and,  hence,  time-variant  relations  of 
the  cross-coupling  terms  relative  +o  the  other  terms. 

7.2  Basic  Flight  Data 

Application  of  many  of  the  simpler  equations  for  determining  derivatives  requires 
an  evaluation  of  the  period  and  damping;  whereas,  application  of  the  time-vector  method 
requires,  in  addition,  the  determination  of  amplitude  and  phase  relationships.  These 
quantities  are  obtained  from  the  free-oscillatlon  portion  of  the  pulse  maneuver,  as 
illustrated  in  Figure  45.  The  spacing  of  the  peaks  of  the  oscillatory  motions  deter¬ 
mines  the  damped  natural  period,  and  a  comparison  of  these  peaks  for  the  different 
oscillatory  quantities  determines  their  phase  rclrtionship.  Determination  of  the 
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phase  relationships  by  an  averaging  process,  typified  by  the  table  in  Figure  45,  has 
provided  more  consistent  data  than  obtained  by  single  readings.  Hie  first  line  of  the 
exaaple  table  lists  the  tioe  of  occurrence  of  consecutive  plus  and  minus  peaks  of  the 
roll  rate  Ap  .  Similarly,  the  second  line  lists  the  plus  and  minus  peaks  of  the  yaw 
rate  Ar  .  The  third  line  lists  the  time  difference  of  the  first  two  lines  in  each 
column.  Since  the  yaw  rate  Ar  is  the  reference  in  this  instance,  the  signs  in  the 
third  line  indicate  that  the  roll  rate  Ap  lags  the  yaw  rate  Ar  .  The  values  in  the 
third  line  are  ave:  aged  and  converted  to  degrees. 

It  will  be  noticed  in  Figure  45  that  a  yawing  divergence  is  evident  in  the  yaw-rate 
record.  To  isolate  the  oscillatory  motions  and  determine  the  time  to  damp  the  oscilla¬ 
tions,  exponential  curves  are  drawn  as  shown.  A  semilog  plot  of  the  double  amplitudes 
included  between  the  exponential  outlines  of  each  motion  versus  time  establishes  the 
time  to  damp  of  the  oscillations  (Fig. 46).  A  comparison  of  the  plotted  double  ampli¬ 
tudes  of  the  variables  determines  the  amplitude  ratios. 

As  stated  earlier,  accuracy  of  measuring  period  and  damping  becomes  rather  poor  for 
damping  ratios  greater  than  about  0.3.  Generally,  configurations  tested  at  moderate 
and  high  altitudes  and  without  damper  augmentation  have  been  rather  lightly  damped  so 
that  free-oscill&tion  methods  of  analysis  can  be  applied  with  good  accuracy. 

The  damping  ratio  £  ,  damping  angle  $d  ,  and  the  undamped  natural  frequency  co n  , 
are  obtained,  for  both  short -period  and  phugoid  free-oscillations,  from  the  following 
relations: 


l  =  sin 


tan 


-i  Z0-6931^ 

V277T1 /2, 


=  tan 


-l 


'0.693P) 


i2ttT 


1/2 / 


col  =  co , 


nd 


+  =  T 


2rr 


/0.693 

VTl/2 


(124) 


(125) 


(126) 


7.3  Determination  of  a  and  0  From  Free  Oscillations 
in  the  Absence  of  or  Questionable  a  and  0  Data 

7.3.1  Longitudinal  Free  Oscillations 

Should  the  a  records  be  unavailable  or  questionable  in  free-oscillation  longi¬ 
tudinal  data  and  the  pitch-rate  records  available,  !Aal/IAq|  and  $  may  be  obtained 
by  using  time-vector  techniques.  Once  these  quantities  are  determined,  it  is  a  simple 
matter  to  plot  a  as  a  fmction  of  time  or,  of  more  immediate  concern,  to  determine 
lAajjl/IAal  for  use  in  determining  CN(X  . 

The  complete  procedure  for  determining  lAal/IAql  ,  <1^  ,  and  lAaJ/IAd  is  shown 
in  Figure  47.  Hie  procedure  involves  the  application  of  the  following  linearized 
auxiliary  equation  to  correct  the  sensed  normal  acceleration,  anl  ,  to  the  center  of 
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gravity  of  the  aircraft,  as  shown  in  Figure  47(a), 


IAantl 
I  Aq  | 


q 


I  J£i!  a. 

g  I  Aq !  « 


(127) 


and  the  vector  application  of  Equation  (56b)  (Table  V)  in  Figure  47(b)  in  the  format 


|Aqi  i  |Aq |  ™  ,  |A,|  <™ 


0 


(128) 


to  solve  for  $aq  and 


V  |Aa| 

jA a|  _  g  |Aq| 

|Aq|  "  V 
a>_  - 


(129) 


which  now  permits  the  determination 


|A«nl 

I  Aa| 


lAaj  • 
!  Aq| 


(130) 


When  the  vector  quantities  Aani  and  Aa,j  are  approximately  in  phase  and  Aq  is 
approximately  90°  out  of  phase  with  Aan  ,  which  is  usually  the  case,  the  vector 
Equation  (127)  may  be  solved  by  the  simple  algebraic  format 


|AaJ  jAa^l  +xw 

|Aq|  |Aq!  g  1 


(131) 


7.3.2  Lateral-Directional  Free  Oscillations 

Should  the  /3  records  be  unavailable  or  questionable  in  free-oscillation  lateral- 
directional  data,  and  yaw-rate  records  available,  |A/?l/IArl  and  $gp  may  be  obtained 
by  using  a  vector  solution  of  the  following  linearized  auxiliary  equation  to  correct 
the  transverse  accelerometer  record  to  the  center  of  gravity  of  the  aircraft, 


lAa.*  | 

]A^Z$atr 


atji 


e  lArl  rr 


+ 


g  lArl  *>r 


(132) 


and  the  application  of  equation  (59)  in  the  format 


IA/3|  |Ai*|  !  Ap  |  lAi//' 1  , , 

IAH  Z^r  =  2T£H ^rr-^iZTl  ZV-CLsinl ZVr- 

-  CL  cos  6  cos  4>  -  CL  ^ 


lArl  “*&tr 


(133) 
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7. 4  Equations  for  Longitudinal  Control  and 
Stability  Derivatives 

The  nature  of  the  input  and  the  ensuing  free  oscillations  of  the  longitudinal -pulse 
maneuver  permit  the  use  of  relatively  simple  methods  of  analysis  in  determining  longi¬ 
tudinal  control  and  stability  derivatives.  These  methods  give  results  comparable  to 
those  from  the  more  complicated  methods  investigated.  Only  the  simple  methods  are 
discussed  at  this  time  and  only  data  from  these  methods  are  presented.  Unless  other¬ 
wise  stated,  it  is  to  be  assumed  that  stability  augmentation  systems  are  not  operational 
during  the  maneuver  and  that  the  aircraft  behaves  similarly  to  a  rigid  structure,  in 
that  its  behavior  can  be  represented  by  the  linearized  smal 1 -per tur’  ati on  equations. 


7.4.1  Control-Effectiveness  Derivative,  Cm5H 

The  control-effectiveness  derivatives  are  determined  from  the  initial  portion, 
approximately  0.2  second,  ox  a  rapid  pulse  maneuver  (fig.  49).  During  this  part  of 
the  maneuver,  the  airplane  response  is  almost  entirely  pitch  acceleration,  with  the 
result  that  the  pitch  control-effectiveness  derivative  can  be  determined  from 


„  _  I,  Aq 

"*•  "  qSc  AS 


In  similar  fashion,  the  change  in  normal-force  coefficient  due  to  elevator  deflection 
can  be  determined  from 


W  Aa 

CN*e  =  tsas-;- 


(135) 


With  the  preceding  restriction  in  mind,  it  is  desirable,  for  accuracy,  to  read  the 
peak  control  input  and  acceleration  response  with  a  disregard  of  the  phase  lag  between 
the  two,  as  shown  in  Figure  49.  It  has  been  found  that  the  time  difference  in  peak 
values  of  control  input  and  acceleration  response  is  primarily  the  result  of  instrument 
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phase  lag  and,  to  a  lesser  extent,  air-mass  Inertia  effects.  Analysis  by  this  method 
requires  instruments  with  flat  response  characteristics  extending  to  relatively  high 
frequencies  (8  c/s). 

Pulses  applied  at  slower  rates,  and  thus  extending  over  a  longer  time  interval, 
may  require  inclusion  of  damping  and  angle-of-attack  terms  in  the  equation,  especially 
a  .  This  may  necessitate  the  inclusion  of  instrument  phase-lag  corrections  for  q 
and  a  . 


7.4.?  Slope  of  the  Normal-Force -Coefficient  Curve 

From  the  short-period  free-oscillation  data  of  the  airplane  with  the  controls  fixed, 
the  variation  of  the  normal  -i'oice  coefficient  with  angle  of  attack  may  be  evaluated 
from 


C"a 


_W  |Aanl 

qS  | Aocl 


(136) 


This  expression  neglects  the  pitching-velocity  and  angle-of-attack-rate  terms  of  the 
short-pericd  form  of  the  normal-force  equation  (Equation  (58),  Table  V).  These  terms 
have  been  found  to  be  negligible,  as  will  be  noticed  in  the  typical  vector  diagram 
(Fig. 50)  of  the  vector  form  of  this  equation  wherein  the  pitch  rate  was  used  as  the 
base  of  the  amplitude  ratios. 

In  instances  where  “free-oscillation  data"  have  inadvertent  inputs  of  the  elevator 
and  the  angle-of-attack  data  have  been  ascertained  as  reliable,  CNfit  may  be  deter¬ 
mined  hy  selecting  those  portions  of  the  time  history  in  which  the  elevator  is  at  its 
steady-state  position  and  plotting  an  versus  a  for  a  numoer  of  data  points  which 
encompass  the  range  of  an  on  the  records.  The  slope  of  the  plotted  points  is 
lAaJ/IAal  .  This  fundamental  technique,  which  involves  some  labor,  may  still  be  the 
simplest  technique  where  a  control-fixed  free  oscillation  is  heavily  damped  and  thus 
precludes  the  determination  of  lAa^/IAal  by  other  means. 

The  derivative  Cri0L  may  be  converted  to  the  effective  lift-curve  slope,  Ci,a  , 
which  includes  the  contribution  of  power,  by  using  Equation  (38).  The  inclusion  or 
exclusion  of  the  power  term  depends  upon  the  influence  of  power.  For  conventional 
low-performance  aircraft,  CL(X  2:  CNa  at  small  angles-of-attack. 


7.4.3.  The  Derivative  (CNq  +  CN^) 

As  explained  in  Section  3.4,  the  phenomenon  involving  a  is  different  from  that 
involving  q  .  The  pairing  of  the  derivatives  as  (CNq  +  Cn^)  is  valid  only  for  longi¬ 
tudinal  small-perturbation,  free-oscillation  maneuvers.  In  this  maneuver,  Aq  and 
Aa  vectors  are  approximately  in  phase  and  |Aa|/|Aq|  ^  1  ,  thus  permitting  the  pairing. 
Determination  of  the  individual  derivatives  CNq  and  CN&  has  thus  far  defied  solution. 

The  determination  of  (CNq  +  CN^)  itself  is  difficult.  It  may  be  readily  deduced 
from  the  vector  diagram  (Fig. 50)  of  the  following  vector  form  of  the  short-period  mode 
of  Equation  (58), 
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(142) 


V 
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and 


(C«q  +  C«^)  - 


The  approximate  form  of  the  Cga 
(Cno/4 hc)  is  omitted,  results  in  a  small  error  of  the  order  of  3%  or  less. 

Attempts  to  determine  CB(1  and  as  individual  quantities  required  a  precision 

of  flight  data  and  analysis  of  these  data  that  is  difficult  to  achieve.  The  difficulty 
arises  primarily  from  the  acuteness  of  the  phase  angle,  ,  which  is  generally  of 
the  order  of  a  few  degrees;  an  error  of  1°  in  this  phase  angle  can  result  in  large 
errors  in  the  solution. 


7.4.5  The  Phugoid  Static  Stability  and  Damping 
Derivatives  CCu  and  Cnu 

Unlike  the  short -period  mode  of  oscillation  in  which  the  velocity  is  essentially 
non-variant  and  the  angle  of  attack  is  variant,  the  long-period  (phugoid)  node  of 
oscillation  involves  velocity  perturbations  and  essentially  constant  angle -of -attack. 
This  implies  that  any  variations  in  aerodynamic  forces  during  the  phugoid  are  primarily 
the  result  of  perturbations  of  the  normal  and  axial  forces  due  to  the  velocity  per¬ 
turbations,  that  is  to  say  that  C0u  and  Cnu  in  Equations  (58a)  and  (58b)  are  the 
only  derivatives  of  concern. 

Upon  dividing  Equations  (56a)  and  (56b)  by  V  and  substituting  these  equations  for 
Aax  and  Aaj,  in  Equations  (58a)  and  (58b),  respectively,  and  neglecting  second-order 
effects,  the  following  approximate  expressions  are  obtained  for  a  phugoid  initiated 
from  steady-state  horizontal  flight; 


+  CCu  —  ^  -  A  O'  =  0 


mV 


(144) 


and 


-  CNll  —  +  Aq  =  0  . 

u  mV 


(145) 


Hie  characteristic  equation  of  the  phugoid  described  by  these  two  equations  is  a  second- 
order  linear  differential  equation  which  takes  the  Laplace  form 


It  is  readily  recognized  that 


c 

C°u  mV 


^ph^ph 


(146) 


(147) 
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-  BOS 

Cnu  v2  "  ^ph 


Transposing  these  two  equations  results  in  the  following  approximate  equations  for 
determining  Cnu  and  Ccn  from  fligLt  data 


CNn  2 


^ph 


r  ^  2mV^ph6Jnph  _  2m^ph"nph 

V/Cn  *«  _  —  _ 


The  flight  values  of  £pb  and  ^nph  are  determined  from  the  phugoid  oscillations  in 


accordance  with  Equations  (124)  and  (126) 


An  interesting  byproduct  of  this  brief  consideration  of  the  phugoid  parameters 


suggests  itself. 


can  be  considered  to  be  similar  to 


then  Equation  (148) 


takes  on  the  approximate  form 


Thus,  the  phugoid  frequency,  cunph  ,  is  approximately  a  function  of  velocity,  V  ,  only. 

7.4.6  Corrections  for  Effects  of  Stability  Augmentation 

System  in  Determining  Derivatives  from  Short-Period  Oscillations 

In  performing  a  pulse  maneuver  with  the  stability  augmentation  system  engaged,  the 
ensuing  transient  short-period  oscillation  of  the  aircraft  will  be  characterized  by  a 
period  of  oscillation  and  a  damping  ratio  which  will  be  different  from  those  obtained 
with  the  pitch  stability  augmentation  system  off  (Pig. 51).  With  the  system  on,  the 
period  will  decrease  with  increasing  damping  provided  by  the  systun;  whereas,  normally, 
the  period  increases  with  increase  in  inherent  unaugmented  damping.  Ibis  is  due  to  the 
system  gain  and  the  time  constant.  Ihus,  the  gain  and  time  constant  are  factors  to  be 
considered  in  equations  for  determining  the  stability  derivatives,  as  is  brought  out 
in  Reference  42.  The  subsequent  discussion  is  based  on  this  reference. 

Ihe  following  procedure  for  determining  CB(X  and  (CmQ  +  Cm&)  from  flight  data  which 
includes  stability  augmentation  effects  has  been  useful  but  is  of  limited  utility. 

Ihe  principal  value  of  the  ensuing  discussion  is  the  Insight  gained  into  the  complica¬ 
tions  which  may  be  encountered  in  data  which  Include  stability  augmentation  effects. 

For  rigid-aircraft  perturbations  about  a  mean  flight  path,  the  Laplace  transformed 
short-period  mode  two-degree-of-freedom  longitudinal  equations  of  motion  may  be  re¬ 
presented  in  approximate,  but  practical,  form  as 

(s  -Mq)Aq  +  (-MjjS  -  Ma)Aa  =  M5/>Se  (152) 

-Aq  +  (s  -  Za)Aa  =  Z?eASe  . 


(153) 
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In  the  absence  of  pilot  input,  the  transfer  function  for  a  damper  with  a  first-order 
time  lag  may  be  represented  by 


ASe(s) 

Aq(s) 


- -  ~  k ( 1  -  r's)  . 

1  +  t' s 


(154) 


Substituting  Equation  (154)  into  Equations  (152)  and  (153)  results  in  the  following 
determinant 


[(1  +  M8ekr')s  +  (-Mq  -  fisek)]  (-MaS  -  Ma) 
[(Z8ekr')s  +  (-1  -  Zsek)]  (s  -  Za) 


(155) 


whose  characteristic  equation  is 

(1  +  MsekT'  +  MaZsekT')s2  + 

+  [-Za  -  Mq  -  M&  -  fisek  -  Z^kr'  -  (fia  -  Ma)Z8ek]s  + 


+  (-Ma  +  ZaMq  +  Zj.Msek  -  MaZ$ek)  =  0 


(156) 


Considering  only  those  terms  in  Equation  (156)  which  thus  far  have  been  shown  to  be 
significant,  the  short -period  longitudinal  frequency  and  d-"Ding  of  the  aircraft  with 
a  first-order  time-lag  pitch  damper  are 


(w')2  ~ 

“  1  +  M8_kr' 


-M, 


^  -(Za  +  Mq  +  Mft  +  M8ek) 
“  1  +  Msekr' 

Solving  these  equations  for  C,„a  and  (Cnq  +  Cn^), 

c-«  "  ($k*  WT')  K>’ 


(157) 


(158) 


(159) 


((miq  +  ^nia) 


(  -  ^ 

mVc 

(l  +  MqkT'j 

- kCms 

Iy  "5e 

(160) 


Prom  the  above,  it  is  seen  that  CB(X  is  readily  determined  for  a  first-order  linear 
pitch -damper  system.  The  determination  of  (C^  +  CB&),  on  the  other  hand,  may  offer 
a  problem,  inasmuch  as  Casq  in  Mq  is  not  readily  determined  by  itself. 

If  the  pitch-damper  system  is  not  a  first-order  linear  system,  which  is  the  case 
for  many  systems,  analytical  solutions  for  C„a  and  (CBq  +  c®a)  are  impractical.  In 
such  instances,  analog  techniques  are  applied  in  attempts  to  extract  these  derivatives. 
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7. 4 . 7  Representati ve  Resu Its 

Typical  time  histories,  the  flight-determined  period  and  damping  ratios,  and  the 
flight-determined  longitudinal  stability  derivatives  of  the  D-558-II  research  airplane 
have  been  reproduced  in  Figures  52,  £3,  and  54  from  Reference  43.  Most  of  the  data 
were  obtained  from  the  all -rocket-powered  version  of  the  airplane;  the  remainder  of 
the  data  is  based  on  the  jet-  and  rocket -powered  version. 

These  data  have  been  used  to  illustrate  representative  results  because  they  show 
the  need  for  a  concentration  of  flight  test  data  in  the  transonic  zone  to  establish 
the  extent  of  any  abrupt  changes  of  the  derivatives  and  to  show  the  influence  of 
altitude  on  this  particular  aircraft.  Because  the  results  did  not  include  control 
effectiveness,  Figure  55  shows  representative  data  from  Reference  42  for  Cg,$  .  All 
data  shown  were  obtained  from  wings-level  pulse  maneuvers  and  are  typical  of  those 
that  can  be  obtained  from  good  flight  techniques  -  which  include  control  of  flight 
variables,  pilot  skill,  and  instrumentation  -  and  careful  application  of  the  methods 
of  analysis  discussed. 

The  maximum  deviation  from  the  faired  value  in  the  stability  derivatives  shown  in 
Figure  54  is  of  the  order  of  5%  for  Cfja  ,  10%  for  C«a  •  and  20%  for  (CBq  +  C^); 
deviation  of  this  order  of  magnitude  occur  in  only  a  minor  portion  of  the  data  analyzed. 
The  maximum  deviation  of  CBg  in  Figure  55  is  difficult  to  assess  because  the  data 
shown  were  obtained  over  a  lafge  range  of  altitudes  and  elevation  trim  settings; 
however,  the  maximum  deviation  from  faired  values  would  be  of  the  order  of  10%,  which 
would  be  representative. 


7.5  Equations  for  Lateral -Directional  Stability  and 
Control  Derivatives 

The  lateral -directional  control  and  stability  derivatives  are  not  as  readily  and 
reliably  determined  by  the  use  of  approximate  equations  as  are  the  longitudinal  de¬ 
rivatives,  because  of  the  more  complex  behavior  of  the  airplane  and  the  larger  number 
of  derivatives  involved.  In  the  following  discussion,  unless  otherwise  stated,  it  is 
again  assumed  that  stability  augmentation  systems  are  not  operational  during  the 
maneuver  and  that  the  aircraft’  s  perturbed  behavior  can  be  represented  by  the  linearized 
perturbation  equations. 


7.5.1  Control-Effectiveness  Derivatives 


The  basic  procedures  for  determining  lateral  and  directional  control  effectiveness 
are  similar  to  those  previously  discussed  for  longitudinal  control  effectiveness.  How¬ 
ever,  the  expressions  for  lateral -directional  control  effectiveness  are  complicated  by 
the  need  to  account  for  the  possible  influence  of  the  inclination  of  the  principal  axis 
as  well  as  the  aerodynamic  terms.  Tests  with  a  conventional  hi  eh -performance  airplane 
utilizing  a  rapid  control  pulBe  or  step  input  showed  that  the  directional  control  de¬ 
rivative,  Cng  ,  could  be  determined  to  good  accuracy  by  considering  only  the  inertia 
term.  For  example. 


C  -  Iz 

Cn*r  -  |_i£ 


Ar  -  -52L 
qSb 


Ap  -  (Cnr  -  Cn^g)  —  Ar  -  Cnp  —  Ap  -  Cn/SA/3 


2V 


J  Sr 


100  = 


98  - 


(161) 


0  + 


2 


0 
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where  the  magnitudes  of  the  individual  tens  are  given  as  percentages  of  the  answer. 
Ibis  slmplifl cation  in  detenining  CDSr  nay  not  be  applicable  to  other  aircraft. 


For  the  roll-control  derivative,  Cjs  ,  consideration  must  be  given  to  the  aero¬ 
dynamic  derivative  tens.  For  example,  using  the  same  high-perfonanco  airplane  and 
a  rapid  aileron  control  input, 


c4  H^-Ap  — *2  Ar  -  C j_  —  Ap  -  c i  —  Ar  -  CjoA/31 

[qSb  5Sb  tp  2V  Ir  2V  J  Sg 


100  = 


73 


-  4 


31 


-  0 


(162) 


The  cross-control  derivatives,  Cns  and  Cis  .can  be  evaluated  by  using  Equations 
(161)  and  (162),  respectively.  Ihe  cross-control  derivatives  are  usually  of  smaller 
magnitude  and  are  therefore  more  difficult  to  detenine.  It  appears  that  all  aero¬ 
dynamic  tens  may  require  consideration,  as  shown  in  the  following  example  of  the 
analysis  for  Cns&  .  The  flight  quantities  were  obtained  from  the  records  as  shown 
in  Figure  56.  The  time  difference  in  the  peaks  of  the  control  input  and  the  accelera¬ 
tions  is  due  to  the  phase  lag  of  the  instruments.  The  acceleration  and  angular-rate 
records  have  essentially  the  correct  phase  relationship  with  respect  to  each  other  in 
this  instance.  The  magnitudes  of  the  Individual  terms  as  percentages  of  the  answer 
are 


^ns. 


■  [» ■ "  -  44  -  «w  -  s ar  -  -  **»]  i 


100  = 


206  -  141 


10 


16 


(163) 


It  will  be  noticed  that  the  produce-of-inertla  term  is  particularly  significant  in 
this  example.  An  error  in  principal-axis  inclination  would  significantly  affect  the 
answer.  FOr  instance,  in  this  example  an  error  of  1/4°  in  the  inclination  of  the 
principal  axis  (3°)  would  result  in  an  error  of  12%  in  CnJa  . 


7.5.2  The  Side-Force  Derivative,  Cy^ 

This  derivative,  which  contributes  to  the  Dutch  roll  mode  of  oscillation  and  is  an 
index  of  the  pilot’ s  ability  to  sense  transverse  accelerations,  can  be  determined  from 
the  equation 


W  iAaJ 

^  -  i  lA#  • 

The  ratio  |Aatl/|A/3|  is  obtained  from  the  control-fixed  transient  oscillations 
resulting  from  a  pulse  maneuver.  If  the  fi  record  is  suspect  or  missing,  the  ratio 
may  be  determined  from  the  at  and  r  records  as  explained  in  Section  7.3.2  and 
Figure  48.  This  indirect  technique  for  obtaining  lAa^.l/IAfll  is  analogous  to  that 
for  obtaining  lAa^/lAal  and  considers  Cyp  .  C,r  .  and  Cy^  as  negligible. 


64 


7.5.3  Du  Directional-Stability  Derivative,  Cap 

The  stetlc  directional -stability  derivative  la  one  of  primary  Importance,  and  good 
accuracy  la  required  In  Its  measurement.  Although  a  number  of  closed-form  equations 
have  been  used,  each  possesses  limitations  which,  If  not  recognized,  can  lead  to 
erroneous  answers.  Several  of  the  equations  are  based  on  various  degrees  of  degradation 
of  the  following  expression,  the  derivation  of  which  was  based  cn  the  solution  of  the 
determinant  of  the  linearized  lateral-directional  small-perturbation  equations  (Equations 
(61a),  (61b),  and  (61c)).  The  expression  Includes  all  but  the  most  negligible  quantities. 


This  equation  shows  that  when  Is  small,  that  Is,  of  the  order  of  0.08  per  radian 
(0.0014  per  deg)  or  less,  the  ordinarily  insignificant  damping  terms  become  Important. 
In  such  Instances,  Cjp  is  particularly  significant. 


When  Cap  Is  of  an  order  higher  than  0.0014,  Equation  (165)  can  be  reduced  to  the 
following  workable  equation 


*/»  -W^  +  aTc^  -  ¥Cl* 


(166) 


This  expression  can  also  be  obtained  by  differentiating  an  approximate  form  of 
Equation  (61a)  to  provide 


=  -Ar  +  aAp+i-cVgA/3 
fflV  r 


and  also  using  Equations  (61a)  and  (61c)  with  the  assumption  that  C{p  ,  Cjr  ,  C \p  , 
and  CDp  ar=  all  equal  to  zero.  Substitutions  result  in  the  following  linear  differ¬ 
ential  equation. 


A^- 


qS  qSb2 

^  ^  ~WZ  (Cnr "  ^ 


A/3  + 


+ 


qSb 

<*  —  °i/3  +  77-  Ci^sb 
*x  IxIz 


A  fi 


0  . 


(167) 


in  which  the  frequency  term  is  identical  to  Equation  (166). 
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The  fact  that  Cn^  is  a  function  of  Cj^  in  Equation  (166)  may  result  in  question¬ 


able  values  of  Cn 


if  C ip  is  estimated  from  wind-tunnel  data  rather  than  flight 
data,  especially  when  flexibility  effects  as  well  as  other  phenomena  may  appreciably 
alter  the  wind-tunnel  values  of  Cj^  . 


An  approximation  of  Equation  (166)  provides  the  following  simple  expression,  which 
is  of  limited  utility: 


C,„  =  -2-  oj2 
r‘3  qSb  n 


(168) 


The  expression  has  been  used  successfully  on  occasions  when  angle -of- attack  and  dihedral 
effects  were  small.  At  low  indicated  airspeeds,  where  these  effects  are  not  small,  the 
discrepancy  can  be  50%  or  more. 

Values  of  CD^  have  alsc  been  obtained  from  constant-heading  sideslip  maneuvers 
using  the  expression 


un/3 


-  “(cnsr^r^  +  Cnj^Sa^) 


(169) 


This  simple  expression  is  obtained  from  Equations  (61b)  and  (61c)  with  the  stipulation 
that  angular  rates  and  accelerations  are  zero  during  the  sideslip  maneuver.  The 
successful  use  of  this  equation  is  dependent  upon  the  accurate  determination  of  the 
apparent  stability  parameters  Sr^  and  as  well  as  the  control -effectiveness 


*0 


derivatives.  The  results  obtained  from  Equation  (169)  have  shown  a  relatively  poor 
consistency  in  the  supersonic  speed  range,  primarily  because  of  the  difficulty  of 
obtaining  sufficient  sideslip  angle  at  supersonic  conditions  to  make  accurate  deter¬ 
mination  of  the  apparent  stability  parameters. 


In  instances  where  the  influence  of  Ixz  and  cnp  is  negligible,  an  accurate 
equation  for  Cn/3  ,  without  the  necessity  of  relying  on  Cip  ,  has  been  derived  from 
the  yawing-moment  equation 


— ”Ar  -  (Cn_  -  CaB)  ~Ar  -  Cn«A/i  =  0 
qSb  r  ■  2V  p 


(170) 


and  the  following  expressions  for  a  transient  oscillatory  sinusoidal  motion: 


.  .  I  Art  ..  .  -C«_t  .  1A  ,  -C»nt  /  v  \ 

Ar  =  t— :  lArle  cos  (£o„dt  +  $rr)  =  o^lArle  cos  [wndt  +  -  +  $d] 
I  Ar  I  \  2  / 


Ar  =  |A: 


■C"„t 


rie 


cos 


(<^n(lt  +  ^rr)  =  lArle  n  008 ^d1 


A/J  =  j^j  lArle'1""* 


COS  (0)ndt  +  $/3r) 


(171) 
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Substituting  expressions  (171)  into  Equation  (170), 
ties,  and  regrouping  results  in 


I  A/3" 

“=**«■*.  cosS  +  cnfl  —  sin  a 
qoo  _  u  •  ur  pr_ 


sin  o;n(]|t  - 


expanding  by  trigonometric  ldenti- 


"  8ln  $d  +  (°nr  ~  <**>  ™ 


M 


1 


The  first  bracketed  quantity  is  a  simulation  of  components  perpendicular  to  the  Ar 
vector;  the  second  is  a  summation  of  components  parallel  to  the  Ar  vector.  Hence 


cos$d 


|  A/31 

+  Cn£  i ah  8ia  9  fir 


0 


(173) 


and 


M 


^ n  sln  $d  +  Mr  ~  W  gf  +  C°fi  C°8  $/9r  =  0 


(174) 


Considering  only  Equation  (173)  at  this  time,  if  the  phase  angle  V  is  of  the  order 
of  90°  and  the  damping  angle  is  small  -  which  are  the  conditions  normally  encountered  - 
then  sin  $£r  and  cos  $d  will  each  be  similar  to  1  and  Equation  (173)  can  be  trans¬ 
posed  to 


_  lz  lArl 

n/3  “  qSb  |A£|Wn 


(185) 


This  equation  provides  accurate  values  of  CD/3  ,  provided  it  is  used  within  the  limita¬ 
tions  imposed  in  its  derivation. 

Table  IX  lists  the  results  of  the  application  of  Equations  (165),  (166),  and  (168) 
to  flight  data  of  the  F-104  and  YF-102.  The  values  of  Cn^  ,  as  determined  by  Equation 
(l'J5),  are  used  as  reference  values.  For  the  P-104,  Equation  (156)  shows  good  corre¬ 
lation  with  the  reference  value  because  of  the  high  value  of  ,  whereas  the  simple 
frequency  equation  (Equation  (168))  shows  poor  agreement.  For  the  YF-102,  which  has  a 
low  value  of  CD/3  for  the  flight  condition  shown,  Equation  (166)  shows  a  significant 
discrepancy  with  reference  CD/3  and  points  up  the  influence  of  the  damping  terms  when 
CD/3  is  small.  For  this  same  case,  it  will  be  observed  that  the  simple  frequency 
equation  is  unworkable. 

A  relative  comparison  of  the  results  obtained  for  the  F-100  airplane  using  Equations 
(166),  (168),  and  (169)  and  the  results  obtained  using  the  more  comprehensive  graphical 
time- vector  method  (to  be  discussed  later)  are  shown  in  Figure  57.  Considering  the 
graphical  time-vector  results  as  most  representative  for  the  airplane,  it  will  be 
observed  that  the  simple  frequency  equation  (Equation  (168))  would  show  poorest  corre¬ 
lation  at  low  subsonic  speeds  due  to  angle-of-attack  and  dihedral  effects  not  accounted 
for  in  the  equation,  whereas  Equation  (169)  shows  poorest  results  in  the  supersonic 
region  because  of  the  difficulty  in  obtaining  accurate  values  of  Sr^  and  SB/3  . 


Table  X  compares  the  values  of  determined  from  analog  matching  of  oscillatory 

maneuvers  of  the  X-1S  airplane  with  values  of  Cn/3  determined  from  Equations  (166) 
and  (175).  The  values  of  Ixz  and  Cnp  are  essentially  equal  to  zero  on  this  vehicle. 
The  agreement  between  analog  values  of  Cn/3  and  the  equations  is  good.  In  Equation 
(166),  the  agreement  is  due  to  the  high  value  of  .  Equation  (175)  would  be  the 
more  desirable  to  use  on  this  airplane  because  it  does  not  depend  on  the  use  of  Cj^ 
for  a  solution. 


7.5.4  The  Effective  Dihedral  Derivative, 

Several  simple  equations  for  C|„  are  available  with  limitations  on  their  utility, 
as  in  the  case  with  most  simplified  equations. 

Values  of  C can  be  obtained  from  the  constant -heading  i,Jdeslip  maneuver  using 
the  expression 


clp  -  ~(clsr^r0  +  cl Sfpip)  .  (176) 

The  derivation  of  this  expression  and  circumstances  limiting  its  accuracy  are  identical 
to  that  brought  cut  for  its  counterpart  (Equation  (169)). 

A  comparison  of  C determined  by  Equation  (176)  and  the  more  comprehensive 
graphical  time-vector  method  is  shown  in  Figure  58  for  the  F-100  airplane.  At  low 
Mach  numbers,  the  results  from  the  sideslip  equation  (Equation  (176))  compare  favorably 
with  the  time-vector  results;  at  high  Mach  numbers,  a  large  discrepancy  exists  between 
the  two  methods.  Even  though  Cj^  is  not  one  of  the  derivatives  determined  most 
accurately  by  the  time-vector  method,  the  vector  method  is  the  most  practical  analytical 
means  available  for  evaluating  this  derivative. 


In  Instances  where  the  influence  of  Izz  is  negligible,  it  is  possible  to  combine 
Equations  (166)  and  (175)  to  obtain 
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The  use  of  this  equation  is  subject  to  the  additional  restriction  that  it  should  not 
be  used  when  Cn^  is  small,  as  was  noted  in  the  discussion  of  Equaticn  (166).  A1  so, 
the  equation  must  be  used  with  caution  when  the  angle-of-attack  is  less  than  about 
3°  or  4°.  When  the  angle-of-attack  is  less,  ( | Ar I /|A/S! ) (l/o>n)  may  approach  1.0  and 
the  error  in  reading  |Ar|/IAfi|  from  the  flight  records  may  result  in  an  error  in 
(|Arl/|A/3|)(l/a>n)  that  may  exceed  the  net  magnitude  of  the  parenthesized  quantity. 

If  the  record  is  the  major  contributor  to  inaccuracy  in  the  amplitude  ratio,  the 
technique  discussed  in  Section  7,3.2  may  be  employed  to  determine  the  ratio  without 
recourse  to  the  actual  /3  record. 


A  final  precaution  regarding  the  use  of  Equation  (177)  is  in  order.  At  very  low 
angles-of-attsck,  the  error  in  the  flight -determined  values  of  a  can  produce  large 
errors  in  the  equation;  also,  as  a  approaches  zero,  the  equation  approaches  an  in¬ 
determinate  form,  inasmuch  as  the  bracketed  quantity  itself  approaches  zero. 
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7.5.5  The  Damping-in-Roll  Derivative,  C ip 

Simple  expressions  for  the  determination  of  C;p  are  dependent  upon  a  roll  maneuver 
initiated  from  win?s-level  flight  by  a  step  input  of  the  ailerons.  The  derivations  of 
the  expressions  impose  the  restrictions  that  yav  due  to  aileron,  Cqs  ,  sideslip  due  to 
the  effective  dihedral,  ,  and  product-of-inertia  effect  are  negligible.  If  these 
highly  restrictive  conditions  are  satisfied,  the  following  relation  can  be  employed 
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In  using  this  equation,  C can  bo  determined  from  the  initial  part  of  the  control 
input  as  .  .scussed  in  Section  7.5.1  snd  Ap1  is  determined  at  some  time  point,  tx  , 
on  the  roll -rate  time  history  where  Ap  is  zero  -  the  region  of  steady-state  roll. 


If  desired,  the  separate  determination  of 


can  be  avoided  by  solving  for 
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and  substituting  this  ratio  into  the  equation 
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resulting  in  the  format 
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In  these  last  two  equations,  the  subscript  2  indicates  that  Ap  and  Ap  were  obtained 
at  a  time  point  2  on  the  roll-rate  time  history,  preferably  at  the  point  of  maximum 
rolling  acceleration. 


Although  the  restrictions  imposed  at  the  beginning  of  this  section  seriously  limit 
the  application  of  these  equations,  the  last  equation  (Equation  (181))  is  interesting 
in  that  it  shows  that  Cjp  can  be  obtained  without  requiring  the  solution  of  . 

7.5.6  The  Effective  Damping- in-Yav  Derivative,  (Cur  -  Cn^) 

It  was  pointed  out  in  Section  3.4.4  that  Cnr  and  CD/g  may  be  combined  as  an 
equivalent  derivative,  (Cj,r  -  Cn^),  only  for  oscillatory  meneuvers.  providing  the 
stability  axis  system  is  being  considered  or  that  the  ansle-of-attack  is  small  if  the 
body  axis  system  is  used.  When  the  body  axis  system  is  employed,  this  is  tantamount 
to  saying  that  when  the  amplitude  ratio  IAyi'l/lA/31  ~  1  ,  at  a  <  3°  or  so,  Cnr  and 
Cn^  may  be  combined  as  an  equivalent  derivative  for  yaw  rntes. 
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The  combined  derivatives  are  frequently  shown  in  the  results  of  analysis  of  oscilla¬ 
tory  motions  relative  to  body  axes,  even  though  this  amplitude-ratio  condition  is 
exceeded,  then  this  is  done,  it  means  that  an  effective  value  of  cnr  has  been  obtained 
which  includes  the  influence  of  Cn^  and  the  results  of  the  analysis  based  on  the  use 
of  the  actual  IAi/i'I/IA/SI  have  produced  an  answer  which  is  equivalent  to  the  net  con¬ 
tribution  of  Ar  and  A/3  to  Acn  in  terms  of  Ar  . 

An  approximate  equation  for  (Cnr  -  C„4)  is  obtained  directly  from  the  damping  term 
of  the  second-order  differential  equation  (Equation  (167)).  Inasmuch  as 
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a  transposition  results  in 
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Considering  the  assumptions  made  in  deriving  Equation  (167),  from  which  Equation  (183) 


was  obtained,  and  the  stipulations  regarding  the  combining  of  cnr  and  C, 
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it  may 


be  stated  that  Equation  (183)  will  provide  better  accuracy  when  |A'/',I/|A/3I  ~  i  and 


as  lAp|/|A/3l  decreases  to  satisfy  the  condition  that 
influence  on  the  equation. 


C [p  and  CDp  have  a  negligible 


An  approximate  equation  for  (CD].  -  CDi)  which  has  been  used  successfully  in  the 
X-15  airplane  flight  test  program  was  derived  from  Equation  (174) 
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This  equation  is  a  summation  of  yawing-moment  components  parallel  to  the  Ar  vector 
during  a  free-osci nation  maneuver  and  is  subject  to  the  restrictions  that  Ixz  and 
CDp  have  a  negligible  influence  on  the  yawing  moment. 


Since  V  generally  varies  only  a  few  degrees  from  90°  for  angles-of -attack  less 
than  about  15°,  and  since  the  damping  angle  is  small,  the  preceding  equation  can  be 
reduced  and  transposed  to 
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Analog  records  of  free-oscillation  maneuvers  of  the  X-15  airplane,  on  which 
and  Ixz  are  essentially  zero,  were  analyzed  for  (Cnr  -  Cnj)  by  using  Equation 
and  (184).  The  results,  presented  in  Tabic  XI,  show  that  the  latter  equation  was  better 
suited  for  determination  of  the  effective  damping-in-yaw  derivative,  for  angles-of- 
attack  up  to  approximately  12°,  than  Equation  (183)  for  this  vehicle. 
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7.5.7  Correlation  for  Effects  of  Stability  Augmentation 

System  in  Determining  Lateral-Directional  Derivatives 
from  Dutch  Roll  Oscillations 

When  lateral  and  directional  stability  augmentation  systems  having  first-order  time 
lags  are  operational  during  a  Dutch  roll  ( free-oscillation)  maneuver,  the  effects  of 
the  augmentation  system  on  the  frequency  and  damping  of  the  oscillations  and  on 
fArl/lAySI  may  be  accounted  for  in  the  Rame  manner  as  was  done  for  the  longitudinal 
mode  of  oscillation  in  Section  7.4.6. 

7.6  The  Graphical  Time-Vector  Technique 

The  graphical  time-vector  method  of  analysis44*47,  the  principles  of  which  were  dis¬ 
cussed  and  applied  in  the  initial  part  of  this  section,  is  the  most  common  manual 
technique  used  (or  determining  the  lateral  and  directional  derivatives.  Successful 
application  is  dependent  upon  availability  of  control -fixed,  Dutch  roll  oscillation 
data  wherein  the  damping  ratio  is  less  than  approximately  0.3  to  permit  definition  of 
the  period  of  oscillations,  the  log  decrement  of  the  damping  of  oscillations,  amplitude 
ratios,  and  phase  angles. 

7.6.1  Advantages 

One  advantage  of  the  method  is  that  the  procedure  is  manual,  and  the  analyst  is 
afforded  a  graphical  presentation  of  various  factors  affecting  the  solution. 

Another  advantage  is  that  it  is  possible  to  obtain  solutions  when  the  /3-vane  records 
are  available,  suspect,  or  when  it  is  desired  to  avoid  applying  corrections  to  these 
records.  Bypassing  the  /3  records  was  discussed  in  Section  7.3.2.  It  was  shown  that 
the  vector  polygon  of  the  transverse-acceleration  equation  is  essential  in  the  solution 
of  the  amplitude  ratio,  |A/3|/|Arl  ,  and  the  phase  angle,  V  .  Both  of  these  quanti¬ 
ties  are  used  in  the  vector  polygons  of  the  rolling-  and  yawing-moment  equations  to 
determine  Cn£  and  C ip  when  the  vector  is  used  as  the  base  for  the  amplitude  ratios 
in  the  equations,  as  in  Figures  59(a)  and  59(b).  The  phase  angle  is  used  in  the 
orientation  of  the  A/3  vector  in  relation  to  the  Ar  vectcr  and  provides  a  more 
accurate  value  of  fgr  than  can  usually  be  obtained  directly  from  flight  records. 

The  amplitude  ratio,  lA/3|/)Arl  ,  is  used  to  extract  Cn^  and  from  the  deter¬ 

mined  values  of  Cn^  |A/3l/IArl  and  C ip  |A/3|/|Ar|  in  Figures  59(a)  and  59(b). 

7.6.2  Disadvantages 

One  disadvantage  is  that  the  development  of  a  definite  technique  is  required  on  the 
part  of  the  analyst  to  minimize  what  would  otherwise  constitute  a  rather  time-consuming 
and  tedious  effort  to  obtain  a  consistent  and  reliable  set  of  results, 

Another  disadvantage  is  that  only  two  of  the  three  derivatives  in  each  of  the  rolling 
and  directional  moment  equations  may  be  determined  by  means  of  the  vector  diagram,  thus 
necessitating  an  estimate  or  a  wind-tunnel  value  of  one  of  the  derivatives  in  each  of 
the  equations.  Since  C„p  and  Cjr  terms  in  the  vector  diagrams  (Figures  59(a)  and 
59(b))  are  the  smallest  vectors,  it  is  customary  to  estimate  these  quantities.  The 
errors  in  the  estimated  values  of  Cnp  will  affect  (Cnr  -  Cn^)  primarily;  the  errors 
in  C;r  will  generally  affect  Cip  primarily,  but  to  a  much  smaller  extent.  For  low 
angles-of-attack,  (Cnr  -  Cn^)  may  be  estimated  by  using  Equations  (183)  or  (184)  within 
the  limits  of  their  applicability. 
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7.6.3  Application  of  the  Graphical  Time-Vector  Technique  to 
the  Determination  of  Cn^  ,  (Cnr  -  Cn^),  C ,  and  C j 

Figures  59(a)  and  59(b)  show  the  application  of  the  graphical  time-vector  technique 
to  the  determination  of  .  (C^  -  C,^).  C ,  and  Cjp  .  The  amplitude  ratio, 
|Ap|/lAr|  ,  and  the  phase  angle,  $pr  ,  were  determined  from  a  semilog  plot  such  as 
that  in  Figure  46.  The  ratio  |A/3|/|Ar|  and  phase  ingle  were  obtained  from  a 

transverse-acceleration  diagram  as  discussed  in  Section  7.3.2.  The  remaining  required 
amplitude  ratios  and  phase  angles  were  determined  as  follows 
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The  derivatives  CDp  and  Cjr  .  which  have  relatively  small  influences  in  this  instance, 
were  obtained  from  wind-tunnel  data.  Assuming  there  is  no  question  of  the  accuracy  of 
the  data,  the  tunnel  data  should  be  based  oi.  oscillatory  tests,  inasmuch  as  the  flight 
data  are  based  on  an  oscillatory  maneuver. 

With  the  various  known  vector  quantities  properly  oriented  in  the  respective  diagrams, 
the  diagrams  were  closed  and  the  unknown  vectors  determined  by  drawing  the  unknown 
vectors  in  their  proper  phase-angle  directions,  $pr  and  $pr  .  The  newly  determined 
vectors,  such  as  -Cn^(|Ay3|/!Arl)  and  (Cnr  -  Cn^)  (b/2V)  ,  were  then  reduced  to  obtain 
cny3  •  (Pnr  "  .  “d  Cip  . 

Figures  60(a)  and  60(b),  from  Reference  43,  show  the  results  of  the  application  of 
the  graphical  time-vector  technique  to  the  rocket-powered  D-558-II  research  airplane. 

An  interesting  aspect  of  the  results  is  the  influence  of  power  on  the  stability 
characteristics  of  this  airplane. 

At  times  there  may  appear  to  be  an  incompatibility  within  wind-tunnel  data  when  the 
data  are  compared  to  flight-determined  derivatives.  It  then  becomes  imperative  to 
resolve  the  discrepancy  within  the  tunnel  data  and  between  the  tunnel  data  and  the 
flight  data.  This  is  illustrated  in  the  following  example  wherein  Cn^  was  relatively 
low. 

Dynamic  model  tests  of  a  relatively  rigid  high-performance  aircraft  at  a  set  Mach 
number  and  a  =  6.6°  showed  that  Cnp  =  0.01  and  (Cnr  -  Cn^)  =  -0.14  .  Tunnel  data 
also  showed  Cn^  to  be  equal  to  0.055  on  the  basis  of  static  tests  and  equal  to  0.0757 
on  the  basis  of  oscillatory  tests.  Flight  data  obtained  from  time  histories  of  con¬ 
vergent  transient  oscillations  of  the  quality  shown  in  Figure  45  indicated  that,  when 
the  wind-tunnel  value  of  Cnp  =  0.01  was  used  in  the  time-vector  solution,  Cn/9  was 
equal  to  0.071  and  (CDr  -  Cn^)  was  equal  to  0.313.  It  was  obvious  that 
(Cnr  ~  =  0.313  was  not  representative  of  the  true  characteristics  of  the  aircraft 

in  the  Dutch  roll  mode,  since  its  positive  value  indicated  an  oscillatory  divergence, 
whereas  flight  data  showed  oscillatory  convergence. 
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A  check  of  the  phase  angle  $  (-104°)  by  several  analysts  showed  agreement  within 

a  few  degrees.  It  was  decided  that  a  reasonable  spread  of  uncertainty  for  the  quality 
of  data  -  corrected  for  instrument  phase  lag  -  would  permit  $pr  to  be  105°  ±5°;  at 
worst,  the  uncertainty  would  be  ±10°.  Accordingly,  solutions  for  Cn/S  and  (C„r  -  Cn^§) 
were  obtained  by  urlng  various  values  of  Cnp  and  $pr  (within  the  spread  of  un¬ 
certainty).  The  results  shown  in  Figure  61  in  the  form  of  a  grid  plot  indicate  the 
sensitivity  of  the  determined  values  of  Cn^  and  (Cnr  -  Cap)  to  Cnp  and  $  and 
the  incompatibility  of  the  wind-tunnel  data.  The  tunnel  data  were  incompatible  even 
when  allowances  were  made  for  unvertainties  in  inertia  characteristics  and  readability 
of  flight  data. 

Use  was  made  of  approximate  Equation  (183)  with  due  consideration  to  the  limitations 
of  the  equation  for  higher  angle-of-attack  conditions  to  aid  in  establishing  the  magni¬ 
tude  of  (C„r  -  Cnjj).  For  the  test  condition  of  an  angle-of-attack  of  6.6°,  the -0.458 
value  of  (Cj,r  -  C,,^)  obtained  by  Equation  (183)  could  be  in  error  to  the  extent  of 
100%  or  so.  Hence,  it  was  estimated  that  the  correct  value  of  (CBr  -  Cn^)  was  closer 
to  -0.20  than  -0.458.  Also,  considering  that  the  state  of  the  art  in  obtaining 
(Cur  ~  Cn^)  from  wind-tunnel  tests  was  more  reliable  than  in  obtaining  cnP  .  the  tunnel 
value  (-0.14)  of  (CDr  -  Cn/j)  was  surmised  to  be  representative  of  the  true  value  of 
this  derivative.  Uncertainties  in  the  inertia  characteristics  required  that  some 
deviation  be  allowed  in  this  value  in  obtaining  (Cnr  -  C,,^)  from  flight  data.  It  was, 
therefore,  concluded  that  the  results  of  the  analysis  should  lie  within  the  shaded 
area  shown  in  Figure  61.  Within  this  area,  the  value  of  Cn/3  (0.054)  compatible  with 
$pr  =  -104°  and  (Cnr  -  C^)  =  -0. 14  was  considered  to  be  a  mean  value  and  wss  used 
as  an  analytical  result.  The  corresponding  value  of  C„p  should  have  been  approximately 
-0.04. 

The  best  accuracy  in  determining  Cnp  and  (Cnr  -  Cnp  is  obtained  when  |Ap|/|Ar| 
is  small,  at  which  time  the  influence  of  Cnp  is  relatively  small.  When  the  roll-to- 
yaw  ratio  is  large,  it  may  be  advantageous  to  estimate  (C„r  -  Cn^)  and  attempt  to  solve 
for  Cnp  .  For  low  angles-of-attack,  (Cnr  -  CDjj)  may  be  estimated  by  using  Equations 
(183)  or  (184)  within  the  limitations  of  their  applicability. 

Hie  best  accuracy  in  determining  C ip  and  Cjp  is  obtained  when  the  roll-to-yaw 
ratio  is  large.  At  this  time,  the  influence  of  C ir  is  relatively  small.  In  either 
case,  the  static  derivatives,  Cn^  and  C ,  are  determined  more  accurately  than  the 
rotary  derivatives,  (Cnr  -  C„^)  and  C[p  . 

It  was  previously  pointed  out  that  the  accuracy  of  analysis  becomes  ratber  poor  for 
damping  ratios  greater  than  0.3.  Although  a  good  approximation  of  the  damping  ratio 
for  heavily  damped  aircraft  may  be  obtained  by  comparing  flight  records  with  records 
of  heavily  damped  motions  -  the  damping  ratio  of  which  is  known  -  it  becomes  difficult 
to  draw  accurately  the  exponential  envelopes  of  the  oscillatory  motions  to  obtain 
reliable  values  of  amplitude  ratios. 

7.7  Other  Analytical  Techniques 

The  preceding  discussions  regarding  determination  of  derivatives  from  flight  data 
have  shown  various  limitations.  The  graphical  time-vector  technique,  although  the 
most  successful,  is  not  usable  for  damping  ratios  in  excess  of  about  0.3,  requires 
control-fixed  transient  oscillation  data,  and  requires  the  assumption  of  some  deriva- 
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tires,  which  may,  at  times,  cause  difficulties  in  solutions.  To  overcome  the  limita¬ 
tions  of  the  preceding  techniques,  a  number  of  methods  have  been  proposed  for  the 
comprehensive  determination  of  derivatives  (References  48-54,  for  example).  Some  have 
been  successful  in  practice;  others  have  not.  In  most  instances,  the  degree  of  sophisti¬ 
cation  involved  in  the  proposals  requires  automatic  data-reduction  equipment  and  the 
time  and  effort  does  not  warrant  their  use  when  analog  equipment  is  available  for 
application  of  analog-matching  techniques.  Several  of  the  methods  are  considered  in 
the  following  sections. 

7.7.1  Least  Squaring  of  the  Equations  of  Motion 

A  logical  and  straightforward  method,  on  the  sophisticated  side,  for  determining 
derivatives  from  flight  data  is  the  application  of  the  least-squaring  technique  to  the 
linearized  equations  of  motion.  Plight  quantities  at  discrete  time  points  are  sub¬ 
stituted  into  the  equations  of  motion.  Many  more  data  points  are  selected  than  the 
number  of  unknowns,  and  a  least-squares  process  is  applied  to  evaluate  the  unknown 
derivatives.  As  logical  and  simple  as  the  approach  may  be,  it  has  not  been  employed 
too  successfully  for  several  reasons,  including:  difficulty  in  properly  conditioning 
the  maneuver,  instrumentation  accuracy,  phase  lags  between  instruments,  insufficient 
amplification  of  recorded  data  to  provide  precise  readability,  noise  in  data  readout, 
and  instrument  alinement. 

One  of  the  more  successful  attempts  to  apply  this  technique  was  reported  in  Reference 
48.  To  excite  all  the  lateral -directional  modes  and  give  measurable  control  inputs 
without  exceeding  the  limits  of  the  linearized  equations  of  motion,  the  following 
control  input  program  was  used: 

"From  trimmed  level  flight,  step  the  rudder  causing  the  airplane  to  yaw  and 
then  roll  due  to  dihedral  effect.  When  the  bank  angle  reaches  approximately  20 
degrees,  apply  a  step  aileron  deflection  such  that  the  airplane  will  roll  toward 
a  level  flight  attitude.  In  order  to  obtain  a  sufficiently  long  record  of  the 
response  to  aileron,  the  airplane  is  allowed  to  roll  to  an  opposite  bank  angle  of 
20  degrees  before  stopping  the  recording  and  initiating  recovery”. 

A  typical  time  history  of  this  maneuver  is  shown  in  Figure  62.  All  instruments  had 
similar  response  characteristics  and  high  recording  sensitivity  which  was  compatible 
with  calibration-sensitivity  spread  and  calibration  spread.  Alinement  of  instruments 
was  within  ±0.3°.  Recorded  data  were  clean.  It  was  found  that  noise  in  the  readout 
data  significantly  affected  the  results.  Twenty  discrete  time  points  used  for  the 
least-squaring  process  were  considered  sufficient. 

The  results,  reproduced  in  Figures  63(a)  and  63(b),  show  the  degree  of  consistency 
obtained  after  the  greater-than-usual  precautions  were  taken  to  provide  conditions 
that  would  be  compatible  with  the  needB  of  the  technique.  The  requirements  for  this 
technique  are  undoubtedly  similar  to  those  necessary  to  make  other  promising  techniques 
workable,  such  as  the  method  of  Reference  49.  This  method  is  also  an  equatlon-of- 
motion  technique  utilizing  the  Fourier  transform,  a  method  function  to  remove  de¬ 
pendence  on  initial  and  end  conditions,  and  a  least -squaring  procedure. 
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7.7.2  Frequency-Response  Method 

Methods  have  been  proposed  (References  50-52  for  example)  to  determine  stability  and 
control  derivatives  by  using  frequency -response  data  obtained  from  flight  tests.  The 
method  of  Reference  52  encompasses  tho  solution  of  all  derivatives  through  &  complex 
procedure.  Other  methods,  such  as  that  of  Reference  51,  provide  only  limited  results 
based  on  various  degrees  of  approximation. 

The  method  of  Reference  52  replaces  the  time  plane  with  the  frequency  plane.  Amplitude 
ratios  and  phase  relationships  of  airplane  response  to  control  input  from  frequency- 
response  analysis  of  a  pulse  maneuver51* 55  provide  real  and  imaginary  quantities.  Ihe 
complex  quantities  at  discrete  frequencies  are  substituted  into  least-squared  equations 
solving  for  the  derivatives  desired.  The  method  is  simple  in  theory;  however,  con¬ 
siderable  core,  work,  and  time  are  involved  in  the  application,  and  some  experience 
is  necessary  in  the  selection  of  discrete  frequencies.  These  factors  minimize  Interest 
in  further  studies  of  the  method,  especially  where  time  is  of  the  essence  in  obtaining 
a  relatively  quick  look  at  the  flight  values.  Automatic  data-reduction  equipment  would 
greatly  expedite  the  frequency-response  analysis  and  would  be  useful  for  the  other 
computations  required. 

7.8  Analog-Matching  Techniques 

When  flight  data  are  of  such  a  nature  as  to  preclude  the  successful  use  of  the 
graphical  time-vector  technique  or  the  approximate  equations,  and  when  time  and  expense 
will  not  permit  the  use  of  an  experimentation  with  more  sophisticated  techniques, 
recourse  is  usually  made  to  the  analog  to  determine  the  derivatives  that  will  provide 
the  best  match  of  the  analog  time  history  with  the  flight  time  history  of  a  maneuver. 

The  use  of  the  analog  should  be  considered  as  a  last  resort,  to  be  used  only  when 
other  techniques  cannot  be  applied.  It  is  not  a  “cure-all”,  for  it  can  produce 
erroneous  answers  under  certain  conditions  and  still  provide  a  good  match  with  the 
flight  time  histoiy  of  a  maneuver. 

7.8.1  Conventional  and  High-Speed  Repetitive  Operation  (REPOP) 

Analog  Matching 

The  mathematical  model  of  the  aircraft  for  the  analog  computer  is  provided  by  the 
airplane  equations  of  motion;  when  attitude  records  (such  as  and  i>)  are  available 
and  used  in  the  matching  process,  transformation  equations  are  included  to  transform 
aircraft  angular  rates  about  the  body  axes  to  angular  rates  about  Euler  axes. 

Generally,  the  simplest  mathematical  model  compatible  with  the  needs  of  an  investi¬ 
gation  is  used  to  reduce  the  number  of  analog  components  required  and  to  expediate 
solutions,  A  five-degree-of-freedom  mathematical  model,  involving  the  general  equations 
of  motion,  is  employed  when  longitudinal  and  lateral -directional  cross-coupling  effects 
are  factors  in  the  responses  of  the  airplane  during  the  maneuver.  When  such  cross¬ 
coupling  effects  are  not  factors  to  be  contended  with,  the  longitudinal  and  lateral- 
directional  motions  can  be  treated  independently  and  as  two  separate  analog  programs 
using  the  linearized  equations.  Under  such  circumstances,  the  longitudinal  program 
is  treated  as  a  two -degree -of -freedom  case  (with  velocity  a  constant)  unless  phugoid 
is  being  considered,  which  is  not  often;  and  the  lateral -directional  program  is 
treated  as  a  three-degree-of-freedom  case.  Small -perturbation  equations  may  be  used 
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to  advantage  in  such  Instances,  particularly  when  datums  of  angular  rates  and  Euler 
attitudes  may  be  suspect  and  angular  accelerations  have  excessive  noise  or  are  not 
available. 

Initial  estimates  of  stability  and  control  derivatives  to  be  used  in  the  mathematical 
model  are  obtained  from  available  theoretical  and/or  wind-tunnel  values.  If  possible, 
flight-determined  derivatives  obtained  through  the  use  of  the  approximate  equations 
are  employed.  In  the  absence  of  the  preceding,  the  best  estimates  possible  are  made. 
Initial  estimates  are  required  to  establish  reasonable  scaling  factors  for  the  manually 
adjusted  derivative  potentiometers  to  save  operational  time. 

Inasmuch  as  errors  in  initial  conditions  shift  the  amplitude  or  rotate  the  response 
time  history,  provisions  are  made  on  the  analog  to  program  initial  conditions  through 
manually  controlled  potentiometers. 

Flight  test  inputs  in  the  form  of  aileron  and  rudder  deflections  are  reproduced  on 
function  generator  components  of  the  analog  in  as  faithful  a  reproduction  as  possible 
within  the  limits  of  the  function  generators,  which  have  a  finite  number  of  breakpoints. 
When  these  inputs  are  introduced  into  the  mathematical  model,  the  analog  computes  a 
response. 

In  conventional  analog-matching,  the  response  is  recorded  by  a  strip  recorder.  The 
recorded  response  is  then  compared  with  the  actual  flight  time  history,  which  is  re¬ 
produced  on  clear  plastic  to  overlay  on  the  analog  time  history.  A  mismatch  indicates 
the  need  to  modify  the  values  of  the  derivatives,  possibly  change  signs  of  several  of 
them,  and  possibly  modify  the  initial  conditions.  These  changes  are  made  by  using  a 
judicious  trial -and-error  process  until  a  Hatch  is  obtained. 

The  conventional  matching  technique  is  laborious  because  of  the  need  to  manually 
match  a  strip  record  with  the  overlay  every  time  a  programed  condition  is  modified  in 
order  to  study  the  effect  of  the  modification  and  assess  the  next  condition  to  be 
modified.  The  conventional  technique  may  require  from  several  days  to  a  week  to 
obtain  a  match. 

High-speed  repetitive  operation  (REPOP)  matching  differs  from  the  conventional  in 
several  basic  aspects56.  The  strip  recorder  is  replaced  by  an  oscilloscope  and  the 
response  to  Inputs  is  projected  onto  the  scope,  which  has  an  overlay  fastened  to  it. 

The  projected  response  appears  as  a  stationary  time  history  as  a  result  of  an  automatic 
high-speed  recycling  of  the  response  computation.  The  maximum  recycling  speed  for 
fidelity  is  governed  by  the  time  span  of  the  time  history  to  be  matched  and  the  frequency- 
response  characteristics  of  the  function  generator.  Where  a  cycling  rate  of  250  cycles 
per  second  may  provide  fidelity  for  a  3-  or  4-second  time  history,  it  may  cause  serious 
distortions  in  projections  onto  the  scope  if  a  10-second  time  history  is  projected. 

High-speed  repetitive  operation  matching  relieves  the  operator  of  manual  matching 
of  the  time  history,  permits  him  to  make  rapid  modifications  of  derivatives  and  initial 
conditions,  and  allows  him  to  observe  effectively  the  Influence  of  a  modification  on 
the  response.  When  an  optimum  match  is  achieved  on  the  scope,  a  strip  record  is  made 
and  matched  with  an  overlay  to  check  t.:e  fidelity  of  the  scope  match  and  to  retain  a 
record  of  the  resulting  match.  A  REFOP  match  can  normally  be  achieved  in  4  to  6  hours. 


76 


7.8.2  Advantages  of  the  Analog-Hatching  Technime 

The  analog -matching  technique  for  derivative  determination,  in  effect,  accomplishes 
what  sophisticated  analytical  techniques  (see  the  preceding  section)  have  attempted. 

It  enables  the  determination  of  derivatives  under  circumstances  where  approximate 
equations  and  the  graphical  time-vector  technique  fail.  It  does  not  rely  upon  definite 
restrictive  maneuvers,  although  there  are  some  maneuvers  that  cannot  be  solved  for. 

Test  data  showing  inadvertent  inputs  and  subsequent  disturbances  may  be  used. 

7.8.3  Limitations  of  the  Ar.alog-Matching  Technique 

The  success  of  every  technique  discussed  for  determining  derivatives  was  contingent 
upon  the  proper  conditioning  of  the  maneuvers  involved.  This  is  no  less  true  of  the 
analog-matching  technique.  A  Dutch  roll  maneuver,  induced  fay  a  control  pulse,  in  which 
no  spiral-  or  roll -subsidence  modes  are  significantly  evident,  is  generally  impossible 
to  natch  with  a  unique  set  of  derivatives.  It  will  be  found  that  any  number  of  com¬ 
binations  of  derivatives  will  provide  a  match.  A  maneuver  involving  continuous  oscilla¬ 
tion  of  the  control  surfaces,  as  would  be  the  case  of  lateral-directional  oscillatory 
aotlors  vlth  the  lateral -directional  stability  augmentation  system  on,  will  also  be 
very  difficult  to  match  with  a  unique  set  of  derivatives. 

A  properly  conditioned  lateral -directional  maneuver  for  use  on  the  analog  to  permit 
determination  of  a  unique  set  of  derivatives  for  a  match  should  excite  the  roll  and 
spiral  modes  as  well  as  Dutch  roll  oscillations.  The  likelihood  of  obtaining  a  unique 
set  of  derivatives  is  Increased  when  the  maneuver  is  conditioned  to  include  a  rudder 
disturbance  of  a  step-like  nature,  a  transient  oscillation,  and  an  aileron  disturbance  - 
not  necessarily  in  this  order  -  as  was  mentioned  in  Section  7.7.1  and  also  illustrated 
in  a  recovery -from-sldesllp  maneuver  which  is  considered  in  the  next  section. 

Also,  as  was  mentioned  in  Section  7.6.3,  better  accuracy  will  be  achieved  in  the 
major  directional  derivatives  when  the  iApj/|Ar!  ratio  of  the  dynamic  characteristics 
is  low  (minimizing  Influence  of  Cnp),  and  in  the  major  lateral  derivatives  when  the 
|Ap|/|Ar|  ratio  is  high  (minimizing  the  influence  of  Cjr).  it  may  be  concluded,  then, 
that  CDp  and  Cjr  are  normally  difficult  to  determine  to  any  respectable  degree  of 
accuracy.  The  possibility  of  determining  Cjr  appears  to  improve  with  increasing 
tendency  of  the  aircraft  to  roll  off  during  a  maneuver. 

7.8. 4  Application  of  the  Analog-Matching  Technique 

Figure  64  shows  the  results  of  an  analog  match  of  a  "recovery-from-sideslip”  maneuver 
at  a  Mach  number  of  1.84  and  an  altitude  of  49,400  ft.  The  match  is  typical  for  this 
aircraft,  which  had  negative  effective  dihedral  and  adverse  aileron  yaw  for  the  match 
shown.  Rigid  wind-tunnel  data  corrected  for  flexibility  effects  on  the  actual  vehicle 
predicted  practically  zero  effective  dihedral  and  proverse  yaw  due  to  aileron.  It  was 
impossible  to  substantiate  the  predicted  values  on  the  analog,  and  only  one  combination 
of  derivatives  would  provide  the  match. 

The  following  procedure  is  typical  of  that  employed  in  arriving  at  the  analog  match 
of  the  flight  data  which  did  net  include  rolling  and  yawing  accelerations: 
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(1)  The  mathematical  model  was  represented  by  three  lateral-directional  small- 
perturbation  equations 

=  ?  sin  (<p0  +  Aqb)  -  Ar  +  o^Ap  +  ^  (0,^/3  +  C,s  A8r  +  (V*  ASa)  -  ^  sin<£0 


(11)  In  addition,  the  following  transformation  was  employed  to  determine  the  change 
in  Euler  roll  angle,  which  attained  magnitudes  of  the  order  of  20°  on  occasions 
in  the  maneuvers  under  consideration 

=  Ap  -  (r0  +  Ar)0„cos  (<£0  +  A£)  -  r„0o  cos<£0 

(ill)  Finally,  the  outputs  of  the  mathematical  model  were  applied  to  the  following 
two  equations  to  modify  analog  values  of  A/3  and  Aat  to  correspond  to  the 
indicated  values  of  the  flight  data: 

A/3t  =  A/3  +  ^Ar--^Ap 


Aatj  t  -  sin  (4>0  +  A<£)  +  sin<£0  +  -  (A/3  +  Ar  -  a^p)  - 


-  yinstr 


/  ~  o 

/-8rAr*'  +  2pAp\ 


Ar 

8 


+  xlnstr  „  _  zinstr 


Ap 

8 


(iv)  Starting  with  the  arbitrarily  selected  time  zero  (as  in  Figure  64)  for  the 
time  history  to  be  matched: 


(a)  Cn£  was  adjusted  for  approximate  frequency  match. 

(b)  Die  control  derivatives  were  adjusted  to  provide  an  initial  rough  match 
in  the  magnitudes  of  r  and  4>  . 

(c)  Potentiometers  for  r0  and  /3„  were  adjusted  to  roughly  aline  r  and 
/3  traces  of  the  analog  with  flight  data;  similarly,  potentiometers  for 
4>0  and  p0  were  adjusted  to  roughly  aline  0  and  p  traces.  These 
actions  involved  the  following  analog  integration 


Ar  =  /(r0  +Af)dt 
Ap  =  /(P0  +  Ap)dt 


A/3  =  J(/30  +  A/?)dt  , 
A 4>  -  J(4>0  +  Ad>)dt  . 
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(d)  Since  adjustment  of  <jf>0  and  p0  modifies  allneaent  of  analog  and  flight 
tracer:  of  r  and  /3  .  step  (c)  was  reiterated  as  many  tiaes  as  necessary 
to  obtain  a  rough  allneaent  of  analog  and  flight  traces  of  r  ,  fi  ,  4> , 
and  p  . 

(e)  Attention  was  then  focused  on  the  r  trace  to  obtain  a  more  refined  natch 

of  this  trace  by  more  cautiously  adjusting  CDr  ,  Cq^  ,  Cnj  ,  and  CqS  . 
This  operation  necessitated  adjustment  of  Cj^  ,  ,  Cjp  ,  and  C,r  rto 

keep  the  p  trace  in  line. 

It  should  be  noted  that  the  preceding  five  steps  (iv)  (a) -(e),  which  constitute  an 
initial  phase  of  operation  to  obtain  an  approximate  match,  in  voire  about  1  hour.  The 
explanation  of  the  procedure  is,  hy  necessity,  brief.  It  will  be  readily  appreciated 
that  the  steps  are  iterative  to  keep  the  frequency  of  disturbance,  the  magnitude  of 
disturbances  of  the  various  traces,  and  allneaent  of  analog  and  flight  tine  histories 
compatible. 

The  second  phase  of  the  analog-matching  process  involved  the  following  operations: 

(f)  With  the  r  trace  roughly  matched,  attention  was  focused  on  the  <f>  and 
p  traces  by  manipulating  the  Cj^  ,  c ip  ,  and  C2r  derivatives  and  the 
lateral -control  derivatives  as  necessary.  During  this  operation,  fine 
adjustments  were  required  and  made  on  the  r  and  /3  traces  (as  per 
step  (iv)  (e)). 

(g)  With  <f>  ,  p  ,  r  ,  and  /8  traces  matched  as  closely  as  possible,  attention 
was  focused  on  the  at  trace.  This  involved  Cjp  ,  C,^  ,  and  Cygr  . 

The  last  phase  of  the  analog -matching  operation  involved  making  fine  adjustments  to 
initial  conditions  (to  compensate  for  probable  errors)  and  fine  adjustments  to  the 
derivatives,  in  essence,  performing  an  iterative  procedure  of  the  preceding  operations. 
The  second  and  last  phase  of  the  analog-matching  process  generally  involved  3  to  4 
tours  and,  at  tiaes,  more. 


7.8,5  Accuracy  of  Results  in  Analog-Matching  of  Flight  Data 

As  mentioned  previously,  the  accuracy  of  the  results  in  analog-matching  of  flight 
data  is  largely  dependent  upon  the  conditioning  of  the  maneuver.  For  the  longitudinal 
derivatives,  results  from  a  pullup-and-relaase  maneuver  of  an  advanced  high-performance 
aircraft  showed  the  following  accuracies  based  on  the  amount  the  derivatives  could  be 
changed  before  a  trend  toward  mismatch  became  evident: 


(1)  For  a  strong  maneuver: 

C«a  10% 

20%  to  30% 

(CNq  +  CNd) 


Cm«  5% 
CmSe  10% 
(CmQ  +  Cm4) 


200%  or  more 


20%  to  30% 
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(2)  For  a  weak  maneuver: 


CN<x 

20% 

C-a 

10% 

cnSe 

100% 

20% 

(Ch,  +  cNd) 

200%  or  more 

(C-q  +  <W 

40% 

The  accuracies  of  the  lateral-directional  derivatives  obtained  from  analog-matching 
of  well -conditioned,  release-from-sideslip  maneuvers  of  the  same  air  craft  are  shown 
in  the  following  tabulation,  along  with  the  factors  which  Influence  the  accuracy: 


Cnp  ~  5% 

Cip  -  5%  to  15% 

Cyp  -  5%  to  20% 

Cnr  -  5%  to  30% 

%  -  5%  to  30% 

CJp  -  5%  to  30% 

Cjr  -  5%  to  50% 

Cn$r  -  5% 

CjSa  -  5%  to  15% 

—  5%  to  30% 

Cl$r  -  10%  to  30% 

Cyjr  -  5%  to  50%  or  more 

Cy*a  -  40%  to  100%  or  more 


True  for  any  rudder  release  involving 
more  than  one  cycle  of  oscillation. 

Depended  upon  oscillatory  characteristics 
of  <f>  and  magnitude  of  /3  ,  after 
rudder  release. 

Depended  on  the  magnitude  of  the  at 
oscillations  and  the  average  slope  of 
/S  from  release  to  steady  value. 

Depended  upon  the  amount  of  transients 
the  aircraft  was  allowed  to  go  through 
before  controls  were  applied  again. 

Depended  upon  the  magnitude  of  the  roll 
rate  during  oscillation.  Higher  roll 
rate  showed  better  accuracy. 

Depended  upon  the  magnitude  of  the  roll 
rate  during  oscillation.  Higher  roll 
rate  showed  better  accuracy. 

Depended  upon  magnitude  and  oscillatory 
characteristics  of  rolloff.  Larger 
rolloff  showed  better  accuracy. 

True  for  any  rapid  rudder  input. 


Depended  upon  the  magnitude  of  the  control 
input. 


These  results  may  be  considered  typical  of  what  may  be  expected  in  analog-matching 
cf  flight  data  obtained  from  properly  conditioned  maneuvers.  The  accuracies  may  well 
be  typical  of  those  that  may  be  expected  when  comprehensive  analytical  techniques  are 
used. 
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8.  APPLICATION  OF  FLIGHT  DERIVATIVES 

If  wind-tunnel  data  and  theory  were  infallible,  it  stands  to  reason  that  there  would 
be  no  need  for  flight  determination  of  derivatives.  However,  such  is  not  the  case. 

As  new  concepts  in  aircraft  were  developed,  either  with  regard  to  physical  geometry 
or  propulsion  systems,  and  as  aircraft  fly  in  new  Mach  and  altitude  regimes,  there  is 
the  need  to  verity  aerodynamic  theory  and  wind-tunnel  data  and  various  influences  of 
aeroelastic  deformations  of  prototype  structures  on  stability  characteristics;  to 
provide  supplementary  information  not  obtained  in  limited  wind-tunnel  studies;  and 
to  uncover  the  source  of  discrepancies  between  predictions  and  actual  flight  behavior. 

The  following  discussions  provide  some  insight  into  several  of  these  areas. 

8.1  Verification  of  Rind -Tunnel  Data  and  Theory 

As  the  Mach  capability  of  the  airplane  increases,  the  technology  in  wind-tunnel 
testing  becomes  more  critical  with  regard  to  model  construction,  support  of  the  model, 
and  interpretation  of  the  tunnel  data.  Whereas  theory  depends  upon  wind-tunnel  data 
for  verification,  or  to  fill  in  gaps  where  theory  fails,  the  wind-tunnel  may  depend 
upon  flight  data,  as  new  regimes  of  flight  unfold,  to  verity  testing  techniques. 

Flight  data  pointed  out  the  need  for  a  greater  concentration  of  test  points  in  the 
transonic  region  to  accurately  define  the  stability  characteristics  in  this  region 
(Fig. 54).  Flight  data  showed  also  that  it  was  not  sufficient  to  use  a  cold  jet  stream 
to  simulate  the  exhaust  of  rocket  engines.  Figure  60  shows  the  effect  of  the  jet 
exhaust  of  the  D-558-II  research  airplane  on  the  lateral -directional  stability  character¬ 
istics  of  the  vehicle  in  the  supersonic  region.  The  destabilizing  influence  of  power 
was  the  result  of  a  pluming  of  the  hot  jet  exhaust  and  consequent  formation  of  a 
lambda  shock  wave  at  the  juncture  of  the  vertical  tail  and  the  fuselage.  During  Dutch 
roll  oscillations,  the  shock  wave  on  the  leeward  side  of  the  vertical  tail  moved 
forward,  while  on  the  windward  side  it  remained  attached  to  the  jet  exit.  This  pheno¬ 
menon  it  not  common;  it  was  the  result  of  overexpansion  of  the  jet  exhaust  and  the 
proximity  of  the  trailing  edge  of  the  vertical  stabilizer  to  the  jet  exhausts. 

Another  illustration  of  discrepancy  between  flight  and  predicted  data  involved 
elevator  setting  for  lg  flight.  A  comparison  of  the  variation  of  predicted  and  flight- 
determined  elevator  settings  with  Mach  number  showed  increasing  discrepancy  with  in¬ 
creasing  Mach  number  for  a  constant  center-of-gravity  position.  In  this  instance, 
involving  aeroelastic  effects,  predictions  showed  reasonably  close  correlation  of 
Coa  with  flight  data;  whereas,  Caa  and  Cm3  showed  a  difference  in  trend  as  well 
as  level.  Preliminary  study  of  the  problem  showed  a  need  to  consider  Cm 0  as  well 
as  Cma  and  Cm3  .  Thus,  the  following  pitching-moment  equation  for  trimmed  un¬ 
accelerated  level  flight,  based  on  Equation  (53b)  (Table  IV),  was  used  and  constituted 
the  major  consideration  in  arriving  at  the  most  likely  causes  for  the  discrepancy 
between  predicted  and  flight  trim  settings  of  the  elevator 

Cje0  +  Cma(a  +  ac.(=o)  +  Cms^e  “  0  .  (186) 

The  angle-of-attack  (a  +  acN=o)  was  replaced  by  its  equivalent 


CN 


«  +  OcN=o 


(187) 
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to  determine  both  predicted  and  flight  values  of  by  the  following  new  fomat 

of  Equation  (186),  which  is  the  slope -intercept  expression  for  solving  ,  or 

c 

C-o  =  ”Cn  -  CwSeSe  (188) 

c"a  e 

Also,  Equation  (186)  was  transposed  to  solve  for 


8,  =  1  ,  (189) 

A  comparison  of  predicted  and  flight  values  of  the  ratios  in  Equation  (189)  showed 
the  values  of  the  ratio  0,^0, ?  to  be  essentially  the  same;  however,  CHo/CH{ 

differed  in  line  with  the  discrepancy  in  Sg  .  Calculation  of  the  static  margin  using 
c»a/cHa  •  "hi0*1  "as  employed  in  determining  Cao  ,  also  showed  a  discrepancy  between 
prediction  and  flight.  In  the  final  analysis,  it  appeared  that  the  major  source  of 
discrepancy  between  predicted  and  flight  longitudinal  trim  elevator  settings  was  due 
primarily  to  the  differences  in  C,0  and  . 

An  illustration  of  a  discrepancy  between  wind-tunnel  and  flight  data  involving 
power  effects  and  aeroelasticity  is  shown  in  Figure  65.  This  instance  concerned  the 
F-100  airplane  (Fig.  16),  which  is  considered  to  be  a  relatively  rigid  aircraft  and 
has  its  air-intake  nozzle  at  the  nose.  As  shown  in  Figure  65,  the  variation  of  the 
wind-tunnel  value  of  CBp  with  Mach  number  has  roughly  the  same  trend  as  the  flight - 
determined  value.  However,  there  is  an  appreciable  difference  in  level  that  is  well 
beyond  the  difference  to  be  expected  due  to  the  values  of  moments  of  inertias;  values 
are  known  to  within  5%  at  best.  The  results  of  an  investigation  to  trace  the  sources 
of  the  discrepancy  showed  appreciable  moment  of  momentum  effects  of  air-intake  flow 
and  aeroelasticity  effects  of  the  vertical  tail.  When  the  basic  rigid  tunnel  data 
were  corrected  for  these  two  factors,  fairly  good  correlation  was  achieved  with  the 
flight  data  (Fig.  65). 

A  technique  in  tracking  down  inconsistencies  in  wind-tunnel  data  involving  Cn>3  , 
Cnp  ,  and  (C„r  -  C^)  was  illustrated  in  Section  7.6.3. 


8.2  Effects  of  Aeroelasticity 

The  effects  of  aeroelastic  deformation  of  the  structural  components  on  the  stability 
and  control  characteristics  of  the  aircraft  are  of  prime  concern,  particularly  in 
large  transport  designs,  as  pointed  out  in  Section  6.3.  The  illustration  of  aero- 
elastic  effects  shown  in  Figure  65  represents  an  intuitive  approach  in  accounting  for 
a  discrepancy  between  wind-tunnel  and  flight  data.  This  approach  presumes  the  basic 
rigid  tunnel  data  to  be  correct.  It  also  presumes  that  aeroelasticity  effects  are 
simple  enough  to  permit  reasonably  reliable  calculation  of  corrections  to  the  data. 


As  aircraft  Increases  in  size  and  slenderness,  and  operate  at  increasing  dynamic 
pressures,  aeroelastic  deformations  of  the  structure  assume  increasing  significance. 
The  influence  of  aeroelastic  deformations  on  the  stability  and  control  characteristics 
is  difficult  to  predict  on  the  basis  of  theory.  The  deformations  of  the  various 
components  of  the  structure  affect  the  shock  patterns  of  the  airflow  which,  in  turn, 
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affect  the  stability  and  control  characteristics  in  a  much  more  complex  manner  than 
the  aeroelastic  deformation  of  one  or  two  surfaces  on  a  relatively  rigid  aircraft. 
Rigid-model  data  may  be  questionable  because  of  the  uncertainties  in  the  true  rigidity 
of  the  model  and  model  supports  and  interference  effects.  Thus,  a  more  positive 
approach  is  required  to  assess  flexibility  effects  to  verify  and  improve  theory  and 
develop  tunnel  techniques. 

A  flight  test  technique  for  determining  aeroelasticity  effects  on  stability  and 
control  characteristics  is  outlined  in  Section  6.3.  Ihe  technique,  as  presented,  is 
somewhat  simplified  in  that  the  lifting  components  of  thrust  is  considered  to  be 
negligible.  This  approximation  simplifies  flight  planning,  monitoring,  and  making 
on-the-spot  changes  in  flight  conditions  of  ff  and  h  for  maneuvers  at  constant  M  , 
approximately  constant  CL  due  to  aerodynamic  lift  alone,  and  constant  center-of- 
gravity.  An  average  of  the  postflight-determined  values  of 


W  -  Tsin0 


for  the  test  points  on  the  “constant  M  ,  CL  ,  and  center-of-gravity  line"  in  Figure 
39  -  such  as  points  1  and  2  -  will  constitute  the  representative  value  of  CL  for 
these  test  points.  The  maximum  deviation  from  actual  CL  is  within  the  experimental 
error  of  the  investigation.  The  stability  and  control  derivatives  of  these  points, 
when  plotted  against  dynamic  pressure,  define  a  curve  which  shows  the  effect  of  aero- 
elasticity  on  the  derivatives.  The  curve  represents  only  one  M  ,  CL  due  to  aero¬ 
dynamic  lift  alone,  and  center  -'f-gravity  condition. 

8.3  Stability  Criteria 

Considerations  of  the  stability  of  an  airplane  include  not  only  its  inherent 
stability,  which  is  its  behavior  without  pilot  inputs  following  an  initial  disturbance, 
but  also  its  behavior  in  response  to  pilot  inputs.  In  general,  the  study  of  the 
stability  of  an  airplane  involves  the  effect  of  derivatives  on  the  increase  or  decrease 
of  the  stability.  It  is  ai  oblective  study.  When  the  stability  of  the  airplane  is 
considered  in  the  light  of  the  degree  of  pilot’ s  acceptance  of  the  airplane,  and  pilot 
ratings  are  introduced,  the  study  becomes  subjective  and  is  referred  to  as  a  handling- 
qualities  study.  As  may  be  :eadily  surmised,  one  study  complements  the  other. 

Any  extensive  discussion  of  handling  qualities,  which  integrates  the  pilot  as  a 
human  servosystem  constituting  a  feedback  loop  in  the  control  system,  is  beyond  the 
scope  of  this  paper.  It  would  involve  the  study  of  human  factors  and  is  affected  by 
the  pilot’ s  technical  background  as  well  as  the  depth  of  piloting  background,  the 
types  of  aircraft  flown,  orientation  and  types  of  displays  in  the  cockpit,  and  general 
cockpit  environment.  The  art  and  science  of  handling -qualities  investigations  is 
covered  extensively  in  the  literature  (References  57-65,  for  example). 

8.3.1  Longiiaiiinal  Short-Period  Oscillation,  o>n 

The  response  of  the  airplane  to  an  elevator  input  or  gust  disturbance  will  normally 
include  a  longitudinal  short-period  oscillation.  An  oscillatory  condition  by  itself 
indicates  a  static  oscillatory  stability.  Positive,  neutral,  or  negative  dynamic 
oscillatory  stability  is  dependent  upon  the  presence  of  positive,  zero,  or  negative 
damping  characteristics,  respectively.  A  study  of  the  longitudinal  characteristics 
involves  both  static  oscillatory  stability  and  damping. 
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The  undamped  natural  frequency  (static  oscillatory  stability)  is  a  measure  of  the 
longitudinal  stiffness  of  the  airplane  -  analogous  to  a  spring-mass  system.  This 
longitudinal  stiffness  is  represented  by 

=  -  <ffa  +  i„za)  ] 


=  ”  +  c*a^Na 


q  \  Sc 


It  will  be  noted  that  for  any  one  mass  distribution  and  configuration  of  the  airplane, 
the  longitudinal  stiffbess  is  a  direct  function  of  Clojj  primarily.  Thus,  the  oscilla¬ 
tory  frequency  of  the  airplane  will  decrease  with  decreasing  Cga  and  decreasing  q  . 

It  should  be  noted  that  when  CB(X  is  zero,  a  degree  of  longitudinal  stiffness 
(static  oscillatory  stability)  will  be  present  as  evidenced  by  the  C*  Cj|a  term  in 
the  equation,  providing  is  negative  -  a  normal  situation.  The  contribution  of 

this  term  to  longitudinal  stiffness  will  increase  with  increase  in  CN(X  ,  decrease  in 
mass-density  parameter,  /xc  ,  and  increase  in  dynamic  pressure,  q  . 

In  maneuvering  flight,  the  pilot  feels  the  effect  of  longitudinal  oscillatory 
stiffness  in  the  stick  force  per  unit  normal  acceleration. 

8.3.2  Longitudinal  Short-Period  Damping 

The  longitudinal  short-period  damping  is  expressed  either  as  the  actual  damping 
coefficient  or  as  a  damping  ratio.  The  damping  coefficient  (ft  lb  sec/rad)  is  dependent 
upon  the  aerodynamic  derivatives  CNa  and  (C,1q  +  CB&),  as  shown  in  the  equation 


2H  =  -  [Za  +  <«,  +  M 


C»<*S  ,c  ,  c  0 

—  <CmQ  +  CBi)  zvi^ 


% PV--  (Cwq  +CBd)pV  — 


A  decrease  in  CN(X  or  the  negative  value  of  (Cmq  +  C^)  will  decrease  the  damping  co¬ 
efficient,  2£o>n  .  it  will  be  noticed  that  the  magnitude  of  the  coefficient  is  also 
dependent  upon  the  mass  density  of  the  air  £  and  airspeed  V  ,  as  well  as  upon  the 
airplane’ s  mass  characteristics  and  configuration. 


The  damping  ratio  £  as  may  readily  be  surmised  from  the  preceding,  is  obtained  from 


t  = 


Zq  +  (Mq  +  M&) 
2^-Wq 
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Thus,  for  any  one  mass  characteristic  and  configuration  of  the  airplane,  the  damping 
ratio  £  is  a  funct'on  of  y/p  ,  CNa  ,  (C>Q  +  C1(i),  and  /(-C^). 

8.3.3  Longitudinal  Short-Period  lead  Ten,  -Za 

The  parameter  -Za  ,  which  is  a  function  of  Ckr  ,  is  a  longitudinal  short -period 
lead  term  which  affects  the  lead  of  the  pitch  rate  q  with  respect  to  the  control 
input  8e  and  angle -of -attack  as  shown  by  the  transfer  functions 


-MSeZa  8  +  — 


s2  +  2£oJns  +  GJ2 


^■5  qSS\  / CNaqS\  /  1 


S2  +  2£w_s  +  w2 


(194a) 


s2  +  2£<u,s  +  <u2 


(194b) 


As  shown  in  Reference  64,  the  time  for  peak  amplitude  of  q  due  to  a  step  input  de¬ 
creases  with  decreasing  -Za  .  If  -Za  becomes  sufficiently  small  in  comparison  to 
<a„  ,  the  response  to  a  step  input  can  be  disconcerting.  It  may  be  characterized  in  a 
tracking  task  by  an  initial  increase  in  pitch  attitude  of  the  airplane  followed  by 
dwell,  possibly  with  the  airplane  aimed  at  the  target;  but,  then,  with  no  further 
control  input,  there  will  be  a  subsequent  increase  in  the  attitude.  This  type  of 
behavior  may  give  the  pilot  the  feeling  that  the  airplane  is  unstable. 

A  low  value  of  -Za  may  cause  the  pilot  to  experience  a  looseness  in  pitch,  pitch - 
rate  overshoot,  lack  of  control  precision,  and  higher  control  forces.  On  the  other 
hand,  a  high  value  of  -Za  may  cause  a  tendency  to  overcontrol,  exceed  normal  g  , 
and,  in  general,  give  the  impression  that  the  control  is  too  sensitive. 

8.3. 4  The  Dutch  Roll  Oscillation,  wn 

The  Dutch  roll  mode  of  oscillation,  represented  by  the  following  equation,  based 
on  an  approximation  of  the  second  equation  in  Equations  (78),  is  a  measure  of  direct¬ 
ional  stiffness 


<4  =  N£-  L£  sin  a  +  (N;  +  L')Y^ 


Insofar  as  derivatives  are  concerned,  Cn;3  and  are  normally  the  only  derivatives 

of  any  consequence  in  defining  the  frequency  of  this  mode  of  oscillation.  Of  these 
two  derivatives,  CDp  is  dominant.  It  should  be  noticed  that  when  the  static  direct¬ 
ional  stability  is  zero  (CD/3  =  0),  there  is  still  some  degree  of  oscillatory  stability, 
providing  the  effective  dihedral  is  positive  (Cj^  =  -  )  and  the  product  of  inertia  is 
negative,  or  vice  versa.  Some  aspects  of  the  controllability  of  the  airplane  when 
Cna  is  ne«r  zero  and  slightly  negative  are  reported  in  Reference  60. 
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8.3.5  Dutch  Roll  Damping  Coefficient,  2 £c>Jn 

The  Dutch  roll  damping  coefficient  represented  by  the  following  equation,  based  on 
an  approximation  of  the  first  equation  in  Equations  (78),  gives  the  measure  of  the 
dynamic  stability  of  the  Dutch  roll  mode 


= 


cn/§)  + 


CJ0 

mV 


qS 


(196) 


This  equation  shows  the  interaction  of  the  more  dominant  derivatives  affecting  the 
damping  ratio.  The  equation  is  -rare  accurate  than  that  shown  as  Equation  (182)  in 
that  it  includes  C jp  . 

8.3.6  Dutch  roll  damping  Ratio,  £ 

On  the  basis  of  Equations  (195)  and  (196),  the  damping  ratio  can  be  approximated  to 
at  least  the  first  degree  of  approximation  by 

r  _  -  %  -  Tj  -  c; 


(197) 


.  An  increase 
£  ,  but  also 
increases  direct¬ 
ional  stiffness  but  also  the  Dutch  roll  damping  ratio,  which  may  be  desirable.  Decreasing 
Ng  increases  the  bank  angle  that  is  induced  by  a  given  amount  of  sideslip  in  the  Dutch 
roll  motion,  a  characteristic  which  could  be  detrimental  to  maneuvering  control  of  the 
airplane.  In  addition,  decreasing  Increases  the  amount  of  Dutch  roll  disturbance 
in  the  roll  mode  response  to  a  step  aileron  input  -  as  reflected  in  the  parameter 
(oj^/ojn)2  to  be  discussed  -  and  can  disturb  and  mislead  the  pilot. 
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The  Dutch  roll  damping  ratio  is  strongly  affected  by  and 
in  the  negative  value  of  not  only  increases  the  damping  ratio, 
improves  the  stability  of  the  spiral  mode.  Increasing  not  only 


8.3.7  Stability  Criteria  for  Aileron-Only  Roll  Control, 

The  roll  parameter,  ,  is  the  roll  numerator  to  Dutch  roll  frequency  ratio  of 

<£/8e  response  function.  It  is  represented  by 


1 


(198) 
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The  parameter  is  a  measure  of  the  amount  by  which  the  Dutch  roll  motion  is  excited  when 
aileron  inputs  (rudder  fixed)  are  made  by  the  pilot.  It  is  particularly  important  in 
the  roll  tracking  task  in  which  the  pilot -airplane  combination  can  exhibit  considerably 
different  lateral -directional  oscillatory  tendencies  than  would  be  exhibited  by  the 
airplane  alone.  It  provides  a  good  index  regarding  the  increase  or  decrease  in  stability 
of  the  airplane  during  the  aileron-alone  roll  tracking  task. 


When  =  1  ,  there  is  no  yaw  due  to  aileron  inputs  and  there  is  little  or  no 

Dutch  roll  motion  in  response  to  aileron  input.  When  <  l  ,  the  pilot-airplane 

combination  in  an  aileron -only  tracking  task  will  exhibit  an  effective  damping  ratio 
in  roll  tracking  tasks  greater  than  the  Dutch  roll  damping  ratio.  When  a>*/a>n  >  1  , 
the  effective  damping  ratio  will  be  less  than  the  Dutch  roll  damping  ratio  and  the  roll 
that  results  from  aileron  input  is  augmented  by  the  roll  due  to  sideslip;  this  can 
cause  stability  problems  in  the  roll  tracking  task,  especially  when  the  Dutch  roll 
damping  ratio  is  small  and  M/I/31  is  large. 

Equation  (198)  shows  significant  interaction  of  stability,  control,  and  inertia 
parameters  affecting  eu^/o^  .  The  interplay  of  Cns  ,  Cip  ,  and  Ixz  is  important, 
inasmuch  as  these  parameters  may  havo  either  plus  or  minus  values.  Normally,  Cns 
and  Cip  are  the  controlling  parameters;  thus,  if  the  effective  dihedral  is  positive 
(Ci^  <  0),  CnS  will  have  to  be  adverse  (C„s  <  0)  tc  assure  <  1  and  a 

stabilizing  action  during  the  roll  tracking  task. 


8.3.8  Dutch  Roll  Stability  Criteria,  MI/I/31 


The  amplitude  ratio  \4> I/I/3I  is  a  characteristic  of  the  Dutch  roll  oscillations 

and  is  thus  independent  of  any  excitations  of  control  inputs.  Its  mathematical  relation 
ship  to  derivatives  is  given  by 
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The  complex  interaction  of  the  derivative  parameters  makes  it  difficult  to  determine 
pilot  sensitivity  to  M I/I/3I  .  However,  if  the  airplane  has  high  directional  stiff¬ 
ness  («n  >  l),  low  MI/I/31  ,  reasonable  £  >  0.1  ,  and  adverse  yaw  due  to  aileron, 
the  pilot  generally  does  not  bother  to  coordinate  turns  by  using  rudder,  Inasmuch  as 
the  lateral -directional  stiffness  keeps  sideslip  small  and  the  low  value  of  !</>!/!/ 3| 
keeps  roll  due  to  sideslip  small  (Ref. 64). 

If  Ml/|/3|  is  large  (of  the  order  of  4  or  more),  rudder  coordination  becomes 
necessary  in  maneuvering  to  keep  sideslip  small  in  order  to  minimize  the  roll  due  to 
sideslip.  If  the  airplane  is  characterized  by  favorable  yaw  due  to  aileron  (Cnj5  >  0) 
as  well  as  high  values  of  Ml/I/3|  ,  the  pilot  uses  a  cross-coordination  of  rudder 
and  aileron  controls  (right  aileron  and  left  rudder)  to  prevent  excessive  rolls  in 
maneuvers  (Ref. 64).  It  is  not  difficult  to  achieve  coordination  of  controls,  pro¬ 
viding  the  airplane  is  not  excited  hy  external  disturbances.  However,  because  this 
cross-coordination  is  unnatural,  the  pilot  is  more  critical  of  favorable  yaw  due  to 
aileron  (CnSa  >  °)  than  adverse  yaw  due  to  aileron  (Cnga  <  0). 


8.3.9  Roll  Subsidence  Root 


,  iAB 


The  roll -subsidence 
shown  in  the  following 


root,  lAB  ,  is  influenced  most  significantly  by  the  parameters 
equation,  which  is  based  on  the  third  equation  of  Equations  (78) 


(200) 


As  shown,  the  roll  subsidence  is  doninated  by  the  damping-ln-roll  derivative,  Cj^  . 


The  roll -subsidence  root  has  a  direct  Influence  on  the  steady-state  roll  rate  in 
response  to  a  specific  aileron  deflection.  When  the  root  is  large,  the  damping  in  roll 
is  high  and  the  pilot  controls  the  bank  angle  by  commanding  and  adjusting  roll  rate. 
Jhen  it  is  small,  the  pilot  controls  bank  angle  by  commanding  and  adjusting  rolling 
acceleration. 


8.3.10  Spiral-Divergence  Root,  lAs 

The  spiral -divergence  root,  1AS  ,  is  affected  primarily  by  the  parameters  shown 
in  the  following  equation,  which  is  based  on  Equation  (83), 
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The  spiral  mode  can  be  convergent,  neutrally  stable,  or  divergent.  Thus,  for  the 
purpose  of  defining  the  spiral  stability  boundary,  the  equation  can  be  shown  as  a 
spiral  stability  criterion 


V*r  “  N/SLr 

>  0 

spirally  convergent 

or,  as  an  approximation, 

=  0 

neutral  spiral  stability 

cl/ 3^nr  "  Cn pC\r  1 

,<0 

spirally  divergent  . 
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It  will  be  noticed  that  spiral  stability  is  dependent  upon  the  interaction  of  four 
derivatives.  Since  Cny8  is  normally  positive  and  cnr  and  C  ijg  are  normally  negative, 
it  is  well  to  have  Cjr  negative.  Under  any  circumstance,  Cnr  should  be  greater 
than  Cn/9C ir  for  spiral  stability. 


A  divergent  spiral  mode  will  result  in  the  airplane  performing  an  increasing  nose- 
down  and  tightening  turn  accompanied  by  an  Increase  in  speed  and  loss  in  altitude. 


8.4  Flight  Guidance 

Research  vehicles  that  incorporate  new  concepts  of  aerodynamic  configuration,  or 
research  vehicles  designed  for  flight  in  previously  unexplored  regions  of  flight  (Mach 
and  altitude),  usually  have  a  considerable  amount  of  wind-tunnel  investigations  per¬ 
formed  on  models  to  check  their  stability  and  control  characteristics.  Despite  the 
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comprehensiveness  of  the  tunnel  tests,  there  will  be  gaps  in  the  data.  In  addition, 
there  is  normally  a  certain  amount  of  reserve  in  placing  complete  confidence  in  the 
data.  As  a  result,  the  flight  envelope  is  built  up  gradually,  using  stability  and 
control  maneuvers  to  obtain  flight -determined  stability  and  control  derivatives  to 
verify  wind-tunnel  data. 

Agreement  in  the  comparisons  results  in  a  more  rapid  buildup  of  the  flight  envelope; 
disagreement  involves  a  slowdown  until  the  flight  data  can  be  reduced  and  cautiously 
extrapolated.  The  most  representative  values  of  the  stability  and  control  character¬ 
istics  are  used  in  stability  criteria  and  are  programed  into  a  flight  simulator,  in 
which  the  pilot  simulates  the  Intended  mission  and  emergency  conditions  to  reduce  the 
amount  of  risk  that  would  otherwise  be  Involved  in  actual  flight.  The  simulator 
normally  uses  the  general  equations  of  motion  for  a  mathematical  model. 

When  roll-coupling  instability  became  a  physical  reality  with  the  loss  of  several 
F-100  airplanes,  considerable  effort  was  expended  at  the  NASA  Flight  Research  Center 
in  flight  and  simulator  studies  of  the  problem66’  67.  Because  of  the  complex  nature 
of  the  motions,  guidance  of  the  flight  program  using  analog  computations  was  desirable. 
In  a  roll  investigation  of  this  type,  a  small  increase  in  aileron  deflection  can  pro¬ 
duce  large  effects  on  airplane  motions.  It  has  been  graphically  demonstrated  on  several 
occasions  that  flight  guidance  based  on  linear  extrapolation  of  flight  data  at  small 
aileron  deflections  can  be  highly  misleading  and  dangerous.  Figure  66  shows  a  repre¬ 
sentative  comparison  of  the  measured  excursions  in  angle-of-attack  and  angle-of-sideslip 
obtained  in  360°  rolls  with  those  predicted  hy  using  flight-determined  derivatives. 

The  good  agreement  has  been  demonstrated  in  most  instances  in  which  flight -determined 
derivatives  have  formed  the  basis  of  calculations.  Consequently,  the  use  of  such 
guidance  in  flight  planning  has  proved  invaluable.  The  use  of  wind-tunnel  and  theoreti¬ 
cal  derivatives  in  analog  studies  has  not  been  as  successful. 


9.  CONCLUDING  REMARKS 

This  paper  has  attempted  to  bring  together  the  various  factors  that  should  be  known 
by  the  engineer  who  is  concerned  with  the  determination  of  stability  and  control  charac¬ 
teristics  from  flight  data  or  the  use  of  these  flight-determined  characteristics  in 
handling-qualities  research. 

The  discussions  have  been  tempered  with  practical  considerations.  The  various 
factors  discussed  and  the  observations  made  are  the  result  of  experience  in  working 
with  flight  data,  developing  techniques,  comparing  the  data  with  predictions,  and 
investigating  the  causes  of  discrepancies, 

The  theoretical  background,  approximations,  and  limitations  of  the  mathematical 
relations  employed  have  been  given  careful  consideration.  The  problems  encountered 
with  several  of  the  more  sophisticated  techniques  have  been  presented  with  the  hope 
that  any  new  comprehensive  technique  that  may  be  proposed  will  take  into  consideration 
some  of  the  practical  problems  with  instrumentation  and  development  of  maneuvers  to 
properly  condition  the  flight  data  for  the  technique. 

The  pulse  maneuver,  properly  executed,  has  been  found  to  be  generally  adequate  in 
exciting  motions  required  for  stability-derivative  analysis  as  well  as  for  determining 
the  characteristics  of  the  oscillatory  modes  if  adequate  instrumentation  and  alinement 
are  provided. 
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For  longitudinal -derivative  analysis,  staple  equations  utilizing  period  and  daaping 
of  the  oscillatory  node  of  the  airplay  were  shown  to  be  as  satisfactory  as  more  com¬ 
prehensive  methods. 

For  lateral -directional  derivative  analysis,  the  graphical  time-vector  method  was 
shown  to  be  the  most  satisfactory  manual  method  of  analysis.  Simple  approximate  methods 
are  useful  if  applied  with  caution. 

Control  effectiveness  can  usually  be  obtained  by  relating  the  peak  acceleration  to 
rapid  control  inputs.  Consideration  must  be  given  to  aerodynamic  contributions  if 
reasonable  accuracy  is  to  be  realized. 

The  analog-matching  technique  for  determining  derivatives  from  flight  data  was 
shown  to  be  a  valuable  method  of  analysis  for  use  in  the  absence  of  data  suitable  for 
analytical  techniques.  However,  the  analog-matching  technique  has  limitations  in  that 
data  must  be  properly  conditioned  in  order  to  obtain  unique  answers.  The  accuracy  of 
the  results  obtained  from  this  technique  and  the  effect  of  the  type  of  maneuver  on  the 
accuracy  may  well  provide  the  clue  to  what  may  be  expected  from  sophisticated  techniques 
that  may  be  proposed. 

The  use  of  flight  data  to  verify  wind-tunnel  results  and  theory  was  discussed  and 
illustrated.  The  possible  inadequacy  of  comparisons  of  flight  data  with  predictions 
for  determining  aeroelastic  effects  was  pointed  out  and  a  flight-planning  technique 
explained  to  permit  determination  of  aeroelastic  effects  from  flight  data  alone. 

Present  instrumentation  and  methods  of  analysis  are  adequate  for  extracting  deriva¬ 
tives  from  flight  data  for  use  in  most  flight-guidance  simulator  studies  and  detection 
of  characteristics  which  have  not  been  predicted  in  the  wind-tunnel. 
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TABLE  I 

Transformation  of  Derivatives  from  Stability  to  Body  Axis 


C„a 

= 

CLacos  a  +  Coasin  a  +  Cc 

Cca 

= 

CDacos  a  -  CLasin  a  -  Cn 

C*a 

= 

(Cma)s 

cn£ 

= 

(CnyS)s  cos  a  +  (Cj^)8  sin  a 

cnr 

= 

(Cnr)8  cos2  a  +  (Cip)8  sin2  a  +  (Cnp 

+  Cir)8  sin  acos  a 

cny§ 

= 

(<^4)30080+  ( C ip) 8  sin  a 

% 

= 

(Cnp)s  cos2  a  -  (CIr)8  sin2  a  -  (Cnr 

-  Cjp)B  sin  acos  a 

Cns 

= 

(cns)8cosa+  (C{j)8  sin  a 

CIy 3 

= 

(Ci^s  cos  a  -  (Cn|S)8sina 

Ch 

= 

(Cir)8  cos2  a  -  (Cnp)8  sin2  a  -  (Cnr 

-  Cip)B  sin  acos  a 

Cl/3 

- 

(Ci^s  cos  a  -  (Cn/3)8  sin  a 

Cb 
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(CIp)8  cos2  a  +  (CDr)8  sin2  a  -  (Cnp 

+  C{r)8  sin  acos  a 

Cis  =  (Cls)8  COS  a  -  (Cns)s  Sin  a 


TABLE  II 

Transformation  of  Derivatives  fro*  Body  to  Stability  Axis 


CLa 
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C|(acos  a  -  Ccasin  a  -  Cd 

C°a 
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Ccaccs  a  +  CNasin  a  +  Cl 

(CaaU 
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«v» 
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(Cn^s 

- 

Cn^g  cos  a.  —  Clp  sin  a 

(Cnp)e 

- 

Qip  cos2  a  -  Czr  sin2  a  +  (Cnr 

-  CZp)  sin  acos  a 

(Cnjis 

= 

Cnj  cos  a  -  Czj  sin  a 

(CZyg)s  =  Cz^gcosa  +  Cn^aina 

(Cir)8  =  Czr  cos2  a  -  Cn^  sin2  a  +  (Cnr  -  Czp)  sin  acos  a 
(Cz^)s  =  Cz^j  cos  a  +  Cn^  sin  a 

(Czp)s  =  Czp  cos2  a  +  Cnr  sin2  a  +  (Czr  +  Cnp)  sin  acos  a 
=  C z 5  cos  a  +  Cnj  sin  a 


(Cz$) s 


Stability  to  Principal 


h0  =  *(*x8  +  *z8)  “  *^z8  “  Ixb)cos2t7+  IXgZg  sin  217 
ho  =  hs 

^zo  =  i^ii  +  Izs)  +  i(*z8  “  Ix8)  cos  2i?  -  Ixszs  Bin  2i? 
Ix0z0  =  0  =  Ix8zs  cos  2i?  -  i(IXB  “  IZg)  sin  2i? 


Principal  to  Body 


^x  =  idio  Izo>  “  *dzo  “  00826 

h  =  ho 

h  =  *dio  +  Iz0>  +  *dZ0  "  I*0)  00826 
Ixz  =  -i(Ix0  -  Iz0)8in2e 
Body  to  Principal 

Ix0  =  tdx  ^  Iz)  “  idz  -  Ix)  cos  2c  ■  Ij2  oln  26 
ho  “  h 

h0  =  *dx  +  I2)  +  tdz  -  Ix)  cos  26  +  Ixz  sin  26 
^XgZg  —  0  =  I xz  cos  26  +  }(IX  —  Iz)  sin  26 


Velocity  Components 

u  =  V  cos /3  cos  a  (50a)  v  =  Vsin/3  (50b)  w  =  V  cos  J3  sin  a  (50c) 

Linear  Acceleration  Equations 


<p  cos  Q  cos  <p  -  Q  sin  <p  (54c)  p  =  (q  sin  #  +  r  cos#)  sec#  (55c) 


Laplace  Transfora  Foraat  of  Snail -Perturbation  Equations  of  Motion 


(g/V)  cos  0  cos  0  (64b) 


1 


TABLE  VII 

Desirable  Characteristics  of  Instruments  for  Free -Oscillation  Maneuver 


Sensitivity 

Undamped 

Damping 

ratio 

Function 

Range 

(per  inch 

natural 

deflection) 

frequency  (c/s) 

a  ,  deg 

±10 

5.0 

8  ov  more 

0.65 

£  .  deg 

±10 

4.0 

8  or  more 

0.65 

q  ,  radian/sec 

±0.2 

0.2 

8  or  more 

0.65 

q  ,  radian/sec2 

±0.5 

0.5 

8  or  more 

0.65 

r  ,  radian/sec 

±0.1 

0.1 

8  or  more 

0.65 

r  ,  radian/sec2 

±0.4 

0.4 

8  or  more 

0.65 

±0.2 

0.2,  rudder  pulses 

8  or  more 

0.65 

p  ,  radian/sec 

±0.6 

0.6,  aileron  pulses 

8  or  more 

0.65 

±0.6 

0.6,  rudder  pulses 

8  or  more 

0.65 

p  ,  radian/sec2 

±6.0 

6.0,  aileron  pulses 

8  or  more 

0.65 

a„  .  g  units 

±1 

1.0 

8  or  more 

0.65 

±0.3 

0.3,  rudder  pulses 

8  or  more 

0.65 

at  ,  g  units 

±0.6 

0.6,  aileron  pulses 

8  or  more 

0.65 
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TABLE  VIII 

i 

i 

i 

!  | 

Fornat  used  by  NASA  Flight  Research  Center  to  Record 
Actual  Conditions  at  Tine  of  Maneuver 

1  | 

:  1 

:  l 

Trace 


a 


f-farwtessf 


iLiiar/rnHB 


Damo.rot*) \Ftt±<  d0  1^* 


i  C.  Cr 


12-0  I  ,*r  | 


Airplane  . 
Configuration 


12  13  14  15  16  IT  18  19 


Period, 

P. 

sec 


-7.5T  2.12  A72  |  7.rtf 
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TABLE  IX 


Airplane 

M 

Altitude, 

a. 

Reference  Cna 
(Eq.(6-53)) 

,  per  radian 

ft 

deg 

Equation  (6~54) 

Equation  (6-56) 

P-104 

0.94 

41,000 

4.9 

0.46 

B 

0.57 

YF-102 

0.74 

40, 000 

6.6 

0.054 

B 

0.106 

TABLE  X 


a/M/h 

15/0.8/60 

3.5/0.8/40 

6.6/1.2/60 

14/1.6/80 

10/2.0/80 

5/2.0/80 

Analog  value 

-0.084 

-0.021 

■ 

BB 

of  C ip 

■HI 

■HU 

B 

B 

Analog  value 

|B 

0.640 

0.445 

0.508 

of  Cn^ 

M 

M 

1 

Equation  (166) 

0.6t>3 

0.434 

0.  498 

Equation  (175) 

mSM 

0.674 

0.451 

0.504 

TABLE  XI 


a/M/h * 

15/0.8/60 

3.5/1.0/40 

18/1.2/80 

6.6/1.2/80 

14/2.0/80 

10/2.0/80 

5/2. 0/80 

Analog  value 
of  (Cnr  -  Cn^g) 

-1.58 

n 

-1.71 

-1.92 

-2.09 

-2.55 

-2.58 

Equation 

(184) 

-2.50 

-1.43 

-2.96 

-1.81 

-4.77 

-2.76 

-2.68 

Equation 

(183) 

-4.22 

-1.647 

-4.69 

-2.36 

-8.46 

-4.41 

-4.39 

*  h  «  altitude/1000 


Pig.  1  Definition  of  body,  stability,  principal,  and  wind  axis  synteas,  control 
deflections,  and  force  and  aoaent  coefficients 


ipo'iol  irons  verse  vertical  plane:  -  SpaHal  lo"«i,udi"al  vertical  plane.. 


(a)  Euler  angle  perturbation  referred  to  the  x,jrrzr  basic  reference  frame 


(b)  Euler  angle  perturbation  referred  to  Xboyb0zlo  axes  servinK 
as  a  secondary  spatial  reference 

t'ig.  4  Several  methods  of  considering  Euler  angle  perturbations 


Effect  of  time  lag  of  modification  of  vortex  flow  about  lifting  surface  on  the 
change  in  Cj,  following  initial  Instant  change  in  a 


Pig. 9  Direct  propulsive  effects  of  propeller 


Fig.  14  Determination  of  pitching  moment  of  inertia 


(b)  Vector  resolution 

Fig.  15  Determination  of  inclination  of  principal  axis  ssd  yawing  moment  of  inertia 


Amplitude  rotio, 


Tangent  of  restoring  spring  ongle,  tan  5 


Pig. 16 


Ph0*°frafh  fh0W4lng  a  general  arrangement  for  determining  inclination  of 
principal  axis  and  yawing  moment  of  inertia.  Springs  attached  to 
mounting  brackets  located  below  wings 


Pig.  17  Amplitude  ratio  lAp|/|Arl 


as  a  function  of  spring  restoring  angle 


View  B-B 

Fig.  18  Details  of  total -pressure  chamber  and  static-pressure  orifices.  Reproduced 

from  Reference  24 


Photograph  of  a  typical  NASA  installation 


boom  and  fuselage  at  low  speeds.  Reproduced  from  Reference  30 


curvature  on  vane  Indications  of  angle-of -attack 


1^*1 


gravity  of  the  aircraft 


; 

i 

\ 

- 

1 

Pig. 29(b)  Detail  of  linkage  and  damping  system 


129 


Fig. 30  Influence  of  interference  angular  velocity  (“q”  rate)  about  spin  reference 

axis  of  a  sensitive  rate  gyro 


correcting  rate  gyro  records  for  instrument  misallnement 


Airborno  components 


\ 
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Pig.  37  Variation  of  the  period  of  an  P-100  series  airplane  as  a  function  of  Mach 
number,  altitude,  angle-of -attack,  and  load  factor  (from  Reference  41) 


90x10* 


Comparison  of  lateral -directional  responses  to  different 


simulation 


l^pical  tine  histories  of  the  longitudinal  response  characteristics  of  the 
test  airplane  resulting  from  abrupt  stabilizer  deflection 


radian/s«c 


rudian/s«c 


P. 

radian/s«c 


0 


Pig. 44  Relation  of  small -perturbation  rolling  velocity  and  acceleration  vectors  to 
small -perturbation  roll-displacement  vector  in  a  transient  oscillation 


Pig. 47  Vector  polution  of  |Aan|/|Aa|  using  pitch  rate  as  base  for  amplitude  ratio 

when  angle-of-attack  records  are  unavailable 


■  1  uslng  yaw  rate  as  base  for  amplitude  ratios 
sideslip  records  are  unavailable 


148 


longitudinal  pulses  performed  on  the  X-15  analog  wit 
tation  system  engaged  and  disengaged  (from  Reference 


Pig. 51  Time  histories  of 
stability  augmen 


Variation  of  lg  lift  coefficient  and  corresponding  trim  angle-of-attack  with 

Mach  number  (from  Reference  43) 


r  2.3.  * 


to: 


2.3,  and  **  C/lin.Jerj  fir—  J 


1,2,3,  and  •* 


Poire  r  off 

Power  on 

V  « 

0 

20,000 

0 

■ 

30,000 

o 

♦ 

4b, tOO 

A 

▲ 

55,000 

tb, 

k 

60.000 

Unflagged  symbols. all- 

rocket  airplay 

Flagged  srymbO- 

rocket 

and  Jet  airplane 

Pig. 54  Variation  of  static  and  dynamic  longitudinal  stability  derivatives  of  the 
D-558-II  airplane  with  Mach  number  (from  Reference  43) 


Pig. 55  Comprrison  of  flight-determined  horizontal -tail  effectiveness  with  wind-tunnel 

results 


Pig.  58  Comparison  of  results  of  determining  by  time-vector  method  and  steady 

sideslip  equations 


||<ar|  rr  fsb\*r)  *fr  t,/>p  *K|  ,«r|  V  I'-zy  * rr' 


(a)  Determination  of  Cn^  and  (Cnr  -  Cn^) 


A  typical  graphical  time-vector  solution  of  yawing  and  rolling  stability 

derivatives  (continued) 


(b)  Determination  of  C ip  and  Cjp 


Pig. 59  A  typical  graphical  time-vector  solution  of  yawing  and  rolling  stability 

derivatives  (concluded) 


M 


M 

(a)  Influence  of  power  on  the  variation  of  lateral  period  and  damping 


Results  of  graphical  time-vector  analysis  of  the  effects  of  power  on  the 
lateral -directional  period,  damping,  and  stability  derivatives  of  the 
D-558-II  research  airplane  (from  Reference  43)  (continued) 


(b)  Influence  of  power  on  the  variation  of  static  and  dynamic  lateral 

stability  derivatives 


Pig, 60  Results  of  graphical  time-vector  analysis  of  the  effects  of  power  on  the 
lateral -directional  period,  damping,  and  stability  derivatives  of  the 
D-558-II  research  airplane  (from  Reference  43)  (concluded) 


p. 

dtg/sec 


Pig.  62 
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Sr,  d*fl  2  \f 


5a,  d«fl  0 


2 

<*•«  1 

0 

/  1  1 

_ »  1 _ 

_ 1 _ 1 _ 

Typical  time  history  of  maneuver  to  determine  derivatives  by  least  squaring 
the  equations  of  motion  (Reference  48) 


Comparison  with  flight  data  of  results  of  analog  simulation  studies  of  360' 
rolls  using  flight -determined  derivatives 


